... beta sum(dnorm(y, mean = linpred, sd = sigma, log = TRUE)) } ### Initialize the search ls.reg <— lm(y ~ X — 1) beta.hat.ls <— coef(ls.reg) sigma.hat.ls ...
確定! 回上一頁