Я пытаюсь выполнить многомерный тест на нормальность некоторых данных плотности с пяти участков, используя mshapiro.test из пакета mvnormtest
. Каждый сайт представляет собой столбец, а значения плотности указаны ниже. Это 5 столбцов и 5 строк, причем верхняя строка является заголовком (названия сайтов). Вот как я загрузил свои данные:
datafilename="/Users/megsiesiple/Documents/Lisa/lisadensities.csv"
data.nc5=read.csv(datafilename,header=T)
attach(data.nc5)`
Данные выглядят так:
B07 B08 B09 B10 M
1 72571.43 17714.29 3142.86 22571.43 8000.00
2 44571.43 46857.14 49142.86 16857.14 7142.86
3 54571.43 44000.00 26571.43 6571.43 17714.29
4 57714.29 38857.14 32571.43 2000.00 5428.57
Когда я вызываю mshapiro.test () для data.nc5, я получаю следующее сообщение: Ошибка в mshapiro.test (data.nc5) :
U [] не является матрицей с числом столбцов (размером выборки) от 3 до 5000
Я знаю, что для выполнения теста Шапиро-Уилка с использованием mshapiro.test ()
данные должны быть в числовой матрице с числом столбцов от 3 до 5000. Однако даже когда я делаю .csv матрицей только с числами (т. е. когда я опускаю имена сайтов), я все равно получаю ошибку. Мне нужно настроить матрицу по-другому? У кого еще возникла эта проблема?
Спасибо!