请教高手指点R报错，如何定义变量？

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

下面代码在运行的时候说找不到prob2，是不是应该事先定义，但不不知道如何定义？
prob1<-scan()
1 2 50 3.2
1 3 45 3.5
1 4 40 2.8
2 4 50 3.5
2 5 46 3.6
2 7 48 3.2
3 8 51 4.6
3 9 5  5.2
3 10 48 5.9

dim(prob1)<-c(4,9)
prob1<-t(prob1)
colnames(prob1)<-c("id","time","wt","dv")
prob1
prob1<-as.data.frame(prob1)
mindv<-aggregate(x=prob1[,c("dv")],by=list(prob1$id),FUN=min)
mindv

i<-numeric()
for(i in 1:length(mindv$x)){
j<-which(prob1$dv==mindv$x[i])
prob2[i,]<-prob1[j,]
}

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：定义变量 请教高手 高手指点 Rob 不知道 如何

我只会定义excel

prob2需要先定义，因为第一次出现R不知道prob2[i,]是什么，你得先告诉R，prob2是个矩阵，看程序的意思，prob2应该和prob1是相同维度的矩阵，你的prob1是9*4的吧？那就把prob2也定义成9*4的吧。
另外i<-numeric()这句话没有用，可以去掉。

prob2 <- NULL
for(i in 1:length(mindv$x)){
j<-which(prob1$dv==mindv$x[i])
prob2 <- rbind(prob2, prob1[j,])
}

做成外部文件，然后读取

zippo414 发表于 2011-12-15 00:37
prob2需要先定义，因为第一次出现R不知道prob2是什么，你得先告诉R，prob2是个矩阵，看程序的意思，prob2应 ...

首先谢谢你的热心帮助！但我的问题还没用解决。我要的结果是在prob1中提取出与mindv中dv相等的各行数据，即
id time wt  dv
1  1    2 50 3.2
3  1    4 40 2.8
6  2    7 48 3.2
7  3    8 51 4.6
我按你说法定义prob2<-matrix(,nrow=4,ncol=4),但是得不到我想要结果，可能是循环那出了问题，望指教！多谢！

ELIANE! 发表于 2011-12-15 09:48
首先谢谢你的热心帮助！但我的问题还没用解决。我要的结果是在prob1中提取出与mindv中dv相等的各行数据， ...

4楼版主的代码应该能解决问题了吧？

qoiqpwqr 发表于 2011-12-15 04:55
prob2

谢谢，但是我不明白prob2 <- rbind(prob2, prob1[j,])这样就得到了我想要的结果，为什么用rbind？能给我简单说说么？

ELIANE! 发表于 2011-12-15 16:29
谢谢，但是我不明白prob2

rbind就是把你的结果按照行一行行的累加起来，一般用在不知道最后的矩阵大小的情况下。

qoiqpwqr 发表于 2011-12-15 23:21
rbind就是把你的结果按照行一行行的累加起来，一般用在不知道最后的矩阵大小的情况下。

就是说rbind是把两个矩阵按照行累加起来，但是其中一个是空的，其实就等于已知的那个矩阵，是么？