2014-11-20 1 views
0

Moran의 I 분석을 실행하여 피쳐 간의 공간적 관계를 찾습니다. 분석은 ncf R 패키지의 correlog 기능을 사용하여 수행되었으며 유전자 데이터에서 생성 된 처음 3 가지 주성분을 사용했습니다. 분석 결과는 다음과 같습니다.동일한 그래프에서 여러 줄의 주위에 신뢰 구간 그리기

distance=c(2.806063,8.208133,14.03604,19.03151,24.44091, 2.806063, 8.208133,14.03604,19.03151,24.44091,2.806063,8.208133,14.03604,19.03151,24.44091)  

correlation=c(-0.006933,0.029481,-0.071406,0.038319,-0.049990,0.006267,0.055945,-0.048551,-0.035062,-0.031578,0.022629,-0.065584,0.000986,-0.052754,0.0424931) 
component=c(PC1,PC1,PC1,PC1,PC1,PC2,PC2,PC2,PC2,PC2,PC3,PC3,PC3,PC3,PC3) 

data1<-data.frame(distance,correlation,component) 

난 후, 그 결과 지금은 주요 구성 요소 각각에 대한 95 % 신뢰 구간을 계산하면된다 좋아하고 ggplots에 그 그릴 것입니다 무엇

library(ggplot2) 
ggplot(data1,aes(x=data1$distance,y=data1$correlation,group=component,colour=component))+theme_classic()+ geom_line(size=1)+geom_point(size=1.5) 

음모 ggplot을 사용 희미한 음영을 사용하여 각 라인 주변의 신뢰 영역을 사용하고 서로 다른 PC를 나타내는 다양한 라인 색상을 유지합니다. 불행히도, 나는 완전히 붙어있어 이것을하는 법을 모릅니다. 어떤 도움도 감사 할거야.

+0

재현 가능한 예제를 포함하면 도움이됩니다. ('PC1','PC2' 등은 볼 수 없으므로'구성 요소 '를 만들 수 없습니다.) – jbaums

+0

구성 요소를 만들 필요가 없습니다. 그 변수는 단순히 제 2 변수에 의해 표현 된 상관 계수들을 그룹화하는 방법으로서 포함된다. – user3260086

+1

예를 들어 'component <- factor (rep (1 : 3, each = 5))'등의 이유는 무엇입니까? 코드가 실행되지 않으면 대답을 얻지 못할 가능성이 있습니다. – jbaums

답변

1

코드가 그대로 실행되지 않으므로 아무도 지난 10 시간 동안 응답하지 않았습니다.

가정 뜻 :

component=c("PC1","PC1","PC1","PC1","PC1","PC2","PC2","PC2","PC2","PC2","PC3","PC3","PC3","PC3","PC3") 

을하고 대 거리의 상관 관계에 대한 95 %의 CL을 원하는이 그것을 제공 할 것입니다 :

library(ggplot2) 
ggplot(data1,aes(x=distance,y=correlation,color=component))+ 
    geom_line(size=1)+ 
    geom_point(size=1.5)+ 
    stat_smooth(aes(fill=component), alpha=.2, 
       method=lm, formula=y~1, se=TRUE, level=0.95)+ 
    theme_classic() 

주요 추가 사항은 stat_smooth(...) 선이며 선형 모델을 사용하여 상관 관계 데이터와 거리 데이터를 매끄럽게합니다. 상수 항만 (즉, 평균)을 갖는다. 이 경우 기본값이 level=0.95이고 기본값은 se=TRUE이므로 이러한 절이 실제로 필요하지 않습니다.

또한, aes(...)에 대한 호출의 표현은 데이터 1의 열을 참조해야합니다 (그래서 x=distance하지x=data1$distance), 및 color=...가 같은 그룹 변수를 사용하는 경우에는 group=... 절을 필요로하지 않습니다.

+0

고마워, 그게 내가 필요한거야. 코드에 대한 다른 유용한 의견을 주셔서 감사합니다. 나는 R에 익숙하지 않고 코드 작성법을 배우기 시작했다. – user3260086