我正在使用数据表来处理大量线性模型,但我无法使用公式界面。 一个例子如下:
require(data.table)
dt <- data.table(grp = c(rep("a", 10), rep("b",10)),
x1 = rnorm(20),
x2 = rnorm(20),
y = rnorm(20))
n <- names(dt)
modFormula <- paste0("y ~ ", paste0(n[!n %in% c("y", "grp")], collapse = "+"))
#This gives me an error - Error in eval(expr, envir, enclos) : object 'y' not found
results <- dt[, as.list(coef(lm(modFormula))), by = .(grp)]
#This works but I ideally don't want to be typing out the variables manually
results <- dt[, as.list(coef(lm(y~x1+x2))), by = .(grp)]
如何动态传递公式/ x,y变量?
我正在通过成千上万的模型进行搅拌,这是迄今为止我发现的最快的方式 - 欢迎任何其他建议。
答案 0 :(得分:3)
我们需要使用<h1 id = "FirstHeader">Coffes here</h1>
<img src="coffee.png" id = "Coffee" alt="coffee image" width="100" height="125"/>
<div id = "bodytext">
<div id = "FirstParagraph">
Duis maximus enim in vestibulum posuere. Nullam tempor, est in rutrum ullamcorper,
Vivamus ut mi orci. Donec ut ipsum in nisl elementum posuere mollis id massa.
Cras mattis felis est, a varius nisl euismod ac. Mauris cursus sed nulla in viverra.
Duis maximus enim in vestibulum posuere. Nullam tempor, est in rutrum ullamcorper,
purus purus sodales purus, sit amet imperdiet urna velit id velit. Proin vel dapibus
est. Nulla vel tortor id ligula consequat fringilla dignissim a sem. Duis eget risus vel
est ullamcorper auctor. Aliquam vel ex in leo suscipit condimentum sed vitae leo. In porta
eget augue eu porta. Pellentesque tellus eros, consequat vel aliquet id, consectetur a ante.
Vestibulum sit amet ligula faucibus ligula auctor pretium non vel arcu.<br>
purus purus sodales purus, sit amet imperdiet urna velit id velit. Proin vel dapibus
Vivamus ut mi orci. Donec ut ipsum in nisl elementum posuere mollis id massa.
Cras mattis felis est, a varius nisl euismod ac. Mauris cursus sed nulla in viverra.
Duis maximus enim in vestibulum posuere. Nullam tempor, est in rutrum ullamcorper,
purus purus sodales purus, sit amet imperdiet urna velit id velit. Proin vel dapibus
est. Nulla vel tortor id ligula consequat fringilla dignissim a sem. Duis eget risus vel
est ullamcorper auctor. Aliquam vel ex in leo suscipit condimentum sed vitae leo. In porta
eget augue eu porta. Pellentesque tellus eros, consequat vel aliquet id, consectetur a ante.
Vestibulum sit amet ligula faucibus ligula auctor pretium non vel arcu.purus purus sodales purus, sit amet imperdiet urna velit id velit. Proin vel dapibus
Vivamus ut mi orci. Donec ut ipsum in nisl elementum posuere mollis id massa.
Cras mattis felis est, a varius nisl euismod ac. Mauris cursus sed nulla in viverra.
Duis maximus enim in vestibulum posuere. Nullam tempor, est in rutrum ullamcorper,
purus purus sodales purus, sit amet imperdiet urna velit id velit. Proin vel dapibus
est. Nulla vel tortor id ligula consequat fringilla dignissim a sem. Duis eget risus vel
est ullamcorper auctor. Aliquam vel ex in leo suscipit condimentum sed vitae leo. In porta
eget augue eu porta. Pellentesque tellus eros, consequat vel aliquet id, consectetur a ante.
Vestibulum sit amet ligula faucibus ligula auctor pretium non vel arcu.
</div>
<div id = "SecondParagraph">
purus purus sodales purus, sit amet imperdiet urna velit id velit. Proin vel dapibus
Vivamus ut mi orci. Donec ut ipsum in nisl elementum posuere mollis id massa.
Cras mattis felis est, a varius nisl euismod ac. Mauris cursus sed nulla in viverra.
Duis maximus enim in vestibulum posuere. Nullam tempor, est in rutrum ullamcorper,
purus purus sodales purus, sit amet imperdiet urna velit id velit. Proin vel dapibus
est. Nulla vel tortor id ligula consequat fringilla dignissim a sem. Duis eget risus vel
est ullamcorper auctor. Aliquam vel ex in leo suscipit condimentum sed vitae leo. In porta
eget augue eu porta. Pellentesque tellus eros, consequat vel aliquet id, consectetur a ante.
Vestibulum sit amet ligula faucibus ligula auctor pretium non vel arcu.
<br/>
purus purus sodales purus, sit amet imperdiet urna velit id velit. Proin vel dapibus
Vivamus ut mi orci. Donec ut ipsum in nisl elementum posuere mollis id massa.
Cras mattis felis est, a varius nisl euismod ac. Mauris cursus sed nulla in viverra.
Duis maximus enim in vestibulum posuere. Nullam tempor, est in rutrum ullamcorper,
purus purus sodales purus, sit amet imperdiet urna velit id velit. Proin vel dapibus
est. Nulla vel tortor id ligula consequat fringilla dignissim a sem. Duis eget risus vel
est ullamcorper auctor. Aliquam vel ex in leo suscipit condimentum sed vitae leo. In porta
eget augue eu porta. Pellentesque tellus eros, consequat vel aliquet id, consectetur a ante.
Vestibulum sit amet ligula faucibus ligula auctor pretium non vel arcu.purus purus sodales purus, sit amet imperdiet urna velit id velit. Proin vel dapibus
Vivamus ut mi orci. Donec ut ipsum in nisl elementum posuere mollis id massa.
Cras mattis felis est, a varius nisl euismod ac. Mauris cursus sed nulla in viverra.
Duis maximus enim in vestibulum posuere. Nullam tempor, est in rutrum ullamcorper,
purus purus sodales purus, sit amet imperdiet urna velit id velit. Proin vel dapibus
est. Nulla vel tortor id ligula consequat fringilla dignissim a sem. Duis eget risus vel
est ullamcorper auctor. Aliquam vel ex in leo suscipit condimentum sed vitae leo. In porta
eget augue eu porta. Pellentesque tellus eros, consequat vel aliquet id, consectetur a ante.
Vestibulum sit amet ligula faucibus ligula auctor pretium non vel arcu.
</div>
<div id = "ThirdParagraph">
eget augue eu porta. Pellentesque tellus eros, consequat vel aliquet id, consectetur a ante.
Vivamus ut mi orci. Donec ut ipsum in nisl elementum posuere mollis id massa.
Cras mattis felis est, a varius nisl euismod ac. Mauris cursus sed nulla in viverra.
Duis maximus enim in vestibulum posuere. Nullam tempor, est in rutrum ullamcorper,
purus purus sodales purus, sit amet imperdiet urna velit id velit. Proin vel dapibus
est. Nulla vel tortor id ligula consequat fringilla dignissim a sem. Duis eget risus vel
est ullamcorper auctor. Aliquam vel ex in leo suscipit condimentum sed vitae leo. In porta
eget augue eu porta. Pellentesque tellus eros, consequat vel aliquet id, consectetur a ante.
Vestibulum sit amet ligula faucibus ligula auctor pretium non vel arcu.purus purus sodales purus, sit amet imperdiet urna velit id velit. Proin vel dapibus
Vivamus ut mi orci. Donec ut ipsum in nisl elementum posuere mollis id massa.
Cras mattis felis est, a varius nisl euismod ac. Mauris cursus sed nulla in viverra.
Duis maximus enim in vestibulum posuere. Nullam tempor, est in rutrum ullamcorper,
purus purus sodales purus, sit amet imperdiet urna velit id velit. Proin vel dapibus
est. Nulla vel tortor id ligula consequat fringilla dignissim a sem. Duis eget risus vel
est ullamcorper auctor. Aliquam vel ex in leo suscipit condimentum sed vitae leo. In porta
eget augue eu porta. Pellentesque tellus eros, consequat vel aliquet id, consectetur a ante.
Vestibulum sit amet ligula faucibus ligula auctor pretium non vel arcu.
<br>
purus purus sodales purus, sit amet imperdiet urna velit id velit. Proin vel dapibus
Vivamus ut mi orci. Donec ut ipsum in nisl elementum posuere mollis id massa.
Cras mattis felis est, a varius nisl euismod ac. Mauris cursus sed nulla in viverra.
Duis maximus enim in vestibulum posuere. Nullam tempor, est in rutrum ullamcorper,
purus purus sodales purus, sit amet imperdiet urna velit id velit. Proin vel dapibus
est. Nulla vel tortor id ligula consequat fringilla dignissim a sem. Duis eget risus vel
est ullamcorper auctor. Aliquam vel ex in leo suscipit condimentum sed vitae leo. In porta
eget augue eu porta. Pellentesque tellus eros, consequat vel aliquet id, consectetur a ante.
Vestibulum sit amet ligula faucibus ligula auctor pretium non vel arcu.
</div>
</div>
<form id = "MainForm">
<label for = "Name">Name: </label>
<input type="text" id = "Name" placeholder = "e.g nirajan basnet" required/>
<label for= "Email">Email</label>
<input type="email" id = "Email" placeholder="sth@sthmail.com" required/>
<label for = "Comment">Your comment: </label>
<br>
<textarea maxlength="100" id = "Comment" placeholder="What do you think about the site?"></textarea>
</form>
data
lm
使用第二种方法检查输出
dt[, as.list(coef(lm(modFormula, data = .SD))), by = .(grp)]
# grp (Intercept) x1 x2
#1: a 0.2185345 -0.7334866 -0.5986349
#2: b 0.4499405 -0.5264809 0.1454797