2015-05-14 86 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),。