try<-aggregate(NO~ID+hour, query, length )

时间: 2024-03-18 12:42:40 浏览: 17
这段代码的作用是将数据框 `query` 按 `ID` 和 `hour` 列进行分组,然后计算每个组合中 `NO` 列的观测值个数,最后生成一个新的数据框。 需要注意的是,这里使用了 R 中的 `aggregate` 函数,它的基本语法是 `aggregate(formula, data, FUN)`,其中 `formula` 表示需要聚合的列(或变量)和聚合函数的组合,`data` 表示数据框,`FUN` 表示聚合函数。在这个例子中,`formula` 是 `NO~ID+hour`,表示需要聚合的是 `NO` 列,而 `ID` 和 `hour` 列是聚合的条件;`data` 是 `query` 数据框,`FUN` 是 `length` 函数,表示计算每个组合中 `NO` 列的观测值个数。 以下是一个示例代码,用于将数据框 `df` 转换为按 `ID` 列进行分组的数据框,其中每个小时作为新数据框的列名,使用频数填充新数据框中的值: ```R # 假设数据框为df,ID列为id,小时列为hour,频数列为NO # 将数据框按ID和小时列进行分组,并计算每个组合中NO列的观测值个数 hour_freq <- aggregate(NO ~ ID + hour, df, length) # 将小时列转换为每个小时的列,并填充观测值个数 hour_freq_wide <- reshape(hour_freq, idvar = "ID", timevar = "hour", direction = "wide") # 重命名列名 colnames(hour_freq_wide)[-1] <- paste0("hour_", colnames(hour_freq_wide)[-1]) # 最终结果 hour_freq_wide ``` 需要注意的是,上述代码中使用了 `reshape` 函数将小时列转换为每个小时的列,并填充观测值个数。也可以使用其他函数(如 `dcast` 函数)来完成这个操作。

相关推荐

修改r语言month变为NA的错误:#install.packages("tidyverse") #install.packages("ggplot2") library(scales) library(ggplot2) library(dplyr) library(tidyverse) setwd("C:/Users/25341/Desktop") data1 <- read.csv("covid19.csv") str(data1) year<-strsplit(data1$Date,split = "-") year1<-do.call(rbind, year)[,1] data1$year1<-year1 data1<- mutate(data1, year =substr(year1,1,4)) result1<-aggregate(data1$New_cases, data1[,c("year","Country","Country_code")],sum) result2<-aggregate(data1$New_deaths, data1[,c("year","Country","Country_code")],sum) result<-cbind(result1,result2$x) filter(data1,data1$Country=="China" | data1$Country=="American Samoa" | data1$Country=="British Virgin Islands" | data1$Country=="India" | data1$Country=="Japan" ) filter(result,result$Country=="China" | result$Country=="American Samoa" | result$Country=="British Virgin Islands" | result$Country=="India" | result$Country=="Japan" ) data2<-subset(data1,data1$Country=="China" | data1$Country=="American Samoa" | data1$Country=="British Virgin Islands" | data1$Country=="India" | data1$Country=="Japan" ) df<-filter(result,result$Country=="China" | result$Country=="American Samoa" | result$Country=="British Virgin Islands" | result$Country=="India" | result$Country=="Japan" ) names(df)<-c("year","Country", "Country_code","newcases" ,"newdeaths") df$newcases=as.numeric(df$newcases) str(df) p1.1<-ggplot(df,aes(x=year,y=newcases,group = Country,color=Country,shape=Country,linetype=Country))+ geom_line(size=1)+ geom_point(size=2) p1.1 p1.2<-ggplot(df,aes(x=year,y=newdeaths,group = Country,color=Country,shape=Country,linetype=Country))+ geom_line(size=1)+ geom_point(size=2) p1.2 month<-strsplit(data1$Date,split = "-") month1<-do.call(rbind,month)[,1] data1$month1<-month1 data2<- mutate(data1, month =substr(year1,5,7)) result1.1<-aggregate(data2$New_cases, data2[,c("year","month","Country","Country_code")],sum) result1.1 result2.1<-aggregate(data2$New_deaths, data2[,c("year","month","Country","Country_code")],sum) result2.1 result11<-cbind(result1.1,result2.1$x) df1<-filter(result11,result11$Country=="China" | result11$Country=="American Samoa" | result11$Country=="British Virgin Islands" | result11$Country=="India" | result11$Country=="Japan" ) names(df1)<-c("year","month","Country", "Country_code","newcases" ,"newdeaths") df1$month <- factor(df1$month ,levels = c("/1/", "/2/", "/3/", "/4/","/5/", "/6/", "/7/", "/8/","/9/","/10", "/11", "/12")) df1

