实验数据处理R语言报错-经管之家官网！

扫码加入统计交流群

library("reshape")

library("readr")

data1<-read.table("expb.txt","\t",header= T)

inteval<-100;totaltime<-5000 #设定总观察时间和时间分段

subid<-unique(data1$sub)

data2<-subset(data1,start<totaltime)

for (i in 1:length(subid)) {

datatmp<-subset(data2,sub==subid[1])

for(l in 1:length(datatmp[,1])){

if(datatmp[l,"end"]>=totaltime)

datatmp[l,"end"]<-totaltime-1

}

dataYN<-data.frame(Subject=NA,Item=NA,Condition=NA,Fixseq=NA,area1=NA,area2=NA,area3=NA,area4=NA,area5=NA)

for(t in 1:48){

datatmptmp<-subset(datatmp,trialid==t)

dataY<-data.frame(subject=NA,item=NA,condition=NA,fixseq=NA,area1=NA,area2=NA,area3=NA,area4=NA,area5=NA)

for(m in 1:max(datatmptmp[,"end"]))

{dataY[m,c(1,2,3,4)]<-c(datatmptmp[1,c("sub","trialid","cond")],as.integer(m/inteval)+1)}

for(k in 1:length(datatmptmp[,1]))

{

for(n in datatmptmp[k,"start"]:datatmptmp[k,"end"])

{dataY[n,4+datatmptmp[k,"stimuli"]]<-as.integer(n>=datatmptmp[k,"start"]&&n<=datatmptmp[k,"end"])}

}

dataY <- aggregate(dataY, by=list(Subject = dataY$subject,Item =dataY$item,Condition = dataY$condition,Fixseq = dataY$fixseq), FUN=sum,na.rm=TRUE)

dataY$item <- dataY$condition <- dataY$subject<-dataY$fixseq<- NULL

dataYN<-rbind(dataYN,dataY)

}

dataYN<-subset(dataYN,!is.na(Subject))

dataYNtmp<-dataYN[,c(5:8)];

dataYNtmp<-log((dataYNtmp+0.5)/(inteval-dataYNtmp+0.5));

#dataYN$datalgd<-dataYNtmp[,2]-dataYNtmp[,4];

dataYN[,c(9:12)]<-dataYNtmp;

write_csv(dataYN,"datalin.csv",append= T)

}

出现报错：Error in sub == subid[1] : 只能比较(1)基元或串列种类

「经管之家」APP：经管人学习、答疑、交友，就上经管之家！
免流量费下载资料----在经管之家app可以下载论坛上的所有资源，并且不额外收取下载高峰期的论坛币。
涵盖所有经管领域的优秀内容----覆盖经济、管理、金融投资、计量统计、数据分析、国贸、财会等专业的学习宝库，各类资料应有尽有。
来自五湖四海的经管达人----已经有上千万的经管人来到这里，你可以找到任何学科方向、有共同话题的朋友。
经管之家（原人大经济论坛），跨越高校的围墙，带你走进经管知识的新世界。
扫描下方二维码下载并注册APP