if (FALSE) {
## (this example can take quite a while to run)
##
## Fitting two instances of a multiplicative interaction (i.e. a
## two-component interaction)
yield.scaled <- wheat$yield * sqrt(3/1000)
treatment <- factor(paste(wheat$tillage, wheat$summerCrop, wheat$manure,
wheat$N, sep = ""))
bilinear2 <- gnm(yield.scaled ~ year + treatment +
instances(Mult(year, treatment), 2),
family = gaussian, data = wheat)
}
Run the code above in your browser using DataLab