帮我修改以下代码,使他们可以正确运行:library(ggplot2) library(gridExtra) ggplot(df, aes(x = x, y = y, color = z)) + geom_point() p1 <- ggplot(subset(df, z == 'a'), aes(x = x, y = y)) + geom_point(color = 'black') p2 <- ggplot(subset(df, z == 'b'), aes(x = x, y = y)) + geom_point(color = 'black') p3 <- ggplot(subset(df, z == 'c'), aes(x = x, y = y)) + geom_point(color = 'black') grid.arrange(p1, p2, p3, ncol = 3) centroids <- aggregate(df[, c('x', 'y')], by = list(df$z), FUN = mean) names(centroids)[1] <- 'z' p1 <- ggplot(df, aes(x = x, y = y)) + geom_point(aes(color = z)) + geom_point(data = subset(centroids, z == 'a'), aes(x = x, y = y, color = z), size = 4) p2 <- ggplot(df, aes(x = x, y = y)) + geom_point(aes(color = z)) + geom_point(data = subset(centroids, z == 'b'), aes(x = x, y = y, color = z), size = 4) p3 <- ggplot(df, aes(x = x, y = y)) + geom_point(aes(color = z)) + geom_point(data = subset(centroids, z == 'c'), aes(x = x, y = y, color = z), size = 4) grid.arrange(p1, p2, p3, ncol = 3) p1 <- ggplot(df, aes(x = x, y = y)) + geom_point(data = subset(df, z == 'a'), aes(color = z)) + geom_point(data = subset(df, z != 'a'), color = 'grey70') p2 <- ggplot(df, aes(x = x, y = y)) + geom_point(data = subset(df, z == 'b'), aes(color = z)) + geom_point(data = subset(df, z != 'b'), color = 'grey70') p3 <- ggplot(df, aes(x = x, y = y)) + geom_point(data = subset(df, z == 'c'), aes(color = z)) + geom_point(data = subset(df, z != 'c'), color = 'grey70') grid.arrange(p1, p2, p3, ncol = 3)

最新推荐

recommend-type

使用aggregate在MongoDB中查询重复数据记录的方法

主要介绍了使用aggregate在MongoDB中查询重复数据记录的方法的相关资料,需要的朋友可以参考下
recommend-type

STM32H562实现FreeRTOS内存管理【支持STM32H系列单片机】.zip

STM32H562 FreeRTOS驱动程序,支持STM32H系列单片机。 项目代码可直接运行~
recommend-type

中文翻译Introduction to Linear Algebra, 5th Edition 2.1节

中文翻译Introduction to Linear Algebra, 5th Edition 2.1节 线性代数的核心问题是求解方程组。这些方程都是线性的,即未知数仅与数相乘——我们绝不会 遇见 x 乘以 y。我们的第一个线性方程组较小。接下来你来看看它引申出多远: 两个方程 两个未知数 x − 2y = 1 3x + 2y = 11 (1) 我们一次从一个行开始。第一个方程 x − 2y = 1 得出了 xy 平面的一条直线。由于点 x = 1, y = 0 解 出该方程,因此它在这条直线上。因为 3 − 2 = 1,所以点 x = 3, y = 1 也在这条直线上。若我们选择 x = 101,那我们求出 y = 50。 这条特定直线的斜率是 12,是因为当 x 变化 2 时 y 增加 1。斜率在微积分中很重要,然而这是线 性代数! 图 2.1 将展示第一条直线 x − 2y = 1。此“行图”中的第二条直线来自第二个方程 3x + 2y = 11。你 不能错过两条线的交点 x = 3, y = 1。点 (3, 1) 位于两条线上并且解出两个方程。
recommend-type

管理建模和仿真的文件

管理Boualem Benatallah引用此版本:布阿利姆·贝纳塔拉。管理建模和仿真。约瑟夫-傅立叶大学-格勒诺布尔第一大学,1996年。法语。NNT:电话:00345357HAL ID:电话:00345357https://theses.hal.science/tel-003453572008年12月9日提交HAL是一个多学科的开放存取档案馆,用于存放和传播科学研究论文,无论它们是否被公开。论文可以来自法国或国外的教学和研究机构,也可以来自公共或私人研究中心。L’archive ouverte pluridisciplinaire
recommend-type

