pr<-read.table(file="C:/VS-model/res638 result/ 1975 - 2013.dat", skip=1, header = TRUE)
View(pr)
pr<-read.table(file="C:/VS-model/res638 result/ rate 1975 - 2013.dat", skip=1, header = TRUE)
View(pr)
cli<-read.table(file="C:/VS-model/res638 result/ rate 1975 - 2013.dat", skip=1, header = TRUE)
View(pr)
cli<-read.table(file="C:/VS-model/res638 result/ rate 1975 - 2013.dat", header = T)
cli<-read.table(file="C:/VS-model/res638 result/ rate 1975 - 2013.dat", header = F)
cli<-read.table(file="C:/VS-model/res638 result/ rate 1975 - 2013.dat", skip=0, header = T)
cli<-read.table(file="C:/VS-model/res638 result/ rate 1975 - 2013.dat", skip=1, header = T)
View(cli)
cli<-read.table(file="C:/VS-model/res638 result/ rate 1975 - 2013.dat", skip=1, header = F)
View(cli)
cum<-unlist(lapply(split(cli,cli$year),function(x) round(cumsum(x[8]),4)))
View(cli)
cum<-unlist(lapply(split(cli,cli$year),function(x) cumsum(x[11])))
cum<-unlist(lapply(split(cli,cli$V1),function(x) cumsum(x[11])))
cum<-unlist(lapply(split(cli,cli$V1),function(x) mean(x[11])))
cum<-lapply(split(cli,cli$V1),function(x) mean(x[11]))
cum<-sapply(split(cli,cli$V1),function(x) mean(x[11]))
cum<-sapply(split(cli,cli$V1),function(x) sum(x[11]))
cum<-sapply(split(cli,cli$V1),function(x) mean(x[11]))
cum<-sapply(split(cli,cli$V1),function(x) sum(x[11])/length(x[11]))
cum<-sapply(split(cli,cli$V1),function(x) sum(x[11])/length(x[,11]))
cum<-lapply(split(cli,cli$V1),function(x) sum(x[11])/length(x[,11]))
## скорости роста по клеточным измерениям
getwd()
setwd("C:/Users/Рита/Desktop")
tab<-read.table(file="p_all.csv", sep=";", header=T) ## читаем таблицу с клеточными измерениями
tab_new<-1/subset(tab, select = -c(1:1)) ## удаляем ненужные года и берем обратное число от каждого элемента
pr<-read.table(file="C:/VS-model/486 result/sira 1970 - 2008.dat", skip=1, header = TRUE)
beg<-cbind(pr$year,pr$BG1)
end<-cbind(pr$year,pr$EG1)
period<-end[,2]-beg[,2] ## считаем период роста
str(period)
ring<-colSums(tab)
cumsum<-lapply(tab,function(x) cumsum(x))
t<- do.call(cbind.data.frame, cumsum)
tt<-t
for(i in 1:nrow(t)){t[i,]<-round(t[i,]*period/ring)}
cum1<-1/(16*t[1,])
cum2<-lapply(t,function(x) diff(x))
cum<-lapply(t,function(x) 1/(16*diff(x)))
rate<-do.call(cbind.data.frame, cum)
rrate<-rbind(cum1,rate)
write.table(round(rrate,5), file = "C:/Users/Рита/Desktop/rate.csv", sep = ";", col.names=TRUE, row.names=FALSE)
## день формирования каждой клетки
day<-beg[,2]
dday<- t(matrix(day,nrow=39,ncol=16))
sum<-t+dday
write.table(sum, file = "C:/Users/Рита/Desktop/day.csv", sep = ",", col.names=F, row.names=FALSE)
## скорости роста по модели
setwd("C:/Users/Рита/Desktop")
cli<-read.table(file="C:/VS-model/486 result/sira rate 1970 - 2008.dat", skip=1, header = F)
pr<-read.table(file="C:/VS-model/486 result/sira 1970 - 2008.dat", skip=1, header = TRUE)
t<-t-rbind(0,t[-c(16),])
t0<-rbind(0,cumsum(t[-c(16),]))
t1<-t
t2<-list(1:16)
t3<-list(1:16)
t4<-list(1:16)
t5<-list(1:16)
t6<-list(1:16)
t7<-list(1:16)
t8<-list(1:16)
t9<-list(1:16)
t10<-list(1:16)
t11<-list(1:16)
t12<-list(1:16)
t13<-list(1:16)
t14<-list(1:16)
t15<-list(1:16)
t16<-list(1:16)
t17<-list(1:16)
t18<-list(1:16)
t19<-list(1:16)
t20<-list(1:16)
t21<-list(1:16)
t22<-list(1:16)
t23<-list(1:16)
t24<-list(1:16)
t25<-list(1:16)
t26<-list(1:16)
t27<-list(1:16)
t28<-list(1:16)
t29<-list(1:16)
t30<-list(1:16)
t31<-list(1:16)
t32<-list(1:16)
t33<-list(1:16)
t34<-list(1:16)
t35<-list(1:16)
t36<-list(1:16)
t37<-list(1:16)
t38<-list(1:16)
t39<-list(1:16)
t40<-list(1:16)
scli<-split(cli, cli$V1)
names(scli) <- c(1:39)
for (i in 1:16){
t2[i]<-sum(scli$"1"[(beg[[1,2]]+t0[i,1]):((beg[[1,2]]+t0[i,1])+(t1[i,1]-1)),10])/t[i,1]
t3[i]<-sum(scli$"2"[(beg[[2,2]]+t0[i,2]):((beg[[2,2]]+t0[i,2])+(t1[i,2]-1)),10])/t[i,2]
t4[i]<-sum(scli$"3"[(beg[[3,2]]+t0[i,3]):((beg[[3,2]]+t0[i,3])+(t1[i,3]-1)),10])/t[i,3]
t5[i]<-sum(scli$"4"[(beg[[4,2]]+t0[i,4]):((beg[[4,2]]+t0[i,4])+(t1[i,4]-1)),10])/t[i,4]
t6[i]<-sum(scli$"5"[(beg[[5,2]]+t0[i,5]):((beg[[5,2]]+t0[i,5])+(t1[i,5]-1)),10])/t[i,5]
t7[i]<-sum(scli$"6"[(beg[[6,2]]+t0[i,6]):((beg[[6,2]]+t0[i,6])+(t1[i,6]-1)),10])/t[i,6]
t8[i]<-sum(scli$"7"[(beg[[7,2]]+t0[i,7]):((beg[[7,2]]+t0[i,7])+(t1[i,7]-1)),10])/t[i,7]
t9[i]<-sum(scli$"8"[(beg[[8,2]]+t0[i,8]):((beg[[8,2]]+t0[i,8])+(t1[i,8]-1)),10])/t[i,8]
t10[i]<-sum(scli$"9"[(beg[[9,2]]+t0[i,9]):((beg[[9,2]]+t0[i,9])+(t1[i,9]-1)),10])/t[i,9]
t11[i]<-sum(scli$"10"[(beg[[10,2]]+t0[i,10]):((beg[[10,2]]+t0[i,10])+(t1[i,10]-1)),10])/t[i,10]
t12[i]<-sum(scli$"11"[(beg[[11,2]]+t0[i,11]):((beg[[11,2]]+t0[i,11])+(t1[i,11]-1)),10])/t[i,11]
t13[i]<-sum(scli$"12"[(beg[[12,2]]+t0[i,12]):((beg[[12,2]]+t0[i,12])+(t1[i,12]-1)),10])/t[i,12]
t14[i]<-sum(scli$"13"[(beg[[13,2]]+t0[i,13]):((beg[[13,2]]+t0[i,13])+(t1[i,13]-1)),10])/t[i,13]
t15[i]<-sum(scli$"14"[(beg[[14,2]]+t0[i,14]):((beg[[14,2]]+t0[i,14])+(t1[i,14]-1)),10])/t[i,14]
t16[i]<-sum(scli$"15"[(beg[[15,2]]+t0[i,15]):((beg[[15,2]]+t0[i,15])+(t1[i,15]-1)),10])/t[i,15]
t17[i]<-sum(scli$"16"[(beg[[16,2]]+t0[i,16]):((beg[[16,2]]+t0[i,16])+(t1[i,16]-1)),10])/t[i,16]
t18[i]<-sum(scli$"17"[(beg[[17,2]]+t0[i,17]):((beg[[17,2]]+t0[i,17])+(t1[i,17]-1)),10])/t[i,17]
t19[i]<-sum(scli$"18"[(beg[[18,2]]+t0[i,18]):((beg[[18,2]]+t0[i,18])+(t1[i,18]-1)),10])/t[i,18]
t20[i]<-sum(scli$"19"[(beg[[19,2]]+t0[i,19]):((beg[[19,2]]+t0[i,19])+(t1[i,19]-1)),10])/t[i,19]
t21[i]<-sum(scli$"20"[(beg[[20,2]]+t0[i,20]):((beg[[20,2]]+t0[i,20])+(t1[i,20]-1)),10])/t[i,20]
t22[i]<-sum(scli$"21"[(beg[[21,2]]+t0[i,21]):((beg[[21,2]]+t0[i,21])+(t1[i,21]-1)),10])/t[i,21]
t23[i]<-sum(scli$"22"[(beg[[22,2]]+t0[i,22]):((beg[[22,2]]+t0[i,22])+(t1[i,22]-1)),10])/t[i,22]
t24[i]<-sum(scli$"23"[(beg[[23,2]]+t0[i,23]):((beg[[23,2]]+t0[i,23])+(t1[i,23]-1)),10])/t[i,23]
t25[i]<-sum(scli$"24"[(beg[[24,2]]+t0[i,24]):((beg[[24,2]]+t0[i,24])+(t1[i,24]-1)),10])/t[i,24]
t26[i]<-sum(scli$"25"[(beg[[25,2]]+t0[i,25]):((beg[[25,2]]+t0[i,25])+(t1[i,25]-1)),10])/t[i,25]
t27[i]<-sum(scli$"26"[(beg[[26,2]]+t0[i,26]):((beg[[26,2]]+t0[i,26])+(t1[i,26]-1)),10])/t[i,26]
t28[i]<-sum(scli$"27"[(beg[[27,2]]+t0[i,27]):((beg[[27,2]]+t0[i,27])+(t1[i,27]-1)),10])/t[i,27]
t29[i]<-sum(scli$"28"[(beg[[28,2]]+t0[i,28]):((beg[[28,2]]+t0[i,28])+(t1[i,28]-1)),10])/t[i,28]
t30[i]<-sum(scli$"29"[(beg[[29,2]]+t0[i,29]):((beg[[29,2]]+t0[i,29])+(t1[i,29]-1)),10])/t[i,29]
t31[i]<-sum(scli$"30"[(beg[[30,2]]+t0[i,30]):((beg[[30,2]]+t0[i,30])+(t1[i,30]-1)),10])/t[i,30]
t32[i]<-sum(scli$"31"[(beg[[31,2]]+t0[i,31]):((beg[[31,2]]+t0[i,31])+(t1[i,31]-1)),10])/t[i,31]
t33[i]<-sum(scli$"32"[(beg[[32,2]]+t0[i,32]):((beg[[32,2]]+t0[i,32])+(t1[i,32]-1)),10])/t[i,32]
t34[i]<-sum(scli$"33"[(beg[[33,2]]+t0[i,33]):((beg[[33,2]]+t0[i,33])+(t1[i,33]-1)),10])/t[i,33]
t35[i]<-sum(scli$"34"[(beg[[34,2]]+t0[i,34]):((beg[[34,2]]+t0[i,34])+(t1[i,34]-1)),10])/t[i,34]
t36[i]<-sum(scli$"35"[(beg[[35,2]]+t0[i,35]):((beg[[35,2]]+t0[i,35])+(t1[i,35]-1)),10])/t[i,35]
t37[i]<-sum(scli$"36"[(beg[[36,2]]+t0[i,36]):((beg[[36,2]]+t0[i,36])+(t1[i,36]-1)),10])/t[i,36]
t38[i]<-sum(scli$"37"[(beg[[37,2]]+t0[i,37]):((beg[[37,2]]+t0[i,37])+(t1[i,37]-1)),10])/t[i,37]
t39[i]<-sum(scli$"38"[(beg[[38,2]]+t0[i,38]):((beg[[38,2]]+t0[i,38])+(t1[i,38]-1)),10])/t[i,38]
t40[i]<-sum(scli$"39"[(beg[[39,2]]+t0[i,39]):((beg[[39,2]]+t0[i,39])+(t1[i,39]-1)),10])/t[i,39]
}
mrate<-t(do.call(rbind.data.frame, list(t2,t3,t4,t5,t6,t7,t8,t9,t10,t11,t12,t13,t14,t15,t16,t17,t18,t19,t20,t21,t22,t23,t24,t25,t26,t27,t28,t29,t30,t31,t32,t33,t34,t35,t36,t37,t38,t39,t40)))
colnames(mrate) <- unique(cli[,1])
rownames(mrate) <- c(1:16)
write.table(mrate, file = "C:/Users/Рита/Desktop/mrate.csv", sep = ",", col.names=TRUE, row.names=FALSE)
##преобразование модельных скоростей с помощью метода натменьших квадратов
y<-rrate/mrate
write.table(y, file = "C:/Users/Рита/Desktop/y.csv", sep = ",", col.names=TRUE, row.names=FALSE)
y<-do.call(data.frame,lapply(y, function(x) replace(x, is.infinite(x),0)))
y<-do.call(data.frame,lapply(y, function(x) replace(x, is.nan(x),0)))
x<-matrix(rep(1:16,each=39),nrow=16,ncol=39,byrow=TRUE)
write.table(x, file = "C:/Users/Рита/Desktop/x.csv", sep = ",", col.names=TRUE, row.names=FALSE)
xy <- x*y
xm<-mean(colMeans(x, na.rm =  TRUE))
ym<-mean(colMeans(y, na.rm =  TRUE))
xym<-mean(colMeans(xy, na.rm =  TRUE))
x2m<-mean(colMeans((x^2), na.rm =  TRUE))
b<-(xym-xm*ym)/(x2m-xm^2)
a<-ym-b*xm
v0<-b*x+a
v<-mrate*v0
write.table(round(v,5), file = "C:/Users/Рита/Desktop/rate_mm.csv", sep = ";", col.names=TRUE, row.names=FALSE)
View(v)
View(t)
View(pr)
View(rate)
View(rrate)
## скорости роста по клеточным измерениям
getwd()
setwd("C:/Users/Рита/Desktop")
tab<-read.table(file="p_all.csv", sep=";", header=T) ## читаем таблицу с клеточными измерениями
tab_new<-1/subset(tab, select = -c(1:1)) ## удаляем ненужные года и берем обратное число от каждого элемента
pr<-read.table(file="C:/VS-model/456 result/sira 1970 - 2008.dat", skip=1, header = TRUE)
beg<-cbind(pr$year,pr$BG1)
end<-cbind(pr$year,pr$EG1)
period<-end[,2]-beg[,2] ## считаем период роста
str(period)
ring<-colSums(tab)
cumsum<-lapply(tab,function(x) cumsum(x))
t<- do.call(cbind.data.frame, cumsum)
tt<-t
for(i in 1:nrow(t)){t[i,]<-round(t[i,]*period/ring)}
cum1<-1/(16*t[1,])
cum2<-lapply(t,function(x) diff(x))
cum<-lapply(t,function(x) 1/(16*diff(x)))
rate<-do.call(cbind.data.frame, cum)
rrate<-rbind(cum1,rate)
write.table(round(rrate,5), file = "C:/Users/Рита/Desktop/rate.csv", sep = ";", col.names=TRUE, row.names=FALSE)
## день формирования каждой клетки
day<-beg[,2]
dday<- t(matrix(day,nrow=39,ncol=16))
sum<-t+dday
write.table(sum, file = "C:/Users/Рита/Desktop/day.csv", sep = ",", col.names=F, row.names=FALSE)
## скорости роста по модели
setwd("C:/Users/Рита/Desktop")
cli<-read.table(file="C:/VS-model/456 result/sira rate 1970 - 2008.dat", skip=1, header = F)
pr<-read.table(file="C:/VS-model/456 result/sira 1970 - 2008.dat", skip=1, header = TRUE)
t<-t-rbind(0,t[-c(16),])
t0<-rbind(0,cumsum(t[-c(16),]))
t1<-t
t2<-list(1:16)
t3<-list(1:16)
t4<-list(1:16)
t5<-list(1:16)
t6<-list(1:16)
t7<-list(1:16)
t8<-list(1:16)
t9<-list(1:16)
t10<-list(1:16)
t11<-list(1:16)
t12<-list(1:16)
t13<-list(1:16)
t14<-list(1:16)
t15<-list(1:16)
t16<-list(1:16)
t17<-list(1:16)
t18<-list(1:16)
t19<-list(1:16)
t20<-list(1:16)
t21<-list(1:16)
t22<-list(1:16)
t23<-list(1:16)
t24<-list(1:16)
t25<-list(1:16)
t26<-list(1:16)
t27<-list(1:16)
t28<-list(1:16)
t29<-list(1:16)
t30<-list(1:16)
t31<-list(1:16)
t32<-list(1:16)
t33<-list(1:16)
t34<-list(1:16)
t35<-list(1:16)
t36<-list(1:16)
t37<-list(1:16)
t38<-list(1:16)
t39<-list(1:16)
t40<-list(1:16)
scli<-split(cli, cli$V1)
names(scli) <- c(1:39)
for (i in 1:16){
t2[i]<-sum(scli$"1"[(beg[[1,2]]+t0[i,1]):((beg[[1,2]]+t0[i,1])+(t1[i,1]-1)),10])/t[i,1]
t3[i]<-sum(scli$"2"[(beg[[2,2]]+t0[i,2]):((beg[[2,2]]+t0[i,2])+(t1[i,2]-1)),10])/t[i,2]
t4[i]<-sum(scli$"3"[(beg[[3,2]]+t0[i,3]):((beg[[3,2]]+t0[i,3])+(t1[i,3]-1)),10])/t[i,3]
t5[i]<-sum(scli$"4"[(beg[[4,2]]+t0[i,4]):((beg[[4,2]]+t0[i,4])+(t1[i,4]-1)),10])/t[i,4]
t6[i]<-sum(scli$"5"[(beg[[5,2]]+t0[i,5]):((beg[[5,2]]+t0[i,5])+(t1[i,5]-1)),10])/t[i,5]
t7[i]<-sum(scli$"6"[(beg[[6,2]]+t0[i,6]):((beg[[6,2]]+t0[i,6])+(t1[i,6]-1)),10])/t[i,6]
t8[i]<-sum(scli$"7"[(beg[[7,2]]+t0[i,7]):((beg[[7,2]]+t0[i,7])+(t1[i,7]-1)),10])/t[i,7]
t9[i]<-sum(scli$"8"[(beg[[8,2]]+t0[i,8]):((beg[[8,2]]+t0[i,8])+(t1[i,8]-1)),10])/t[i,8]
t10[i]<-sum(scli$"9"[(beg[[9,2]]+t0[i,9]):((beg[[9,2]]+t0[i,9])+(t1[i,9]-1)),10])/t[i,9]
t11[i]<-sum(scli$"10"[(beg[[10,2]]+t0[i,10]):((beg[[10,2]]+t0[i,10])+(t1[i,10]-1)),10])/t[i,10]
t12[i]<-sum(scli$"11"[(beg[[11,2]]+t0[i,11]):((beg[[11,2]]+t0[i,11])+(t1[i,11]-1)),10])/t[i,11]
t13[i]<-sum(scli$"12"[(beg[[12,2]]+t0[i,12]):((beg[[12,2]]+t0[i,12])+(t1[i,12]-1)),10])/t[i,12]
t14[i]<-sum(scli$"13"[(beg[[13,2]]+t0[i,13]):((beg[[13,2]]+t0[i,13])+(t1[i,13]-1)),10])/t[i,13]
t15[i]<-sum(scli$"14"[(beg[[14,2]]+t0[i,14]):((beg[[14,2]]+t0[i,14])+(t1[i,14]-1)),10])/t[i,14]
t16[i]<-sum(scli$"15"[(beg[[15,2]]+t0[i,15]):((beg[[15,2]]+t0[i,15])+(t1[i,15]-1)),10])/t[i,15]
t17[i]<-sum(scli$"16"[(beg[[16,2]]+t0[i,16]):((beg[[16,2]]+t0[i,16])+(t1[i,16]-1)),10])/t[i,16]
t18[i]<-sum(scli$"17"[(beg[[17,2]]+t0[i,17]):((beg[[17,2]]+t0[i,17])+(t1[i,17]-1)),10])/t[i,17]
t19[i]<-sum(scli$"18"[(beg[[18,2]]+t0[i,18]):((beg[[18,2]]+t0[i,18])+(t1[i,18]-1)),10])/t[i,18]
t20[i]<-sum(scli$"19"[(beg[[19,2]]+t0[i,19]):((beg[[19,2]]+t0[i,19])+(t1[i,19]-1)),10])/t[i,19]
t21[i]<-sum(scli$"20"[(beg[[20,2]]+t0[i,20]):((beg[[20,2]]+t0[i,20])+(t1[i,20]-1)),10])/t[i,20]
t22[i]<-sum(scli$"21"[(beg[[21,2]]+t0[i,21]):((beg[[21,2]]+t0[i,21])+(t1[i,21]-1)),10])/t[i,21]
t23[i]<-sum(scli$"22"[(beg[[22,2]]+t0[i,22]):((beg[[22,2]]+t0[i,22])+(t1[i,22]-1)),10])/t[i,22]
t24[i]<-sum(scli$"23"[(beg[[23,2]]+t0[i,23]):((beg[[23,2]]+t0[i,23])+(t1[i,23]-1)),10])/t[i,23]
t25[i]<-sum(scli$"24"[(beg[[24,2]]+t0[i,24]):((beg[[24,2]]+t0[i,24])+(t1[i,24]-1)),10])/t[i,24]
t26[i]<-sum(scli$"25"[(beg[[25,2]]+t0[i,25]):((beg[[25,2]]+t0[i,25])+(t1[i,25]-1)),10])/t[i,25]
t27[i]<-sum(scli$"26"[(beg[[26,2]]+t0[i,26]):((beg[[26,2]]+t0[i,26])+(t1[i,26]-1)),10])/t[i,26]
t28[i]<-sum(scli$"27"[(beg[[27,2]]+t0[i,27]):((beg[[27,2]]+t0[i,27])+(t1[i,27]-1)),10])/t[i,27]
t29[i]<-sum(scli$"28"[(beg[[28,2]]+t0[i,28]):((beg[[28,2]]+t0[i,28])+(t1[i,28]-1)),10])/t[i,28]
t30[i]<-sum(scli$"29"[(beg[[29,2]]+t0[i,29]):((beg[[29,2]]+t0[i,29])+(t1[i,29]-1)),10])/t[i,29]
t31[i]<-sum(scli$"30"[(beg[[30,2]]+t0[i,30]):((beg[[30,2]]+t0[i,30])+(t1[i,30]-1)),10])/t[i,30]
t32[i]<-sum(scli$"31"[(beg[[31,2]]+t0[i,31]):((beg[[31,2]]+t0[i,31])+(t1[i,31]-1)),10])/t[i,31]
t33[i]<-sum(scli$"32"[(beg[[32,2]]+t0[i,32]):((beg[[32,2]]+t0[i,32])+(t1[i,32]-1)),10])/t[i,32]
t34[i]<-sum(scli$"33"[(beg[[33,2]]+t0[i,33]):((beg[[33,2]]+t0[i,33])+(t1[i,33]-1)),10])/t[i,33]
t35[i]<-sum(scli$"34"[(beg[[34,2]]+t0[i,34]):((beg[[34,2]]+t0[i,34])+(t1[i,34]-1)),10])/t[i,34]
t36[i]<-sum(scli$"35"[(beg[[35,2]]+t0[i,35]):((beg[[35,2]]+t0[i,35])+(t1[i,35]-1)),10])/t[i,35]
t37[i]<-sum(scli$"36"[(beg[[36,2]]+t0[i,36]):((beg[[36,2]]+t0[i,36])+(t1[i,36]-1)),10])/t[i,36]
t38[i]<-sum(scli$"37"[(beg[[37,2]]+t0[i,37]):((beg[[37,2]]+t0[i,37])+(t1[i,37]-1)),10])/t[i,37]
t39[i]<-sum(scli$"38"[(beg[[38,2]]+t0[i,38]):((beg[[38,2]]+t0[i,38])+(t1[i,38]-1)),10])/t[i,38]
t40[i]<-sum(scli$"39"[(beg[[39,2]]+t0[i,39]):((beg[[39,2]]+t0[i,39])+(t1[i,39]-1)),10])/t[i,39]
}
mrate<-t(do.call(rbind.data.frame, list(t2,t3,t4,t5,t6,t7,t8,t9,t10,t11,t12,t13,t14,t15,t16,t17,t18,t19,t20,t21,t22,t23,t24,t25,t26,t27,t28,t29,t30,t31,t32,t33,t34,t35,t36,t37,t38,t39,t40)))
colnames(mrate) <- unique(cli[,1])
rownames(mrate) <- c(1:16)
write.table(mrate, file = "C:/Users/Рита/Desktop/mrate.csv", sep = ",", col.names=TRUE, row.names=FALSE)
##преобразование модельных скоростей с помощью метода натменьших квадратов
y<-rrate/mrate
write.table(y, file = "C:/Users/Рита/Desktop/y.csv", sep = ",", col.names=TRUE, row.names=FALSE)
y<-do.call(data.frame,lapply(y, function(x) replace(x, is.infinite(x),0)))
y<-do.call(data.frame,lapply(y, function(x) replace(x, is.nan(x),0)))
x<-matrix(rep(1:16,each=39),nrow=16,ncol=39,byrow=TRUE)
write.table(x, file = "C:/Users/Рита/Desktop/x.csv", sep = ",", col.names=TRUE, row.names=FALSE)
xy <- x*y
xm<-mean(colMeans(x, na.rm =  TRUE))
ym<-mean(colMeans(y, na.rm =  TRUE))
xym<-mean(colMeans(xy, na.rm =  TRUE))
x2m<-mean(colMeans((x^2), na.rm =  TRUE))
b<-(xym-xm*ym)/(x2m-xm^2)
a<-ym-b*xm
v0<-b*x+a
v<-mrate*v0
write.table(round(v,5), file = "C:/Users/Рита/Desktop/rate_mm.csv", sep = ";", col.names=TRUE, row.names=FALSE)
View(v)
getwd()
setwd("C:/VS-model/jakutsk_clim")
temp <- list.files(pattern="*.cli")
myfiles <- lapply(temp, read.delim, sep=" ", header=F)
df <- do.call(rbind.data.frame, myfiles)
write.table(df, file = "C:/Users/Рита/Desktop/all_years.csv", sep = ";", col.names=FALSE, row.names=FALSE)
## скорости роста по клеточным измерениям
getwd()
setwd("C:/Users/Рита/Desktop")
tab<-read.table(file="p_all.csv", sep=";", header=T) ## читаем таблицу с клеточными измерениями
tab_new<-1/subset(tab, select = -c(1:1)) ## удаляем ненужные года и берем обратное число от каждого элемента
pr<-read.table(file="C:/VS-model/463 result/sira 1970 - 2008.dat", skip=1, header = TRUE)
beg<-cbind(pr$year,pr$BG1)
end<-cbind(pr$year,pr$EG1)
period<-end[,2]-beg[,2] ## считаем период роста
str(period)
ring<-colSums(tab)
cumsum<-lapply(tab,function(x) cumsum(x))
t<- do.call(cbind.data.frame, cumsum)
tt<-t
for(i in 1:nrow(t)){t[i,]<-round(t[i,]*period/ring)}
cum1<-1/(16*t[1,])
cum2<-lapply(t,function(x) diff(x))
cum<-lapply(t,function(x) 1/(16*diff(x)))
rate<-do.call(cbind.data.frame, cum)
rrate<-rbind(cum1,rate)
write.table(round(rrate,5), file = "C:/Users/Рита/Desktop/rate.csv", sep = ";", col.names=TRUE, row.names=FALSE)
## день формирования каждой клетки
day<-beg[,2]
dday<- t(matrix(day,nrow=39,ncol=16))
sum<-t+dday
write.table(sum, file = "C:/Users/Рита/Desktop/day.csv", sep = ",", col.names=F, row.names=FALSE)
## скорости роста по модели
setwd("C:/Users/Рита/Desktop")
cli<-read.table(file="C:/VS-model/463 result/sira rate 1970 - 2008.dat", skip=1, header = F)
pr<-read.table(file="C:/VS-model/463 result/sira 1970 - 2008.dat", skip=1, header = TRUE)
t<-t-rbind(0,t[-c(16),])
t0<-rbind(0,cumsum(t[-c(16),]))
t1<-t
t2<-list(1:16)
t3<-list(1:16)
t4<-list(1:16)
t5<-list(1:16)
t6<-list(1:16)
t7<-list(1:16)
t8<-list(1:16)
t9<-list(1:16)
t10<-list(1:16)
t11<-list(1:16)
t12<-list(1:16)
t13<-list(1:16)
t14<-list(1:16)
t15<-list(1:16)
t16<-list(1:16)
t17<-list(1:16)
t18<-list(1:16)
t19<-list(1:16)
t20<-list(1:16)
t21<-list(1:16)
t22<-list(1:16)
t23<-list(1:16)
t24<-list(1:16)
t25<-list(1:16)
t26<-list(1:16)
t27<-list(1:16)
t28<-list(1:16)
t29<-list(1:16)
t30<-list(1:16)
t31<-list(1:16)
t32<-list(1:16)
t33<-list(1:16)
t34<-list(1:16)
t35<-list(1:16)
t36<-list(1:16)
t37<-list(1:16)
t38<-list(1:16)
t39<-list(1:16)
t40<-list(1:16)
scli<-split(cli, cli$V1)
names(scli) <- c(1:39)
for (i in 1:16){
t2[i]<-sum(scli$"1"[(beg[[1,2]]+t0[i,1]):((beg[[1,2]]+t0[i,1])+(t1[i,1]-1)),10])/t[i,1]
t3[i]<-sum(scli$"2"[(beg[[2,2]]+t0[i,2]):((beg[[2,2]]+t0[i,2])+(t1[i,2]-1)),10])/t[i,2]
t4[i]<-sum(scli$"3"[(beg[[3,2]]+t0[i,3]):((beg[[3,2]]+t0[i,3])+(t1[i,3]-1)),10])/t[i,3]
t5[i]<-sum(scli$"4"[(beg[[4,2]]+t0[i,4]):((beg[[4,2]]+t0[i,4])+(t1[i,4]-1)),10])/t[i,4]
t6[i]<-sum(scli$"5"[(beg[[5,2]]+t0[i,5]):((beg[[5,2]]+t0[i,5])+(t1[i,5]-1)),10])/t[i,5]
t7[i]<-sum(scli$"6"[(beg[[6,2]]+t0[i,6]):((beg[[6,2]]+t0[i,6])+(t1[i,6]-1)),10])/t[i,6]
t8[i]<-sum(scli$"7"[(beg[[7,2]]+t0[i,7]):((beg[[7,2]]+t0[i,7])+(t1[i,7]-1)),10])/t[i,7]
t9[i]<-sum(scli$"8"[(beg[[8,2]]+t0[i,8]):((beg[[8,2]]+t0[i,8])+(t1[i,8]-1)),10])/t[i,8]
t10[i]<-sum(scli$"9"[(beg[[9,2]]+t0[i,9]):((beg[[9,2]]+t0[i,9])+(t1[i,9]-1)),10])/t[i,9]
t11[i]<-sum(scli$"10"[(beg[[10,2]]+t0[i,10]):((beg[[10,2]]+t0[i,10])+(t1[i,10]-1)),10])/t[i,10]
t12[i]<-sum(scli$"11"[(beg[[11,2]]+t0[i,11]):((beg[[11,2]]+t0[i,11])+(t1[i,11]-1)),10])/t[i,11]
t13[i]<-sum(scli$"12"[(beg[[12,2]]+t0[i,12]):((beg[[12,2]]+t0[i,12])+(t1[i,12]-1)),10])/t[i,12]
t14[i]<-sum(scli$"13"[(beg[[13,2]]+t0[i,13]):((beg[[13,2]]+t0[i,13])+(t1[i,13]-1)),10])/t[i,13]
t15[i]<-sum(scli$"14"[(beg[[14,2]]+t0[i,14]):((beg[[14,2]]+t0[i,14])+(t1[i,14]-1)),10])/t[i,14]
t16[i]<-sum(scli$"15"[(beg[[15,2]]+t0[i,15]):((beg[[15,2]]+t0[i,15])+(t1[i,15]-1)),10])/t[i,15]
t17[i]<-sum(scli$"16"[(beg[[16,2]]+t0[i,16]):((beg[[16,2]]+t0[i,16])+(t1[i,16]-1)),10])/t[i,16]
t18[i]<-sum(scli$"17"[(beg[[17,2]]+t0[i,17]):((beg[[17,2]]+t0[i,17])+(t1[i,17]-1)),10])/t[i,17]
t19[i]<-sum(scli$"18"[(beg[[18,2]]+t0[i,18]):((beg[[18,2]]+t0[i,18])+(t1[i,18]-1)),10])/t[i,18]
t20[i]<-sum(scli$"19"[(beg[[19,2]]+t0[i,19]):((beg[[19,2]]+t0[i,19])+(t1[i,19]-1)),10])/t[i,19]
t21[i]<-sum(scli$"20"[(beg[[20,2]]+t0[i,20]):((beg[[20,2]]+t0[i,20])+(t1[i,20]-1)),10])/t[i,20]
t22[i]<-sum(scli$"21"[(beg[[21,2]]+t0[i,21]):((beg[[21,2]]+t0[i,21])+(t1[i,21]-1)),10])/t[i,21]
t23[i]<-sum(scli$"22"[(beg[[22,2]]+t0[i,22]):((beg[[22,2]]+t0[i,22])+(t1[i,22]-1)),10])/t[i,22]
t24[i]<-sum(scli$"23"[(beg[[23,2]]+t0[i,23]):((beg[[23,2]]+t0[i,23])+(t1[i,23]-1)),10])/t[i,23]
t25[i]<-sum(scli$"24"[(beg[[24,2]]+t0[i,24]):((beg[[24,2]]+t0[i,24])+(t1[i,24]-1)),10])/t[i,24]
t26[i]<-sum(scli$"25"[(beg[[25,2]]+t0[i,25]):((beg[[25,2]]+t0[i,25])+(t1[i,25]-1)),10])/t[i,25]
t27[i]<-sum(scli$"26"[(beg[[26,2]]+t0[i,26]):((beg[[26,2]]+t0[i,26])+(t1[i,26]-1)),10])/t[i,26]
t28[i]<-sum(scli$"27"[(beg[[27,2]]+t0[i,27]):((beg[[27,2]]+t0[i,27])+(t1[i,27]-1)),10])/t[i,27]
t29[i]<-sum(scli$"28"[(beg[[28,2]]+t0[i,28]):((beg[[28,2]]+t0[i,28])+(t1[i,28]-1)),10])/t[i,28]
t30[i]<-sum(scli$"29"[(beg[[29,2]]+t0[i,29]):((beg[[29,2]]+t0[i,29])+(t1[i,29]-1)),10])/t[i,29]
t31[i]<-sum(scli$"30"[(beg[[30,2]]+t0[i,30]):((beg[[30,2]]+t0[i,30])+(t1[i,30]-1)),10])/t[i,30]
t32[i]<-sum(scli$"31"[(beg[[31,2]]+t0[i,31]):((beg[[31,2]]+t0[i,31])+(t1[i,31]-1)),10])/t[i,31]
t33[i]<-sum(scli$"32"[(beg[[32,2]]+t0[i,32]):((beg[[32,2]]+t0[i,32])+(t1[i,32]-1)),10])/t[i,32]
t34[i]<-sum(scli$"33"[(beg[[33,2]]+t0[i,33]):((beg[[33,2]]+t0[i,33])+(t1[i,33]-1)),10])/t[i,33]
t35[i]<-sum(scli$"34"[(beg[[34,2]]+t0[i,34]):((beg[[34,2]]+t0[i,34])+(t1[i,34]-1)),10])/t[i,34]
t36[i]<-sum(scli$"35"[(beg[[35,2]]+t0[i,35]):((beg[[35,2]]+t0[i,35])+(t1[i,35]-1)),10])/t[i,35]
t37[i]<-sum(scli$"36"[(beg[[36,2]]+t0[i,36]):((beg[[36,2]]+t0[i,36])+(t1[i,36]-1)),10])/t[i,36]
t38[i]<-sum(scli$"37"[(beg[[37,2]]+t0[i,37]):((beg[[37,2]]+t0[i,37])+(t1[i,37]-1)),10])/t[i,37]
t39[i]<-sum(scli$"38"[(beg[[38,2]]+t0[i,38]):((beg[[38,2]]+t0[i,38])+(t1[i,38]-1)),10])/t[i,38]
t40[i]<-sum(scli$"39"[(beg[[39,2]]+t0[i,39]):((beg[[39,2]]+t0[i,39])+(t1[i,39]-1)),10])/t[i,39]
}
mrate<-t(do.call(rbind.data.frame, list(t2,t3,t4,t5,t6,t7,t8,t9,t10,t11,t12,t13,t14,t15,t16,t17,t18,t19,t20,t21,t22,t23,t24,t25,t26,t27,t28,t29,t30,t31,t32,t33,t34,t35,t36,t37,t38,t39,t40)))
colnames(mrate) <- unique(cli[,1])
rownames(mrate) <- c(1:16)
write.table(mrate, file = "C:/Users/Рита/Desktop/mrate.csv", sep = ",", col.names=TRUE, row.names=FALSE)
##преобразование модельных скоростей с помощью метода натменьших квадратов
y<-rrate/mrate
write.table(y, file = "C:/Users/Рита/Desktop/y.csv", sep = ",", col.names=TRUE, row.names=FALSE)
y<-do.call(data.frame,lapply(y, function(x) replace(x, is.infinite(x),0)))
y<-do.call(data.frame,lapply(y, function(x) replace(x, is.nan(x),0)))
x<-matrix(rep(1:16,each=39),nrow=16,ncol=39,byrow=TRUE)
write.table(x, file = "C:/Users/Рита/Desktop/x.csv", sep = ",", col.names=TRUE, row.names=FALSE)
xy <- x*y
xm<-mean(colMeans(x, na.rm =  TRUE))
ym<-mean(colMeans(y, na.rm =  TRUE))
xym<-mean(colMeans(xy, na.rm =  TRUE))
x2m<-mean(colMeans((x^2), na.rm =  TRUE))
b<-(xym-xm*ym)/(x2m-xm^2)
a<-ym-b*xm
v0<-b*x+a
v<-mrate*v0
write.table(round(v,5), file = "C:/Users/Рита/Desktop/rate_mm.csv", sep = ";", col.names=TRUE, row.names=FALSE)
getwd()
setwd("C:/VS-model/jakutsk_clim")
temp <- list.files(pattern="*.cli")
myfiles <- lapply(temp, read.delim, sep=" ", header=F)
df <- do.call(rbind.data.frame, myfiles)
write.table(df, file = "C:/Users/Рита/Desktop/all_years.csv", sep = ";", col.names=FALSE, row.names=FALSE)
getwd()
setwd("C:/VS-model/jakutsk_clim")
temp <- list.files(pattern="*.cli")
myfiles <- lapply(temp, read.delim, sep=" ", header=F)
df <- do.call(rbind.data.frame, myfiles)
write.table(df, file = "C:/Users/Рита/Desktop/all_years.csv", sep = ";", col.names=FALSE, row.names=FALSE)
getwd()
setwd("C:/VS-model/jakutsk_clim")
temp <- list.files(pattern="*.cli")
myfiles <- lapply(temp, read.delim, sep="  ", header=F)
df <- do.call(rbind.data.frame, myfiles)
write.table(df, file = "C:/Users/Рита/Desktop/all_years.csv", sep = ";", col.names=FALSE, row.names=FALSE)
getwd()
setwd("C:/VS-model/jakutsk_clim")
temp <- list.files(pattern="*.cli")
myfiles <- lapply(temp, read.delim, sep="\s", header=F)
df <- do.call(rbind.data.frame, myfiles)
write.table(df, file = "C:/Users/Рита/Desktop/all_years.csv", sep = ";", col.names=FALSE, row.names=FALSE)
getwd()
setwd("C:/VS-model/jakutsk_clim")
temp <- list.files(pattern="*.cli")
myfiles <- lapply(temp, read.delim, sep="   ", header=F)
df <- do.call(rbind.data.frame, myfiles)
write.table(df, file = "C:/Users/Рита/Desktop/all_years.csv", sep = ";", col.names=FALSE, row.names=FALSE)
getwd()
setwd("C:/VS-model/jakutsk_clim")
temp <- list.files(pattern="*.cli")
myfiles <- lapply(temp, read.delim, sep=" ", header=F)
df <- do.call(rbind.data.frame, myfiles)
write.table(df, file = "C:/Users/Рита/Desktop/all_years.csv", sep = ";", col.names=FALSE, row.names=FALSE)
getwd()
setwd("C:/VS-model/shira_clim")
temp <- list.files(pattern="*.cli")
myfiles <- lapply(temp, read.delim, sep=" ", header=F)
df <- do.call(rbind.data.frame, myfiles)
write.table(df, file = "C:/Users/Рита/Desktop/all_years.csv", sep = ";", col.names=FALSE, row.names=FALSE)
getwd()
setwd("C:/VS-model/sira_clim")
temp <- list.files(pattern="*.cli")
myfiles <- lapply(temp, read.delim, sep=" ", header=F)
df <- do.call(rbind.data.frame, myfiles)
write.table(df, file = "C:/Users/Рита/Desktop/all_years.csv", sep = ";", col.names=FALSE, row.names=FALSE)
