我有一些测量功率,风速和风向的数据。功率和风速之间存在二次多项式相关性,并且功率和风向之间可能存在正弦相关。我需要对此进行调查并获取Beta值和p值,并了解模型的适用程度。
所以我有两个可能影响力量的变量。该模型需要是广义线性模型。
我是R的新手,但我尝试了不同的方法,却找不到。我有一个矩阵,收集所有数据。任何人都可以帮助我吗?
答案 0 :(得分:1)
Generalized Linear Model
(glm)的公式例如是power ~ wind speed + wind direction
。
您可以使用以下代码执行此glm:
lmfit <- glm(power ~ wind_speed + wind_direction, dataset)
要分析模型,您可以summary(lmfit)
。有关分析数据的更多信息,请查看this page
有关glm的更多信息,请查看this链接,了解有关glm
功能的更多信息。或者搜索有关如何在R中执行线性模型的一些教程,例如this one
编辑:后续问题:
我个人会做以下事情:
formula <- power ~ wind_speed + wind_direction + windspeed2 + sinwd
lmfit <- glm(formula, dataset)
如果您执行以下操作,它也应该有效:
windspeed2 <- windspeed^2
sinwd <- sin(wind_direction)
formula <- power ~ wind_speed + wind_direction + windspeed2 + sinwd
lmfit <- glm(formula, dataset)