Функция SAS PRO C SCORE в R - PullRequest
Новая
       39

Функция SAS PRO C SCORE в R

0 голосов
/ 12 февраля 2020

Я пытаюсь запустить командный файл для применения канонического дискриминантного анализа без размера от SAS до R [Size-Free CDA 1 . Но я понятия не имею, как интерпретировать PRO C SCORE в r

Это мой проблемный c шаг в SAS 

proc princomp cov n=1 data=sashelp.iris out=prin;
var SepalLength SepalWidth PetalLength PetalWidth;

proc reg data=prin outest=coefregs;
r1: Model SepalLength = prin1;
r2: Model SepalWidth = prin1;
r3: Model PetalLength = prin1;
r4: Model PetalWidth = prin1;

proc score type= parms residual data=prin score=coefregs out= resid;
var SepalLength SepalWidth PetalLength PetalWidth prin1;

И это мой скрипт в R для этого шаг, без PRO C SCORE 

pc <- princomp(iris[1:4], cor = F)
scores <- pc[["scores"]]
prin1 <- as.data.frame(scores[,1])

coefregs <- lm(cbind(Sepal.Length, Sepal.Width, Petal.Length, Petal.Width) ~ prin1[,1], 
data=iris) 

glht???

Я до сих пор не понимаю, как эта функция работает в SAS, поэтому я был бы признателен, если бы вы могли помочь мне с некоторым кодом в r, который mimi c PRO C SCORE

Примечание: данные диафрагмы SAS указаны в мм, а R - см

...