# transform from triangle into long format
chrisl <- expand.grid(dimnames(chris))
chrisl$value <- as.vector(chris)
chrisl <- na.omit(chrisl)
lm(log(value) ~ origin + dev, data=chrisl)
model.matrix(log(value) ~ dev + origin-1, data=chrisl)
co=summary(lm(log(value) ~ origin + dev, data=chrisl))$coef[,1]
c(co[1], co[1]+co[2], co[1]+co[3], co[1]+co[4], co[5:7])
Run the code above in your browser using DataLab