Displaying 4 results from an estimated 4 matches for "fitm".
Did you mean:
fit
2018 Jul 20
3
Should there be a confint.mlm ?
It seems that confint.default returns an empty data.frame for objects of
class mlm. For example:
```
nobs <- 20
set.seed(1234)
# some fake data
datf <-
data.frame(x1=rnorm(nobs),x2=runif(nobs),y1=rnorm(nobs),y2=rnorm(nobs))
fitm <- lm(cbind(y1,y2) ~ x1 + x2,data=datf)
confint(fitm)
# returns:
2.5 % 97.5 %
```
I have seen proposed workarounds on stackoverflow and elsewhere, but
suspect this should be fixed in the stats package. A proposed
implementation would be:
```
# I need this to run the code, but stats does n...
2018 Jul 20
0
Should there be a confint.mlm ?
...;mlm' objects are also 'lm' objects, and so
it is confint.lm() which is called here and fails.
> For example:
>
> ```
> nobs <- 20
> set.seed(1234)
> # some fake data
> datf <-
> data.frame(x1=rnorm(nobs),x2=runif(nobs),y1=rnorm(nobs),y2=rnorm(nobs))
> fitm <- lm(cbind(y1,y2) ~ x1 + x2,data=datf)
> confint(fitm)
> # returns:
> 2.5 % 97.5 %
> ```
>
> I have seen proposed workarounds on stackoverflow and elsewhere, but
> suspect this should be fixed in the stats package.
I agree.
It may be nicer to tweak confint.lm() ins...
2008 Aug 22
2
WinBUGS with R
...t;- list("y","X","n","m") #My variables
inits.beta <- rep(0,K)
inits.beta0 <- 0
inits <- function(){list(beta=inits.beta,beta0=inits.beta0,taueps=1.0E-3)}
parameters <- list("sigma","beta","beta0","y.star")
fitm <- bugs(data,inits,parameters,model.file="model.bug",
n.chains=3, n.iter=n.iter, n.burnin=n.burnin, n.thin = n.thin,
debug=FALSE,DIC=FALSE,digit=5,codaPkg=FALSE,
bugs.directory="C:/Program Files/WinBUGS14/")
but I always get the following error:
Error in F...
2011 Mar 23
1
how to add in interaction terms in gamm
I want to use gamm to generate smoothed trend line for three groups
identified by dummy variable genea and geneb. My question is how to add in
an interaction term between the time and another dummy variable such as
gender?
fitm<-gamm(change_gfr~
genea+geneb+s(timea_n,bs="ps")+s(timeb_n,bs="ps")+s(timec_n,bs="ps"),data=mm,random=list(time_n=~1|PID))
--
View this message in context: http://r.789695.n4.nabble.com/how-to-add-in-interaction-terms-in-gamm-tp3400243p3400243.html
Sent from the R...