У меня есть следующий код с помощью R:вычисление коэффициента GLM матрица, диагональные и среднеквадратическая ошибка
x0 <- matrix(rnorm(100,1))
x <- as.matrix(cbind("Intercept"=1, x0))
n <- dim(x0)[[1]]
z <- cbind(rep(1,n),x0)
p <- dim(x0)[[2]]+1
for(i in 1:n) {
gstart <- glm(y~x0,family=binomial)$coef
}
Я хочу, чтобы вычислить оценки предыдущей обобщенной линейной модели в n
образцах и создать матрицу оценок для в n
случаи, то вычисления bias
и mean square error
, где матрица параметров задаются следующим кодом:
n=100 #is the number of samples
parameter.mat<-cbind(rep(2,n),rep(2,n))
Если 'y' быть определен где-то, и' x' и 'p' не используются и 'n' определяется дважды? Это немного запутывает – user20650