나는 시계열이 여러 개 있는데, 각 시계열은 식물 종을 나타냅니다. 우디 밀도에 따라 패턴이 있다고 생각합니다. 높은 우디 밀도 종은 비가 오는 기간에 꽃이 피 웁니다. 낮은 기간 동안 우디 밀도가 높은 종입니다.여러 시계열의 차이를 테스트하는 방법 R
많은 종의 시계열과 우디 밀도 측정치를 사용하여이 패턴을 보여주기 위해 이것을 R으로 모델링하는 방법은 무엇입니까?
#Woody Density
set.seed(69)
wden<-round(c(rnorm(5,mean=50),rnorm(5,mean=90)))
names(wden)<-c(paste("sp",1:10,sep=""))
wden
#Chuva
rain<-c(150,100,50,40,20,20,30,50,70,100,150,200,150,100,50,30,20,20,40,50,70,100,150,200)
#Flowering measures
ydet<-c(10,10,10,10,20,40,50,40,20,10,10,10)
#2 years for 5 low woody density and 5 high density species
flowering<-matrix(NA,nrow=24, ncol=10,dimnames=list(paste("month",1:24,sep=""),paste("sp",1:10,sep="")))
for (i in 1:5) {
flowering[,i]<-round(c(ydet+rnorm(12,mean=5,sd=5),ydet+rnorm(12,mean=5,sd=5)),digits=2)
}
for (i in 6:10) {
flowering[,i]<-round(c(rnorm(12,mean=30,sd=5),rnorm(12,mean=30,sd=5)),digits=2)
}
#Changing objects to Time series
flowering<-ts(flowering)
#Plot series
plot(flowering)
#Making colors for wood density
cores<-heat.colors(10,alpha=1)
matplot(c(1:24),flowering,type="l",lwd=2,lty=1,xlab="Time",ylab="Flowering",col=cores[order(wden)])
#Plotting Rain Together with time series
bargraph<-barplot(rain/max(rain)*100,xlab="Time",ylab="Rain")
matlines(bargraph,flowering,type="l",lwd=2,lty=1,xlab="Time",ylab="Flowering",col=cores[order(wden)])
axis(1,at=bargraph,labels=1:24)
axis(4,at=seq(0,100,by=10))
'brota'와'dmad' 객체를 정의하는 것을 잊었을 것 같습니다. – nograpes
작전, 코드 포르투갈어, 그때 더 많은 의미를 만들기 위해 일을 번역하려고했는데, 내가 모든 라인에서 이름을 변경하는 것을 잊지 것, 죄송합니다, 지금 바로 수정합니다 :) –
글쎄, 나는 모든 것이 작동한다고 생각합니다. 감사 –