[R] Problem of intercept?
varin sacha
v@r|n@@ch@ @end|ng |rom y@hoo@|r
Tue Feb 21 22:33:54 CET 2023
Dear R-experts,
Here below my R code working with quite a few warnings.
x11 and x12 are dichotomous variable (0=no and 1=yes). I substract 1 to ignore intercept.
I would like not to ignore intercept. How to modify my R code because if I just remove -1 it does not work?
y= c(32,45,65,34,23,43,65,76,87,98,7,867,56,45,65,76,88,34,55,66)
x11=c(0,1,1,0,0,1,1,1,0,0,1,0,0,1,0,0,1,1,0,1)
x12=c(0,1,0,1,0,1,1,0,1,1,0,0,1,1,1,0,0,1,0,0)
Dataset=data.frame(y,x11,x12)
a=lm(y~x11+x12-1,Dataset)$coef
b=NULL
for(i in c(1:2)) {
f=formula(paste('y~',names(Dataset)[i],-1))
b=c(b,lm(f,Dataset)$coef)
}
coef=data.frame(rbind(a,b))
coef$Model=c('Multi','Single')
library(reshape2)
coef.long<-melt(coef,id.vars="Model")
library(ggplot2)
ggplot(coef.long,aes(x=variable,y=value,fill=Model))+
geom_bar(stat="identity",position="dodge")+
scale_fill_discrete(name="Model",
labels=c("Multiple", "Simple"))+
labs(title =paste('La différences des coefficients
entre la régression multiple et simple'),
x="Models",y="Coefficient")+
coord_flip()
More information about the R-help
mailing list