Displaying 2 results from an estimated 2 matches for "denomenator".
2007 Jul 25
0
DF and intercept term meaning for mixed (lme) models
Hi,
I am using the lme package to fit mixed effects models to a set of data.
I am having a difficult time understanding the *meaning* of the numDF (degrees
of freedom in the numerator), denDF (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...
2013 Apr 03
1
Problem with integrate function
...)/(0.1+5*x[i]))^2}
L<-optim(c(0),kopt,method="BFGS")
N[i]<-round(L$par+1)}
return(N)
} # Calculate Kv+1
fx[[2]]<-function(x)((0.2+6*x)/(.3+4*x))^n*(1-(0.2+6*x)/(.3+4*x))/(1-(0.2+6*x)/(.3+4*x))^fx[[1]](x)*fv[[1]](x)
# numerator of f'
gx<-integrate(fx[[2]],0,0.045) #denomenator of f'
hx<-integrate(fx[[2]],0.056,1)
fv[[2]]<-function(x)x*fx[[2]](x)/(gx$value+hx$value) #v.f'
fx[[3]]<-function(x){
NV<-0;n0<-0;
rho<-(0.2+6*x)/(.3+4*x)
i<-0
N<-0
for(i in 1:length(x)){
kopt<-function(K){(x[i]-0.05*(K-((1-rho[i]^K)/(1-rho[i])+(K*rho[...