[問題] 請問如何寫出l<- lm(weight ~ group)?

[ml919 (ml919)]

請教各位!

請問要如何寫出類似這樣的function, 以linear regression為例:
lm.D9 <- lm(weight ~ group)
lm.D90 <- lm(weight ~ group - 1)
anova(lm.D9)
summary(lm.D90)

或是像mgcv package
lr.fit1 <- gam(y~s(x0)+s(x1)+s(x2),family=binomial,
               data=dat,method="REML")

有自變數, 有因變數, 有 ~

請問這種函式該怎麼寫? 有什麼書籍可以參考嗎?
或google搜尋, 關鍵字應該打什麼?

感謝!!!!

--
※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 173.91.250.199 (美國)
※ 文章網址: https://www.ptt.cc/bbs/R_Language/M.1708386943.A.4B7.html
→ andrew43: 我猜你可以看看 help(model.extract) 裡面例子
→ andrew43: 它示範了怎麼從formula+data取出變數
→ locka: 是指想自創用formula方式定義function參數嗎？
→ ml919: 對的, 所以我應該往formula方向做嗎? 謝謝!
→ andrew43: 你可以舉一個例子來說明你的目標是什麼
推 sacidoO: 感覺原po想寫一個function return with 獨立變數Xi as i
→ sacidoO: nputs