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

作者: ml919 (ml919)   2024-02-20 07:55:41
請教各位!
請問要如何寫出類似這樣的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搜尋, 關鍵字應該打什麼?
感謝!!!!
作者: andrew43 (討厭有好心推文後刪文者)   2024-02-20 10:42:00
我猜你可以看看 help(model.extract) 裡面例子它示範了怎麼從formula+data取出變數
作者: locka (locka)   2024-02-20 12:21:00
是指想自創用formula方式定義function參數嗎?
作者: ml919 (ml919)   2024-02-20 22:25:00
對的, 所以我應該往formula方向做嗎? 謝謝!
作者: andrew43 (討厭有好心推文後刪文者)   2024-02-21 17:17:00
你可以舉一個例子來說明你的目標是什麼
作者: sacidoO (阿罵)   2024-03-06 22:20:00
感覺原po想寫一個function return with 獨立變數Xi as inputs

Links booklink

Contact Us: admin [ a t ] ucptt.com