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:
图(火炬松)#一个为每一个病人,因为情节数据被分组
你想要一个阴谋显示所有组或多个阴谋吗? –
多个地块。但是,我需要用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
解决方案:total.new < - #创建一个新的拷贝对象 groupedData(ydatalog〜xdata | ID, data = as .data.frame(total),.... – Carlos