相似文件
换一批
经管之家送您一份
应届毕业生专属福利!
求职就业群
感谢您参与论坛问题回答
经管之家送您两个论坛币!
+2 论坛币
各位前辈请帮个忙,本人实在是0基础的R,assignment1 已经做了好久了,还是卡死在part 3。
附件是specdata,下面是我写的代码,后面附上测试答案。
先谢谢各位了。
corr <- function(directory, threshold = 0) {
files_list <- dir("~/Desktop/specdata", full.names=TRUE)
data <- data.frame()
for (i in 1:332){
data <- rbind(data,read.csv(files_list))
}
cr<-c()
for (i in id){
data_subset<-subset(data,data$ID==i)
data2<-data_subset[,2:3]
cc<-complete.cases(data2)
if (nrow(cc) > threshold ) {
cr <- c(cr, cor(cc$sulfate, cc$nitrate) )
}
}
print(cr)
}
测试答案:
summary(cr)## Min. 1st Qu. Median Mean 3rd Qu. Max. ## -0.2110 -0.0500 0.0946 0.1250 0.2680 0.7630cr <- corr("specdata", 400)head(cr)## [1] -0.01896 -0.04390 -0.06816 -0.07589 0.76313 -0.15783summary(cr)## Min. 1st Qu. Median Mean 3rd Qu. Max. ## -0.1760 -0.0311 0.1000 0.1400 0.2680 0.7630cr <- corr("specdata", 5000)summary(cr)## Min. 1st Qu. Median Mean 3rd Qu. Max. ## length(cr)## [1] 0cr <- corr("specdata")summary(cr)## Min. 1st Qu. Median Mean 3rd Qu. Max. ## -1.0000 -0.0528 0.1070 0.1370 0.2780 1.0000length(cr)## [1] 323
扫码加我 拉你入群
请注明:姓名-公司-职位
以便审核进群资格,未注明则拒绝
|