inst/snippet/seatpos-pc.R

pc=with(seatpos,princomp(cbind(HtShoes,Ht,Seated,Arm,Thigh,Leg),
    scores=T))
summary(pc, loadings=T)
seatpos.lmpc <-lm(hipcenter ~ Age + Weight + pc$scores[,1], seatpos )
###hop:3-9
summary(seatpos.lmpc)
vif(seatpos.lmpc)

Try the fastR package in your browser

Any scripts or data that you put into this service are public.

fastR documentation built on May 1, 2019, 11:18 p.m.