R语言计算生理发育时间、有效积温

# 代码注解：
# 首先加载程序包，若程序包加载失败表示R未安装该程序包。请运行"#"之后代码。
# 生理发育时间从播种那天开始计算
# RET表示相对热效应，PDT表示生理发育时间。
# tb 作物生长发育下限温度
# to 作物生长发育最适温度
# tm 作物生长发育上限温度
#RET = 0 t<tb;
# RET = (t-tb)/(to-tb) tb<=t<=to;
# RET = (tm-t)/(tm-to) to<t<=tm
# RET = 0 t>tm;
# RET=0.5*RET_av+0.25*RET_max+0.25*RET_min
# TT表示有效积温 ,年份='2020'
rm(list = ls())#该代码含义是清除历史记录相当于格式化
library(rio)
#install.packages('rio')
library(tidyverse)
# install.packages("tidyverse")
library(openxlsx)
ret <- import('生理发育.xlsx') %>%
rename(tmin='最低气温',
tmax='最高气温',
tav='平均气温') %>% filter(年=='2020')
#rename()函数修改列名
View(ret)
head(ret)

tb <- 12
to <- 30
tm <- 35
t_o <- 25

n <- ret$年
# d<- as.Date('2020-4-24')-as.Date('2020-1-1')+1
d<- as.Date('2020-4-24')-as.Date('2020-1-1')#4-24表示播种时间
d
ret%>% mutate(doy=1:length(n),
DAS=ifelse(doy>=d,(doy-d),0))->rett
# %>% filter(doy>=d)->ret
view(rett)
# 设定温度三基点
# tb 作物生长发育下限温度
# to 作物生长发育最适温度（设置了25℃、30℃）
# tm 作物生长发育上限温度
t <- ret$tav
t_1 <- ret$tmax
t_2 <- ret$tmin
ret %>% mutate(RET_av=(ifelse(t<12,0,
ifelse(12<=t&t<=30,(t-tb)/(to-tb),
ifelse(30<=t&t<=35,(tm-t)/(tm-to),
ifelse(t>35,0,0))))),
RET_max=(ifelse(t_1<12,0,
ifelse(12<=t_1&t_1<=30,(t_1-tb)/(to-tb),
ifelse(30<=t_1&t_1<=35,(tm-t_1)/(tm-to),
ifelse(t_1>35,0,0))))),
RET_min=(ifelse(t_2<12,0,
ifelse(12<=t_2&t_2<=30,(t_2-tb)/(to-tb),
ifelse(30<=t_2&t_2<=35,(tm-t_2)/(tm-to),
ifelse(t_2>35,0,0))))),
RET_30=0.5*RET_av+0.25*RET_max+0.25*RET_min,
PDT_30=cumsum(RET_30),
RET_a_v=(ifelse(t<12,0,
ifelse(12<=t&t<=25,(t-tb)/(t_o-tb),
ifelse(25<=t&t<=35,(tm-t)/(tm-t_o),
ifelse(t>35,0,0))))),
RET_m_ax=(ifelse(t_1<12,0,
ifelse(12<=t_1&t_1<=25,(t_1-tb)/(t_o-tb),
ifelse(25<=t_1&t_1<=35,(tm-t_1)/(tm-t_o),
ifelse(t_1>35,0,0))))),
RET_m_in=(ifelse(t_2<12,0,
ifelse(12<=t_2&t_2<=25,(t_2-tb)/(t_o-tb),
ifelse(25<=t_2&t_2<=35,(tm-t_2)/(tm-t_o),
ifelse(t_2>35,0,0))))),
RET_25=0.5*RET_a_v+0.25*RET_m_ax+0.25*RET_m_in,
PDT_25=cumsum(RET_25),
TT=cumsum(ifelse(t>12,(t-tb),0))) %>% rename('平均气温'=tav,
'最低气温'=tmin,
'最高气温'=tmax) %>% select(-RET_max,
-RET_av,
-RET_min,
-RET_25,-RET_30,-RET_a_v,
-RET_m_ax,-RET_m_in) -> ret_1

view(ret_1)#查看计算结果
export(ret_1,'生理发育时间.xlsx')#保存数据

R语言计算生理发育时间、有效积温相关推荐

r语言计算模型的rmse_直播丨R语言与作物模型高级应用实战技术应用
随着基于过程的作物生长模型(Process-based Crop Growth Simulation Model)的发展,R语言在作物生长模型和数据分析.挖掘和可视化中发挥着越来越重要的作用.想要成为 ...
R语言计算资本资产定价模型(CAPM)中的Beta值和可视化
原文链接:http://tecdat.cn/?p=22588 今天我们将计算投资组合收益的CAPM贝塔.这需要拟合一个线性模型,得到可视化,从资产收益的角度考虑我们的结果的意义. 简单的背景介绍,资本 ...
R语言计算曼哈顿距离（Manhattan Distance）实战:计算两个向量的曼哈顿距离、dist函数计算矩阵中两两元素的曼哈顿距离
R语言计算曼哈顿距离(Manhattan Distance)实战:计算两个向量的曼哈顿距离.dist函数计算矩阵中两两元素的曼哈顿距离目录 R语言计算曼哈顿距离(Manhattan Distance ...
R语言计算杰卡德相似系数（Jaccard Similarity）实战：自定义函数计算Jaccard相似度、对字符串向量计算Jaccard相似度、将Jaccard相似度转化为Jaccard距离
R语言计算杰卡德相似系数(Jaccard Similarity)实战:自定义函数计算Jaccard相似度.对字符串向量计算Jaccard相似度.将Jaccard相似度转化为Jaccard距离目录 R ...
R语言计算平均值的标准误差（standard error of the mean）：自定义函数计算平均值的标准误差、使用plotrix包的std.error函数计算平均值的标准误差
R语言计算平均值的标准误差(standard error of the mean):自定义函数计算平均值的标准误差.使用plotrix包的std.error函数计算平均值的标准误差目录
R语言计算回归模型每个样本（观察、observation、sample）的DFFITS度量实战：忽略单个观察（样本）时，回归模型所做的预测会发生多大的变化
R语言计算回归模型每个样本(观察.observation.sample)的DFFITS度量实战:忽略单个观察(样本)时,回归模型所做的预测会发生多大的变化目录
R语言计算回归模型每个样本（观察、observation、sample）的杠杆值（leverage）实战：如果一个样本的预测变量比其他样本的预测变量值更极端，那么被认为具有很高的杠杆作用
R语言计算回归模型每个样本(观察.observation.sample)的杠杆值(leverage)实战:如果一个样本的预测变量比其他样本的预测变量值更极端,那么被认为具有很高的杠杆作用目录
R语言计算每个分组的行数并将结果添加到dataframe中实战
R语言计算每个分组的行数并将结果添加到dataframe中实战目录 R语言计算每个分组的行数并将结果添加到dataframe中实战 #仿真数据
R语言计算F1评估指标实战：F1 score、使用R中caret包中的confusionMatrix（）函数为给定的logistic回归模型计算F1得分（和其他指标）
R语言计算F1评估指标实战:F1 score.使用R中caret包中的confusionMatrix()函数为给定的logistic回归模型计算F1得分(和其他指标) 目录

R语言计算生理发育时间、有效积温

R语言计算生理发育时间、有效积温相关推荐

最新文章

热门文章