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