search for: fit_2

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