Displaying 3 results from an estimated 3 matches for "fit_2".
Did you mean:
fit2
2011 Dec 17
2
Problem with reproducing log likelihood estimated with ghyp package
...$nu,0)))
#The two log likelihood estimates is approximatly equal, and the parameters
are
#sensible.
#When I check series 2, i get a very different result. The estimate for mu
and nu
#is still sensible, but the log likelihood is very different from what i get
with
#dt(x, df, ncp, log = FALSE)
fit_2=fit.tuv(series_2, silent=T, symmetric=T)
fit_2
sum(log(dt(series_2,coef(fit_2)$nu,0)))
#This is very different
--
View this message in context: http://r.789695.n4.nabble.com/Problem-with-reproducing-log-likelihood-estimated-with-ghyp-package-tp4207833p4207833.html
Sent from the R help mailing l...
2009 Jun 22
1
Problem with storing a sequence of lmer() model fit into a list
...lem.
#the following code can store glm() model fit into a list
---------------------------------------------------
x1<-runif(200)
x2<-rnorm(200)
y<-x1+x2
testdf<-data.frame(y=y, x1=x1, x2=x2)
indepvec<-c("x1","x2")
fit.list<-NULL
fit_1<-glm(y~x1,data=testdf)
fit_2<-glm(y~x2,data=testdf)
fit.list[[paste('fit_',indepvec[1],sep='')]]<-fit_1
fit.list[[paste('fit_',indepvec[12],sep='')]]<-fit_2
----------------------------------------------------
so why cannot I store lmer() model fit in a list?
Would someone kindly explai...
2005 Nov 03
4
nlme questions
Dear R users;
Ive got two questions concerning nlme library 3.1-65 (running on R 2.2.0 /
Win XP Pro). The first one is related to augPred function. Ive been working
with a nonlinear mixed model with no problems so far. However, when the
parameters of the model are specified in terms of some other covariates,
say treatment (i.e. phi1~trt1+trt2, etc) the augPred function give me the
following