Question

我正在尝试使用R Studio运行重复测量方差分析。我拥有12年内从9个地点收集的平均沙粒尺寸。从海滩的四个部分（前岸，后岸，中海滩，沙丘）收集了样本，并以.csv格式保存。我正在测试，以查看根据地点和岸边位置，平均晶粒度是否随时间变化。例如：我想对所有9个地点之间的“前滨”进行平均粒度的重复测量方差分析。我不太确定我是否真的在正确编写此代码，并在可能的情况下寻求帮助。Sand Data p

这是我的代码行，未指定海滩的一部分

```{r}
Years <- c(1,2,3,4,5,6,7,8,9,10,11,12)
YearFactor <- as.factor(Years)
YearFrame <- data.frame(YearFactor)
YearBind <- cbind(SandData$X2001, SandData$X2002, SandData$X2004, SandData$X2006, SandData$X2007, SandData$X2010, SandData$X2012, SandData$X2014, SandData$X2015, SandData$X2016, SandData$X2017, SandData$X2018)
library(car)
SandDataModel <- lm(YearBind ~ 1)
summary(SandDataModel)
```

这是我用来进行相同分析的代码，但仅测试了前滨的平均谷物值。

```{r}
SandData_Foreshore <- subset(SandData, ShorePosition == "Foreshore", select= c("Site", "ShorePosition", "X2001", "X2002", "X2004", "X2006", "X2007", "X2010", "X2012", "X2014", "X2015", "X2016", "X2017", "X2018"))
SandData_Foreshore
YearBind_Foreshore <- cbind(SandData_Foreshore$X2001, SandData_Foreshore$X2002, SandData_Foreshore$X2004, SandData_Foreshore$X2006, SandData_Foreshore$X2007, SandData_Foreshore$X2010, SandData_Foreshore$X2012, SandData_Foreshore$X2014, SandData_Foreshore$X2015, SandData_Foreshore$X2016, SandData_Foreshore$X2017, SandData_Foreshore$X2018)
SandData_Foreshore <- lm(YearBind_Foreshore ~ 1)
analysis_Foreshore <- Anova(SandData_Foreshore, idata = YearFrame, idesign = ~YearFactor)
summary(analysis_Foreshore)
```

运行上面的代码块时，我收到以下错误消息：本征错误（qr.coef（SSPE.qr，x $ SSPH），对称= FALSE）：'x'中的值无穷或缺失

我知道问题与数据集中缺少的值有关。你们都建议采取什么措施来减轻此错误？

基本上，我需要了解以下三件事。 1.）对于这种类型的分析，我的代码看起来正确吗？ 2）如果我的代码不正确，我需要采取什么步骤才能使我的分析正确？ 3.）如何在第二段代码中处理该错误？

预先感谢

使用平均粒度对r进行重复测量方差分析

0 个答案: