Displaying 2 results from an estimated 2 matches for "fit21".
Did you mean:
fit1
2023 Aug 31
1
simulating future observations from heteroscedastic fits
...t to simulate future observations from fits to heteroscedastic
data. A simple example is as follows:
(DF3_2 <- data.frame(y=c(1:3, 10*(1:3)),
gp=factor(rep(1:2, e=3))))
# I want to fit 4 models
# and simulate future observations from all 4:
fit11 <- lm(y~1, DF3_2)
fit21 <- lm(y~gp, DF3_2)
library(nlme)
(fit12 <- lme(y~1, data=DF3_2, random=~1|gp))
(fit22 <- lme(y~gp, data=DF3_2, random=~1|gp))
library(lme4)
(fit12r <- lmer(y~1+(1|gp), data=DF3_2, REML=FALSE))
(fit22r <- lmer(y~gp+(1|gp), data=DF3_2, REML=FALSE))
# I can simulate what I want for fi...
2006 Mar 08
1
function gdist, dist and vegdist in mvpart
...n+land_use+ca_old,
data=ba12, xv="p")
This works fine. Now I would like to use other dissimilarity measures as
can be found in the function dist (STATS) or vegdist (VEGAN).
De'Ath notes that gdist should be interchangeable with dist - but I
receive following error message:
> fit21 <- mvpart(dist (ba12[,18:29], meth="minkowski", diag=T,
upper=T p=2) ~
beers + slope_dem + elev_dem+ plc_dem + pr_curv+
curv+max_depth+doc_rocks+ abandon+land_use+ca_old,
data=ba12, xv="p")
Error in model.frame(formula, rownames, variables, varnames, extras,
extran...