jeudi 11 janvier 2018

stablitty else if in r

I have value of Kor and want to record the order of the value of Kor based on certain criteria and write its order of index

Kor<-c(0.91,0.58,0.86,1.76,0.67,0.55,1.83,1.69,2.21,1.19,1.61,1.36,1.66,1.31,0.63,2.10,1.01,0.81,0.78,0.48,0.82,0.59,1.88,0.81,1.33,1.65,3.60,0.59,0.55,0.61,0.72,2.07,1.10,1.29,0.58,0.98,1.41,2.14,1.08,1.65,0.80,0.71,0.94,1.46,1.51,2.57,0.57,1.21,1.11,0.62,0.98,1.90,6.30,0.68,2.50,2.75,2.05,1.36,1.12,0.75,0.87,0.46,1.53,0.85,2.58,2.49,3.29,3.18,6.11,2.15,2.19,1.35,0.68,1.10,2.16,1.39,2.86,1.68,1.04,0.63,0.52,1.23,0.62,2.34,2.39,2.71,2.33,4.70,3.81,1.16,0.76,0.54,0.70,1.44,1.04,0.59,3.82,29.04,8.72,1.76,0.97,0.40,2.56,3.22,1.40,0.86,2.21,1.77,2.23,0.74,0.88,0.84,1.31,1.29,0.59,1.89,4.99,0.84,1.11,1.84,0.99,0.64,0.76,0.57,0.76,2.93,0.87,7.32,1.72,2.73,2.12,0.89,1.56,2.42,2.49,4.73,8.97,5.72,13.14,17.28,6.50,1.60,1.99,2.76,39.56,6.25,26.56,2.79,1.16,0.83,1.74,0.42,1.12,0.53,1.40,1.61,2.08,0.75,0.58,0.89,1.27,1.10,2.03,3.03,5.74,4.25,2.20,2.51,3.62,0.95,10.94,2.33,5.54,1.44,0.88,2.13,6.95,6.67,11.89,8.11,81.49,12.54,379.47,102.71,71.31,214.85,80.17,Inf,43.88,10.43,2.00,75.05,256.82,179.64,150.71,167.95,10.17,11.67,17.63,32.48,15.05,16.52,8.93,9.41,14.77,2.43,1.67,3.45,1.98,3.79,4.37,3.10,8.29,13.41,11.02,15.77,127.00,25.18,121.84,58.99,13.38,20.55,37.52,12.15,4.40,76.13,33.45,17.45,234.88,57.35,54.64,43.40,98.77,58.00,155.62,21.75,36.71,24.07,1.84,3.94,89.79,111.31,45.55,60.36,65.40,139.66,44.38,226.70,37.83,35.94,26.32,78.31,108.83,36.15,17.47,23.30,115.86,19.60,38.11,72.62,144.56,21.33,4.62,17.28,27.82,45.76,17.77,18.48,21.82,30.63,8.21,12.97,16.70,28.97,25.55,29.35,238.34,7.86,4.48,12.07,12.89,33.19,4.62,16.47,21.61,33.33,28.44,157.04,52.54,72.30,71.84,78.05,165.81,70.56,139.12,47.53,193.11,4.02,5.72,3.42,2.03,2.19,2.51,1.47,4.37,5.25,9.63,2.64,3.37,2.98,2.14,3.39,1.58,0.65,1.18,1.71,5.61,2.11,1.38,1.90,0.89,1.23,2.00,8.57,3.44,17.32,8.05,5.53,5.37,0.81,2.27,1.88,2.04,7.56,7.43,0.92,0.95,0.80,0.63,0.88,1.58,1.45,1.59,1.42,0.83,0.69,2.08,1.10,0.61,0.70,0.56,0.83,1.15,1.32,2.37,2.13,0.91,0.51,0.80,4.59,2.43,1.47,1.50,2.39,1.87)

Then i make this code

Korr<-cbind(c(seq(1,365)),c(Kor[1:365]))
Korr<-as.data.frame(Korr)
id_0<-NULL
id_2<-NULL
id_5<-NULL
id_10<-NULL
id_20<-NULL
id_30<-NULL
for(i in 1:365){
    if(Korr[i,2]>=0 &Korr[i,2]<2) {id1_0=i}
    else if(Korr[i,2]>=2 & Korr[i,2]<5) {id1_2=i}
    else if(Korr[i,2]>=5 & Korr[i,2]<10) {id1_5=i}
    else if(Korr[i,2]>=10 & Korr[i,2]<20) {id1_10=i}
    else if(Korr[i,2]>=20 & Korr[i,2]<30) {id1_20=i}
    else {id1_30=i}
    id_0<-c(id_0,id1_0)
    id_2<-c(id_2,id1_2)
    id_5<-c(id_5,id1_5)
    id_10<-c(id_10,id1_10)
    id_20<-c(id_20,id1_20)
    id_30<-c(id_30,id1_30)
}

but it can't run. How to solve it

Aucun commentaire:

Enregistrer un commentaire