2015-05-14 46 views
1

我有100个患者数据(1个带有rbind的总文件)和2列(时间,var1)。 我需要组数据(100组),以便能够绘制100个病人的数据。按ID或种子分组数据

我在最后添加了一列patientID到整个数据,但应该分组以便为每个患者生成100个图。由于我需要这些数据,我不想分割总数据。

我尝试这样做,也没有工作:

Data1$ID <- 1 # third column now has the ID 
total <- rbind(Data1,Data2,....) 
total <- grouped.data(xdata = xdata, ydatalog = ydatalog, seed = total$ID) 
plot(total) 

举个例子,从读该数据正常工作:

str(Loblolly) 
Classes ‘nfnGroupedData’, ‘nfGroupedData’, ‘groupedData’ and 'data.frame':  84 obs. of 3 variables: 

图(火炬松)#一个为每一个病人,因为情节数据被分组

+0

你想要一个阴谋显示所有组或多个阴谋吗? –

+0

多个地块。但是,我需要用nmle分析整个数据(所以应该按照patiend ID进行分组)。例如:fm1 <-nlme(height_asym +(5 * R0-Asym)* exp(-exp(lrc)* age), data = Loblolly, fixed = Asym + R0 + lrc_1, random = Asym 〜1, start = c(Asym = 103,R0 = -8.5,lrc = -3.3)) – Carlos

+0

解决方案:total.new < - #创建一个新的拷贝对象 groupedData(ydatalog〜xdata | ID, data = as .data.frame(total),.... – Carlos

回答

0

解决方案:total.new < - #创建一个新的拷贝对象groupedData(ydatalog〜xdata | ID,data = as.data.frame(total),。