using GLM using DataFrames auto = readtable("auto.csv"); showcols(auto) olsest = glm(price~mpg+weight,auto,Normal(),IdentityLink()) coef(olsest) stderr(olsest) [coef(olsest) stderr(olsest)] deviance(olsest) pricehat = predict(olsest) [convert(Array,auto[:,:price]) pricehat] vcov(olsest) [sqrt(diag(vcov(olsest))) stderr(olsest)] logitest = glm(foreign~mpg+weight,auto,Binomial(),LogitLink()) [convert(Array,auto[:,:foreign]) predict(logitest)] probitest = glm(foreign~mpg+weight,auto,Binomial(),ProbitLink()) [convert(Array,auto[:,:foreign]) predict(probitest)] olsest2 = glm(price~mpg+weight+rep78,auto,Normal(),IdentityLink()) pool!(auto,:rep78); levels(auto[:rep78]) olsest3 = glm(price~mpg+weight+rep78,auto,Normal(),IdentityLink()) auto2 = readtable("auto.csv",makefactors=true);