dev/fittingAssociationFlipped.R

###
#first test lm with fixed intercept
x<-c(1:100)
a<-10
b<-20
y<-a*x+b
y<-y+rnorm(length(y),0,0.5)

fit.full<-lm(y~x)

##try fixed int
fit.fixIntercept<-lm(y~x+1, offset=rep(20,length(y)))
ffeng23/ADASPR documentation built on July 13, 2019, 1:15 p.m.