Estoy tratando de ejecutar un MANOVA de medidas repetidas en R, con dos variables independientes dentro de los sujetos y tres variables dependientes.
He convertido el df en formato largo y ahora estoy intentando ejecutar el análisis utilizando el paquete RM.MANOVA. No hay valores perdidos en los datos.
head(dflong)
ID IV1 IV2 DV1 DV2 DV3
1 1 BuiltSec OffSec 2.716667 2.504202 2.403361
2 2 BuiltSec OffSec 3.158333 3.134454 3.369748
3 3 BuiltSec OffSec 3.008333 2.361345 2.983193
4 4 BuiltSec OffSec 2.683333 1.865546 2.689076
5 5 BuiltSec OffSec 2.791667 2.798319 2.747899
6 6 BuiltSec OffSec 3.033333 2.680672 3.142857
> str(dflong)
'data.frame': 264 obs. of 6 variables:
$ ID: Factor w/ 44 levels "1","2","3","4",..: 1 2 3 4 5 6 7 8 9 10 ...
$ IV1: Factor w/ 3 levels "BuiltSec","IntervSec",..: 1 1 1 1 1 1 1 1 1 1 ...
$ IV2: Factor w/ 2 levels "NoffSec","OffSec": 2 2 2 2 2 2 2 2 2 2 ...
$ DV1: num 2.72 3.16 3.01 2.68 2.79 ...
$ DV2: num 2.5 3.13 2.36 1.87 2.8 ...
$ DV3: num 2.4 3.37 2.98 2.69 2.75 ...
model.rm <- RM(cbind(dv1, dv2, dv3) ~ IV1*IV2, data = dflong, subject = "ID", no.subf=2)
Obtengo un "Error en WTS_out[i, ] <- round(results$WTS, dec) : el número de elementos a reemplazar no es un múltiplo de la longitud del reemplazo"
Parece que hay algún problema con el cálculo del índice WTS, pero no entiendo cómo resolverlo. Sabe usted qué significa este error y cómo solucionarlo?