Displaying 3 results from an estimated 3 matches for "primula".
2010 Nov 09
2
new column from column in another df
If I have a data frame where a species occupies several rows with different
phases such as (both col's ar factors):
species,phase
Populus tremula,1
Populus tremula,2
Populus tremula,3
Calluna vulgaris,1
Calluna vulgaris,2
Betula alba,1
Betula alba,2
Betula alba,3
Primula veris,1
Primula veris,2
and another df where each species only have one row:
species,growth_form
Populus tremula,tree
Acer platanoides,tree
Ribes rubrum,shrub
Calluna vulgaris,dwarf_shrub
Betula alba,tree
Primula veris,herb
...how can I create a new column in the first data frame where growth for...
2007 Feb 13
1
lme4/lmer: P-Values from mcmc samples or chi2-tests?
...e between different p-value estimates.
How can I proceed, left with these uncertainties in the estimations of
the p-values?
Below is the corresponding R code with some output so that you can see
it all for your own:
##
m1<-lmer(number_pollinators~logpatch+loghab+landscape_diversity+(1|site),primula,poisson,method="ML")
m2<-lmer(number_pollinators~logpatch+loghab+landscape_diversity+(1|site),primula,quasipoisson,method="ML")
summary(m1);summary(m2)
#m1: [...]
Fixed effects:
Estimate Std. Error z value Pr(>|z|)
(Intercept) -0.40302 0.57...
2007 Feb 12
1
lmer and estimation of p-values: error with mcmcpvalue()
...(std * std)
sum(sqdist[-1] > sqdist[1])/nrow(samp) }
m1<-lmer(number_pollinators~logpatch+loghab+landscape_diversity+(1|site),quasipoisson)
Generalized linear mixed model fit using Laplace
Formula: number_pollinators ~ logpatch + loghab + landscape_diversity +
(1 | site)
Data: primula
Family: quasipoisson(log link)
AIC BIC logLik deviance
84.83 93.75 -37.42 74.83
Random effects:
Groups Name Variance Std.Dev.
site (Intercept) 0.036592 0.19129
Residual 1.426886 1.19452
number of obs: 44, groups: site, 15
Fixed effects:...