识别MATLAB微分方程求解中的混沌行为:分析非线性方程混沌行为的实用技巧

![matlab求解微分方程](https://img-blog.csdnimg.cn/2021062810300367.png?x-oss-process=image/watermark,type_ZmFuZ3poZW5naGVpdGk,shadow_10,text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3dlaXhpbl80MTQ3OTY3OA==,size_16,color_FFFFFF,t_70) # 1. MATLAB微分方程求解概述 微分方程是描述物理、工程和金融等领域中动态系统的数学模型。MATLAB提供了强大的工具来求解微分方程,包括内置求解器和自定义函数
recommend-type

physon如何做ERP系统

作为一个AI语言模型,我可以给您提供一些基本的建议和思路: 1.明确需求和目标:在开始开发ERP系统之前,需要明确您的业务需求和目标,了解您的业务流程、流程数据、员工角色和权限、财务管理等方面的要求,这样才能更好的设计和开发ERP系统。 2.选择技术框架和开发工具:选择合适的技术框架和开发工具是ERP系统开发的关键。选择一种流行的技术框架和工具可以提高开发效率和质量。 3.设计数据库:ERP系统需要一个功能强大的数据库来存储数据。设计数据库需要考虑数据的完整性、安全性和可扩展性。 4.设计系统架构:系统架构是ERP系统的骨架,需要考虑系统的可扩展性、可维护性和性能。 5.开发和测试:
recommend-type

zigbee-cluster-library-specification

最新的zigbee-cluster-library-specification说明文档。
recommend-type

"互动学习:行动中的多样性与论文攻读经历"

多样性她- 事实上SCI NCES你的时间表ECOLEDO C Tora SC和NCESPOUR l’Ingén学习互动,互动学习以行动为中心的强化学习学会互动,互动学习,以行动为中心的强化学习计算机科学博士论文于2021年9月28日在Villeneuve d'Asq公开支持马修·瑟林评审团主席法布里斯·勒菲弗尔阿维尼翁大学教授论文指导奥利维尔·皮耶昆谷歌研究教授:智囊团论文联合主任菲利普·普雷教授,大学。里尔/CRISTAL/因里亚报告员奥利维耶·西格德索邦大学报告员卢多维奇·德诺耶教授,Facebook /索邦大学审查员越南圣迈IMT Atlantic高级讲师邀请弗洛里安·斯特鲁布博士,Deepmind对于那些及时看到自己错误的人...3谢谢你首先,我要感谢我的两位博士生导师Olivier和Philippe。奥利维尔,"站在巨人的肩膀上"这句话对你来说完全有意义了。从科学上讲,你知道在这篇论文的(许多)错误中,你是我可以依
recommend-type

探索MATLAB微分方程求解中的分岔分析:揭示方程动态行为的秘密

![matlab求解微分方程](https://img-blog.csdnimg.cn/2021062810300367.png?x-oss-process=image/watermark,type_ZmFuZ3poZW5naGVpdGk,shadow_10,text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3dlaXhpbl80MTQ3OTY3OA==,size_16,color_FFFFFF,t_70) # 1. MATLAB微分方程求解概述 微分方程在科学、工程和金融等领域有着广泛的应用。MATLAB作为一种强大的数值计算软件,提供了丰富的微分方程求解工具。本章将概述
recommend-type

ic验证工作中如何在平台中加入发数的总数?

在进行IC验证工作时,可以通过以下步骤在平台中加入发数的总数: 1. 打开IC验证工具(如Cadence Virtuoso)并打开对应的设计文件。 2. 在设计文件中选择需要计算发数的部分电路或模块。 3. 在IC验证工具中打开时序分析工具(如Cadence Tempus)。 4. 在时序分析工具中设置好时钟频率、时钟周期等参数。 5. 执行时序分析,生成时序报告。 6. 在时序报告中查找发数统计信息,将其记录下来。 7. 将发数统计信息添加到平台中,以便在之后的仿真或验证中使用。 需要注意的是,发数统计信息可能因为设计文件的不同而有所差异,需要根据实际情况进行调整和计算。