Displaying 3 results from an estimated 3 matches for "fit_1".
Did you mean:
fit1
2011 Dec 17
2
Problem with reproducing log likelihood estimated with ghyp package
...l to sqrt(2)
series_2=series_1/sqrt(2)*5
#To get student t distributed variabler with standard
#deviation equal to 5 I rescale the first series
#When i check the first series with the ghyp package, the result coincides
with
#the sum of log likelihood calculated with dt(x, df, ncp, log = FALSE)
fit_1=fit.tuv(series_1, silent=T, symmetric=T)
fit_1
sum(log(dt(series_1,coef(fit_1)$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 li...
2009 Jun 22
1
Problem with storing a sequence of lmer() model fit into a list
...s glm() model fit without a
problem.
#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 l...
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