x <- read.csv("数量化１.csv")
x

attach(x)
s1 <- lm(y~factor(x1)+factor(x2),data=x)
summary(s1)

#標準偏回帰係数（β）を求める
lm.Beta3 = function(object) {
  d = model.matrix(object$terms, eval(object$model, parent.frame()))
  object$coefficients[-1] * apply(d[, -1, drop=FALSE], 2, sd) / sd(object$model[,1])
}

lm.Beta3(s1)