Displaying 3 results from an estimated 3 matches for "modvar".
2011 Jun 27
1
group interaction in a varying coeff. model (mgcv)
Dear UseRs,
I built varying coefficient models (in mgcv) for two groups separately, with one explanatory and one moderator variable (see the example below).
# -------
# Example:
# ------
# generate moderator variable (can the same for both groups)
modvar <- c(1:1000)
# generate group1 values
x1 <- rnorm(1000)
y1 <- scale(cbind(1,poly(modvar,2))%*%c(1,2,1)*x1 +
rnorm(1000,0,0.3))
# generate group2 values
x2 <- rnorm(1000)
y2 <- scale(cbind(1,poly(modvar,2))%*%c(-1,0.5,-1)*x2 +
rnorm(1000,0,0.3))
# separate models for each group
m...
2001 Aug 21
2
Problem using GLM in a loop
...message="Null deviance in initialisation step")
endformula_ strsplit(formula,"~")[[3]]
formulavar_paste("d","~")
formulavar_as.formula(paste(formulavar,endformula))
#
# Modeling loop
#
while ((iter <= itercrit)| (vardev > devcrit)) {
d_(response - mu)^2
modvar_glm(formula=formulavar,family=Gamma(link=log),data=data)
sig <- predict(modvar,type="response")
weights <- 1/sig
modesp <- glm(formula=formula,
family=gaussian,data=data,weights=weights)
mu <- predict(modesp,type="response")
devold <- dev
dev <- modesp$deviance...
2001 Aug 21
0
Re: [R] Problem using GLM in a loop (fwd)
...step")
> endformula_ strsplit(formula,"~")[[3]]
> formulavar_paste("d","~")
> formulavar_as.formula(paste(formulavar,endformula))
> #
> # Modeling loop
> #
> while ((iter <= itercrit)| (vardev > devcrit)) {
> d_(response - mu)^2
> modvar_glm(formula=formulavar,family=Gamma(link=log),data=data)
> sig <- predict(modvar,type="response")
> weights <- 1/sig
> modesp <- glm(formula=formula,
> family=gaussian,data=data,weights=weights)
> mu <- predict(modesp,type="response")
> devold <-...