希望将每个 ID 的唯一 ID 数量和时间段输出到 .txt 文件中。
现在我有：
#Store N and T_Max (ie. Maximum number of periods per individual)
library(plyr)
mydata$id <- as.numeric(as.factor(mydata$ncountry))
N_max <- max(mydata$id)
mydata$count<-do.call(c, lapply(unique(mydata$id), function(curf){
seq(sum(mydata$id==curf))
}))
T_max <- max(mydata$count)
mydata$T_max <- rep(T_max,nrow(mydata)) # make new column
mydata$N_max <- rep(N_max,nrow(mydata)) # make new column
myvars1 <- c("T_max","N_max")
unique(myvars1)
newdata1 <- mydata[myvars1]
write.table(newdata1, "NT_input.txt", col.names=FALSE, row.names=FALSE, sep=" ")
这给出了 .txt 文件的以下样式：
194 15
194 15
194 15
194 15
.
.
.
我只想要一行输出，即：
194 15
任何建议表示赞赏。谢谢你。