Повторные измерения Anova в r с использованием средних размеров зерна - PullRequest
Новая
       73

Повторные измерения Anova в r с использованием средних размеров зерна

0 голосов
/ 21 ноября 2018

Я пытаюсь запустить повторные измерения anova, используя r studio.В моем распоряжении средние размеры зерен песка, которые были собраны в течение 12 лет на девяти участках.Образцы были собраны с 4 участков пляжа: Foreshore, Backshore, MidBeach, Dune и сохранены в формате .csv.Я проверяю, чтобы увидеть, есть ли различия между средними размерами зерна в зависимости от места и положения берега.Пример: я хочу провести повторные измерения анова средних размеров зерна для "Foreshore" между всеми девятью участками.Я не слишком уверен, правильно ли я пишу код, и хотел бы получить некоторую помощь, если это возможно. Песочные данные 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.) Как бороться с этой ошибкой во втором фрагменте кода?

Заранее спасибо

...