Buenas,
estoy consultando un artículo para estimas la correlación entre dos variables, ambas con réplicas, al final la autora pone un ejemplo con SAS, y me gustaria trasladarlo a código de R, es el siguiente:
options nocenter ls = 80 ps = 50 nodate nonumber;
data a;
input patient mvar response replicate;
lines;
1 1 6.68 1
1 2 3.97 1
1 1 6.53 2
1 2 4.12 2
. . . .
. . . .
8 1 6.77 7
8 2 6.09 7
8 1 6.82 8
8 2 5.58 8
;
proc mixed method = ml covtest;
class patient mvar replicate;
model response = mvar / solution ddfm = kr;
random mvar / type=un subject = patient v vcorr;
repeated mvar replicate/ type = un@ar(1) subject = patient r rcorr;
run;
Gracias por la ayuda