我正在R中执行回归分析,该分析看起来如下:
lm_carclass_mod <- lm(log(count_faves+1)~log(views+1)+dateadded+group_url+license+log(precontext.nextphoto.views+1)+log(precontext.prevphoto.views+1)+log(oid.Bridge+1)+log(oid.Face+1)+log(oid.Quail+1)+log(oid.Sky+1)+log(oid.Car+1)+log(oid.Auditorium+1)+log(oid.Font+1)+log(oid.Lane+1)+log(oid.Bmw+1)+log(oid.Racing+1)+log(oid.Wheel+1),data=flickrcar_wo_country)
confint(lm_carclass_mod,level=0.95)
summary(lm_carclass_mod)
在整个分析过程中,因变量以及一些自变量都是相当可变的,这就是为什么我想继续手动插入它们的原因。
但是,我正在寻找一种用单个函数替换所有“ oid。...”变量的方法。
到目前为止,我已经提出了以下:
g <- paste("log(",variables,"+1)", collapse="+")
不幸的是,这在lm()函数内部不起作用。这样的公式也不会:
g <- as.formula(
paste("log(",variables,"+1)", collapse="+")
)
向量变量中包含以下元素:
variables <- ("oid.Bridge", "oid.Face", "oid.Quail", "oid.Off-roading", "oid.Sky", "oid.Car", "oid.Auditorium", "oid.Font", "oid.Lane", "oid.Bmw", "oid.Racing", "oid.Wheel")
在结束中,我的回归模型应如下所示:
lm_carclass_mod <- lm(log(count_faves+1)~log(views+1)+dateadded+group_url+license+log(precontext.nextphoto.views+1)+log(precontext.prevphoto.views+1)+g,data=flickrcar_wo_country)
confint(lm_carclass_mod,level=0.95)
summary(lm_carclass_mod)
感谢您的帮助!
答案 0 :(得分:1)
您需要将两个部分都转换为字符串,然后制定公式:
#the manual bit
manual <- "log(count_faves+1)~log(views+1)+dateadded+group_url+license+log(precontext.nextphoto.views+1)+log(precontext.prevphoto.views+1)"
#the variables:
oid_variables <- c("oid.Bridge", "oid.Face", "oid.Quail", "oid.Off-roading", "oid.Sky", "oid.Car", "oid.Auditorium", "oid.Font", "oid.Lane", "oid.Bmw", "oid.Racing", "oid.Wheel")
#paste them together
g <- paste("log(", oid_variables, "+1)", collapse="+")
#make the formula
myformula <- as.formula(paste(manual, '+', g))
然后将公式添加到lm
中:
lm_carclass_mod <- lm(myformula, data=flickrcar_wo_country