noInt1

test_noint1.jl

using LinearRegressionKit, StatsModels, DataFrames, CSV

y = [130, 131, 132, 133, 134, 135, 136, 137, 138, 139, 140]
x = [60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70]

df = DataFrame(x= x, y= y)

f = @formula(y ~ 0 + x )
lm= regress(f, df, req_stats=["default"])


delta(a,b) = if isapprox(abs(a-b), 0., atol=eps()) return 0. else return abs(a-b) end

println("NoInt1 - LRK")

println(string("delta B1 ", delta(2.07438016528926, lm.coefs[1])))
println(string("delta std err B1 ", delta(0.165289256198347E-01, lm.stderrors[1])))
println(string("delta resid std dev ", delta(3.56753034006338, lm.RMSE)))
println(string("delta R2 ", delta(0.999365492298663, lm.R2)))

# using GLM
println("NoInt1 - GLM")
lr = GLM.lm(f, df)

println(string("delta B1 ", delta(2.07438016528926, coef(lr)[1])))
println(string("delta std err B1 ", delta(0.165289256198347E-01, stderror(lr)[1])))
println(string("delta resid std dev ", delta(3.56753034006338, deviance(lr)/dof_residual(lr))))
println(string("delta R2 ", delta(0.999365492298663, r2(lr))))