Displaying 2 results from an estimated 2 matches for "sdrats".
Did you mean:
drats
2010 Mar 22
0
using lmer weights argument to represent heteroskedasticity
...## lme, heteroskedastic model with variance weights estimated
res[["lme.esthet"]] <- sumLME( lme(fixed = y ~ x, data = d, random =
~1 | groupid, weights = varIdent(form = ~1|verrorf)) )
## lme, heteroskedastic model with true fixed weights
v <- c(v0.5 = 0.5, v1 = 1, v2 = 2, v4 = 4)
sdrats <- sqrt(v / v[1])[-1]
res[["lme.fixedhet"]] <- sumLME( lme(fixed = y ~ x, data = d, random =
~1 | groupid, weights = varIdent(form = ~1|verrorf, fixed = sdrats)) )
detach("package:nlme")
library(lme4)
## lmer, ignoring heteroskedasticity
res[["lmer.nohet"]] <...
2013 Feb 20
0
Bayesian mixing model
Fellow R users,
I'm using the BCE {BCE} function to run a Bayesian sediment mixing model. The aim is to find the optimum % contribution from each of the 4 source areas that can yield the target geochemistry.
I have geochemistry for 4 source areas called Rat:
Rat<-read.table(text="CaO MgO Na2O Al2O3
Topsoils 2.511250 0.7445500 0.7085500 14.10375
ChannelBanks