Displaying 2 results from an estimated 2 matches for "soil_id".
Did you mean:
mail_id
2007 Jan 28
1
extra panel arguments to plot.nmGroupedData {nlme}
Greetings,
I have a groupedData (nmGroupedData) object created with the following syntax:
Soil <- groupedData(
ksat ~ conc | soil_id/sar/rep,
data=soil.data,
labels=list(x='Solution Concentration', y='Saturated Hydraulic Conductivity'),
units=list(x='(cmol_c)', y='(cm/s)')
)
the original data represents longitudinal observations in the form of:
'data.frame': 1197 obs. of 5 variabl...
2007 Jul 25
0
DF and intercept term meaning for mixed (lme) models
...DF (DF in the denomenator), as well as the
Intercept term in the output.
For example:
I have a groupedData object called 'Soil', and am fitting an lme model as
follows:
## fit a simple model
# errors partitioned among replicates
fit1 <- lme(
? ?log(ksat) ~ log(conc) + ordered(sar) + soil_id ,
? ?random = ~ 1 | rep,
? ?data=Soil
)
## check significance of model terms
anova(fit1)
? ? ? ? ? ? ? numDF denDF ?F-value p-value
(Intercept) ? ? ?1 ?1253 64313.21 ?<.0001
log(conc) ? ? ? ?1 ? 597 ? 173.34 ?<.0001
ordered(sar) ? ? 2 ? 597 ? ?13.87 ?<.0001
soil_id ? ? ? ? 29 ? 597 ? ?54...