저는 두 가지 다른 시각적 지각 테스트에서 대규모 참가자 샘플을 테스트했습니다. 이제 두 테스트의 성능이 어느 정도 상관 관계가 있는지보고 싶습니다.
상관 관계를 시각화하기 위해 R을 사용하여 산점도를 플로팅 ggplot()
하고 회귀선을 적합합니다 (사용 stat_smooth()
). 그러나 my x
및 y
variable은 모두 성능 측정이므로 회귀선 stat_smooth(method="lm")
을 맞출 때 두 가지를 모두 고려해야합니다. 따라서 간단한 선형 회귀 (사용 )를 사용할 수 없지만 직교 회귀 (또는 Total 최소 제곱). 어떻게하면 되나요?
formula
에서 지정할 수 stat_smooth()
있다는 것을 알고 있지만 어떤 수식을 사용해야할지 모르겠습니다. 내가 이해 한 바에 따르면 사전 설정 방법 ( lm, glm, gam, loess, rlm
)은 적용 할 수 없습니다.
여기에 표시된대로 (x, y)의 주성분 분석에서 기울기와 절편을 추출 할 수 있습니다 . 이 기본 R에 조금 간단하고 실행하고, 사용하는 동일한 결과 제공 Deming(...)
에를 MethComp
.
# same `x and `y` as @user20650's answer
df <- data.frame(y, x)
pca <- prcomp(~x+y, df)
slp <- with(pca, rotation[2,1] / rotation[1,1])
int <- with(pca, center[2] - slp*center[1])
ggplot(df, aes(x,y)) +
geom_point() +
stat_smooth(method=lm, color="green", se=FALSE) +
geom_abline(slope=slp, intercept=int, color="blue")
이 기사는 인터넷에서 수집됩니다. 재 인쇄 할 때 출처를 알려주십시오.
침해가 발생한 경우 연락 주시기 바랍니다[email protected] 삭제
몇 마디 만하겠습니다