speedglm

    3熱度

    1回答

    我試圖用speedglm實現比glm更快的GLM估計,但爲什麼它更慢? set.seed(0) n=1e3 p=1e3 x=matrix(runif(n*p),nrow=n) y=sample(0:1,n,replace = T) ptm <- proc.time() fit=glm(y~x,family=binomial()) print(proc.time() - ptm)

    1熱度

    1回答

    我想運行加權數據的線性迴歸。 當使用speedlm時,當數據中存在缺失值時,我收到錯誤消息。 library(speedglm) sampleData <- data.frame(w = round(runif(12,0,1)), target = rnorm(12,100,50), predictor = c(NA, rnorm(10, 40, 10