安装.包( " 斯托克实验室 " )
# install.packages(“devtools”)
开发工具 :: 安装github( " HiDef-航空测量/stochLAB " )
停止 <- 斯托克实验室 :: stoch_crm(斯托克_克莱姆)( ... ) 停止 $ 碰撞 $ opt1(操作1) # 对于随机碰撞风险模型选项1的输出
停止 $ 碰撞 $ 操作2 # 对于随机碰撞风险模型选项2的输出
停止 $ 碰撞 $ 选项3 # 对于随机碰撞风险模型选项3的输出
图书馆( 斯托克实验室 ) # ------------------------------------------------------
# 预先设置一些必需的输入
b_dens(密度) <- 数据帧 ( 月 = abb月 , 意思是 = 运行( 12 , 0.8 , 1.5 ), 标准偏差 = 运行( 12 , 0.2 , 0.3 )) # Johnson等人(2014)针对一个物种的通用FHD自举
fhd_靴子 <- 通用_fhd_bootstraps [[ 1 ]] # 风速vs转速vs螺距
风_ rtn _ patch <- 数据帧 ( 风速 = seq_len( 30 ), rtn速度 = 10 / ( 30 : 1 ), bld_间距 = c(代表( 90 , 4 ),代表( 0 , 8 ), 5 : 22 )) # 风力可用性
风力发电机 <- 数据帧 ( 月 = abb月 , 多氯联苯(pctg) = 运行( 12 , 85 , 98 )) # 维护停机时间
载重吨 <- 数据帧 ( 月 = abb月 , 意思是 = 运行( 12 , 6 , 10 ), 标准偏差 = 代表( 2 , 12 )) # 季节规范
海_dt <- 数据帧 ( 季节_id = c(c)( " 一 " , " b条 " , " c(c) " ), 开始_月 = c(c)( " 简 " , " 五月 " , " 10月 " ), 月末 = c(c)( " 四月 " , " 九月 " , " 12月 " )) # ----------------------------------------------------------
# 运行随机CRM,处理转子半径、气隙和
# 叶片宽度作为固定参数(即非随机) stoch_crm(斯托克_克莱姆)( 模型选项(_O) = c(c)( 1 , 2 , 三 ), n比特 = 1000 , flt速度pars = 数据帧 ( 意思是 = 7.26 , 标准偏差 = 1.5 ), 车身_零件 = 数据帧 ( 意思是 = 0.39 , 标准偏差 = 0.005 ), 翼展杆 = 数据帧 ( 意思是 = 1.08 , 标准偏差 = 0.04 ), 避免bsc_pars = 数据帧 ( 意思是 = 0.99 , 标准偏差 = 0.001 ), 避免发送条形码 = 数据帧 ( 意思是 = 0.96 , 标准偏差 = 0.002 ), noct_act_pars = 数据帧 ( 意思是 = 0.033 , 标准偏差 = 0.005 ), prop_crh_pars = 数据帧 ( 意思是 = 0.06 , 标准偏差 = 0.009 ), 鸟巢角 = " tnorm(表格) " , 鸟鸣 = b_dens(密度) , 航班类型 = " 拍打 " , prop_upwind(上风) = 0.5 , 发电机_fhd_boots = fhd_靴子 , n个叶片 = 三 , rtr半径pars = 数据帧 ( 意思是 = 80 , 标准偏差 = 0 ), # sd=0,转子半径固定
空气间隙棒 = 数据帧 ( 意思是 = 36 , 标准偏差 = 0 ), # sd=0,气隙固定
bld宽度条 = 数据帧 ( 意思是 = 8 , 标准偏差 = 0 ), # sd=0,叶片宽度固定
rtn俯仰角 = " 风速Reltn " , windspd_pars(挡风玻璃面板) = 数据帧 ( 意思是 = 7.74 , 标准偏差 = 三 ), rtn俯仰风pd dt = 风_rtn_ptch , trb_wind_avbl = windavb公司 , trb_downtime_pars = 载重吨 , wf_n_标签 = 200 , wf_width(宽度) = 15 , wf_纬度 = 56.9 , 潮汐_偏移 = 2.5 , lrg_错误_错误 = 真的 , 冗长的 = 真的 , 种子 = 1234 , 超出格式(_F) = " 摘要 " , 输出采样条 = 真的 , 非周期 = " 季节 " , 季节性规格 = 海_dt , 日志文件(_F) = 粘贴0(getwd(), " scrm_示例.log " ) ) # >
# >——随机CRM——
# >
# > ℹ 检查输入 ✔ 检查输入[87ms]
# > ℹ 正在准备数据 ✔ 准备数据[152ms]
# > ℹ 采样参数 ✔ 采样参数[423ms]
# >计算碰撞|3/1000次迭代计算碰撞|117/1000次迭代计算冲突|272/1000次迭代 计算碰撞|1000/1000次迭代 ✔ 计算碰撞|1000/1000次迭代[1.6s]
# > ℹ 排序输出 ✔ 分拣输出[761ms]
# > ✔ 工作完成了!
# >$碰撞
# >$碰撞$opt1
# >#A台:3×10
# >季节id周期平均sd中位数pctl_2.5 pctl_25 pctl_75 pctl_9
# ><chr><chr><dbl><dbl><dbl><dbl
# >1 a一月_四月28.0 12.7 28.8 7.15 18.2 36.4 52.0 67.3
# >2 b五月_第56.4 25.3 58.1 15.4 37.2 73.5 105页。 142
# >3 c 10月_ 12月19.3日8.92 19.5 5.13 12.2 25.5 36.6 45.4
# >#…带缩写变量名¹pctl_97.5
# >
# >$碰撞$opt2
# >#A台:3×10
# >季节周期平均值sd中值pctl_2.5 pctl_25 pctl_75 pctl_9…¹pctl_99
# ><chr><chr><dbl><dbl><dbl><dbl
# >1 a 1月_Apr 0.728 1.28 0.393 0.0676 0.230 0.655 5.05 11.6
# >2 b五月_第1.47步2.64 0.804 0.135 0.469 1.30 10.8 21.8
# >3 c 10月12日0.499 0.888 0.274 0.0448 0.158 0.445 3.43 8.06
# >#…带缩写变量名¹pctl_97.5
# >
# >$碰撞$opt3
# >#A台:3×10
# >季节周期平均值sd中值pctl_2.5 pctl_25 pctl_75 pctl_9…¹pctl_99
# ><chr><chr><dbl><dbl><dbl><dbl
# >1 a一月_Apr 0.350 0.790 0.143 0.0337 0.0913 0.256 3.05 6.66
# >2 b五月_步骤0.710 1.61 0.294 0.0687 0.185 0.487 5.55 15.2
# >3 c 10月_ 12月0.240 0.546 0.0995 0.0223 0.0636 0.177 1.93 4.73
# >#…带缩写变量名¹pctl_97.5
# >
# >
# >$sampled_pars
# >$sampled_pars$空气间隙
# >#A台:1×5
# >平均sd中位数pctl2.5pctl97.5
# ><dbl><dbl><dbl><dbl><dbl>
# > 1 36 0 36 36 36
# >
# >$sampled_pars$bld_width
# >#A台:1×5
# >平均sd中位数pctl2.5pctl97.5
# ><dbl><dbl>
# > 1 8 0 8 8 8
# >
# >$sampled_pars$body_lt
# >#A兼容:1×5
# >平均sd中位数pctl2.5pctl97.5
# ><dbl><dbl>
# > 1 0.390 0.00499 0.390 0.380 0.400
# >
# >$sampled_pars$flt_speed(采样部分$flt_speed)
# >#A台:1×5
# >平均sd中位数pctl2.5pctl97.5
# ><dbl><dbl>
# > 1 7.28 1.47 7.28 4.30 10.0
# >
# >$sampled_pars$noct_actv
# >#A台:1×5
# >平均sd中位数pctl2.5pctl97.5
# ><dbl><dbl>
# > 1 0.0333 0.00498 0.0333 0.0241 0.0436
# >
# >$sampled_pars$rtr_radius
# >#A台:1×5
# >平均sd中位数pctl2.5pctl97.5
# ><dbl><dbl>
# > 1 80 0 80 80 80
# >
# >$sampled_pars$翼展
# >#A台:1×5
# >平均sd中位数pctl2.5pctl97.5
# ><dbl><dbl>
# > 1 1.08 0.0398 1.08 1.00 1.16
# >
# >$采样条$hub_height
# >#A台:1×5
# >平均sd中位数pctl2.5pctl97.5
# ><dbl><dbl>
# > 1 116 0 116 116 116
# >
# >$sampled_pars$dens_mth
# >#A tibble:12×6
# >周期平均sd中位数pctl2.5 pctl97.5
# ><chr><dbl><dbl><dbl><dbl
# >1月1日0.849 0.286 0.847 0.315 1.42
# >2月2日0.833 0.224 0.826 0.415 1.28
# >2012年3月3日0.225 1.12 0.665 1.53
# >2011年4月4日0.242 1.11 0.647 1.59
# >5月5日1.22 0.270 1.21 0.700 1.75
# >6月6日0.835 0.252 0.838 0.299 1.35
# >2007年7月7日0.247 1.07 0.600 1.55
# >8月8日1.04 0.279 1.05 0.490 1.61
# >2013年9月9日0.277 1.13 0.601 1.66
# >2014年10月10日0.222 1.04 0.616 1.48
# >11月11日0.965 0.267 0.959 0.460 1.50
# >2011年12月12日0.199 1.21 0.813 1.59
# >
# >$sampled_pars$prop_oper_mth
# >#A台:12×6
# >周期平均sd中位数pctl2.5 pctl97.5
# ><chr><dbl><dbl><dbl><dbl
# >1月1日0.793 0.0207 0.792 0.752 0.829
# >2月2日0.769 0.0207 0.769 0.729 0.808
# >3月3日0.850 0.0204 0.850 0.807 0.890
# >4月4日0.807 0.0196 0.807 0.769 0.845
# >5月5日0.847 0.0195 0.847 0.809 0.884
# >6月6日0.866 0.0198 0.866 0.827 0.905
# >7月7日0.869 0.0199 0.869 0.829 0.910
# >8月8日0.806 0.0194 0.806 0.768 0.844
# >9月9日0.879 0.0205 0.879 0.838 0.917
# >10月10日0.820 0.0200 0.820 0.779 0.858
# >11月11日0.858 0.0201 0.858 0.819 0.898
# >12月12日0.910 0.0201 0.911 0.872 0.947
# >
# >$sampled_pars$停机
# >#A台:12×6
# >周期平均sd中位数pctl2.5 pctl97.5
# ><chr><dbl><dbl><dbl><dbl
# >2009年1月1日2.07 6.12 2.39 10.1
# >2008年2月2日2.07 8.99 5.12 12.9
# >3月3日8.54 2.04 8.54 4.52 12.8
# >4月4日8.57 1.96 8.61 4.79 12.4
# >5月5日7.83 1.95 7.81 4.12 11.7
# >6月6日7.64 1.98 7.61 3.78 11.6
# >7月7日9.43 1.99 9.42 5.37 13.5
# >8月8日6.35 1.94 6.31 2.55 10.1
# >9月9日7.46 2.05 7.44 3.62 11.6
# >2008年10月10日2.00 6.10 2.28 10.1
# >11月11日7.56 2.01 7.54 3.58 11.4
# >12月12日7.03 2.01 6.94 3.32 10.8
# >
# >$sampled_pars$wind_speed
# >#A台:1×5
# >平均sd中位数pctl2.5pctl97.5
# ><dbl><dbl><dbl><dbl><dbl>
# > 1 7.73 3.00 7.73 2.15 13.7
# >
# >$sampled_pars$rtn_speed
# >#A台:1×5
# >平均sd中位数pctl2.5pctl97.5
# ><dbl><dbl>
# > 1 0.428 0.0574 0.417 0.345 0.556
# >
# >$取样棒$bld_pitch
# >#A台:1×5
# >平均sd中位数pctl2.5pctl97.5
# ><dbl><dbl>
# > 1 0.329 0.635 0 0 1.57
# >
# >$sampled_pars$avoid_bsc
# >#A台:1×5
# >平均sd中位数pctl2.5pctl97.5
# ><dbl><dbl>
# > 1 0.990 0.000985 0.990 0.988 0.992
# >
# >$sampled_pars$避免文本
# >#A台:1×5
# >平均sd中位数pctl2.5pctl97.5
# ><dbl><dbl>
# > 1 0.960 0.00203 0.960 0.956 0.964
# >
# >$sampled_pars$prop_crh
# >#A台:1×5
# >平均sd中位数pctl2.5pctl97.5
# ><dbl><dbl>
# > 1 0.0605 0.00904 0.0601 0.0441 0.0800
# >
# >$sampled_pars$gen_fhd
# >#A桥:500×6
# >身高平均sd中位数pctl2.5 pctl97.5
# ><dbl><dbl>
# > 1 0 0.163 0.0188 0.166 0.109 0.187
# > 2 1 0.136 0.0129 0.138 0.0967 0.152
# > 3 2 0.114 0.00855 0.115 0.0863 0.124
# > 4 3 0.0950 0.00530 0.0963 0.0769 0.100
# > 5 4 0.0794 0.00296 0.0803 0.0686 0.0816
# > 6 5 0.0664 0.00146 0.0669 0.0606 0.0670
# > 7 6 0.0556 0.00116 0.0558 0.0530 0.0567
# > 8 7 0.0465 0.00166 0.0466 0.0439 0.0490
# > 9 8 0.0390 0.00216 0.0389 0.0357 0.0432
# > 10 9 0.0327 0.00250 0.0324 0.0290 0.0386
# >#…还有490行
-
建议如何在上一版本中使用输入参数数据集 可以重新调整实现以适应 stoch_crm() 的接口。 建议的代码也与多个场景的上下文相关 应用程序,因为这些数据集的宽表结构是 可能是用户编译输入参数的首选格式 在不同的场景下。 -
提出要运行的功能性编程框架 stoch_crm() 对于 多品种和风机/涡轮机功能。
图书馆( 斯托克实验室 ) # --------------------------------------------------------- #
# ----重塑为列表-列数据帧----
# --------------------------------------------------------- #
#
# ---鸟类特征
小鸟球 <- bird_pars_wide_示例 % > % dplyr公司 :: 重新安置( 航班 , .之后 = dplyr公司 :: last_col())% > % 第三年 :: 枢轴加长器( 回避基本 : 属性CRH_ObsSD ) % > % dplyr公司 :: 变异( 标准 = dplyr公司 :: if_else(grepl( " SD|SD|SD " , 名称 ), " 标准偏差 " , " 意思是 " ), 特征 = gsub公司( " 标准差|标准差|标准差 " , " " , 名称 )) % > % dplyr公司 :: 选择( - 名称 ) % > % 第三年 :: 枢轴_加宽( 名称来自 = 标准 , 值_来自 = 价值 ) % > % 第三年 :: 巢( 部分 = c(c)( 意思是 , 标准偏差 )) % > % 第三年 :: 枢轴_加宽( 名称来自 = 特征 , 值_来自 = 部分 ) % > % 易怒的 :: add_column(添加列)( prop_upwind(上风) = 0.5 ) # ---鸟类密度:提供为截断正常下限的平均值和sd参数
# 以0为界
dens_pars(密集阵列) <- 密度_标称_宽度_示例 % > % 易怒的 :: add_column(添加列)( 密度(_O) = 代表( " tnorm(表格) " ,当前( . )), .之后 = 1 ) % > % 第三年 :: 枢轴加长器( 简 : DecSD公司 ) % > % dplyr公司 :: 突变( 标准 = dplyr公司 :: if_else(grepl( " SD|SD|SD " , 名称 ), " 标准偏差 " , " 意思是 " ), 月 = gsub公司( " 标准差|标准差|标准差 " , " " , 名称 )) % > % dplyr公司 :: 选择( - 名称 ) % > % 第三年 :: 枢轴_加宽( 名称来自 = 标准 , 值_来自 = 价值 ) % > % 第三年 :: 巢( mth密度 = c(c)( 月 , 意思是 , 标准偏差 )) # ---Johnson等人(2014年)关于分析物种的FHD数据
发电机_fhd_boots <- generic_fhd_脚踏板 [ 鸟帕斯 $ 物种 ] # ---季节定义(合成)
季节_dt <- 列表 ( 弧形_弧形 = 数据帧 ( 季节id = c(c)( " 育种 " , " 喂养 " , " 迁移 " ), 开始_月 = c(c)( " 五月 " , " 九月 " , " 简 " ), 月末 = c(c)( " 八月 " , " 12月 " , " 四月 " )), 黑头鸥 = 数据帧 ( 季节id = c(c)( " 育种 " , " 喂养 " , " 迁移 " ), 开始_月 = c(c)( " 简 " , " 五月 " , " 10月 " ), 月末 = c(c)( " 四月 " , " 九月 " , " 12月 " )), 黑腿小猫 = 数据帧 ( 季节id = c(c)( " 育种 " , " 喂养 " , " 迁移 " ), 开始_月 = c(c)( " 12月 " , " 3月 " , " 九月 " ), 月末 = c(c)( " 二月 " , " 八月 " , " 十一月 " ))) # ---汽轮机参数
# #首先处理操作参数
trb_opr_pars <- turb_pars_wide_example(涡轮机_零件_宽度_示例) % > % dplyr公司 :: 选择( 涡轮机模型 , 1月 : DecOpSD公司 ) % > % 第三年 :: 枢轴加长器( JanOp公司 : DecOpSD公司 ) % > % dplyr公司 :: 突变( 月 = 子(substr)( 名称 , 1 , 三 ), 标准 = dplyr公司 :: 案例_时间( 格雷普( " SD|SD|SD " , 名称 ) ~ " 标准偏差 " , 格雷普( " 平均值|平均值|平均值 " , 名称 ) ~ " 意思是 " , 真的 ~ " 多氯联苯(pctg) " )) % > % dplyr公司 :: 选择( - 名称 ) % > % 第三年 :: 枢轴_加宽( 名称来自 = 标准 , 值_来自 = 价值 ) % > % 第三年 :: 巢( 风_矢量 = c(c)( 月 , 多氯联苯(pctg) ), trb_dwntm公司 = c(c)( 月 , 意思是 , 标准偏差 )) # #解决涡轮机特性并随后合并操作参数
trb_pars公司 <- turb_pars_wide_example(涡轮机_零件_宽度_示例) % > % dplyr公司 :: 选择( 涡轮机模型 : 风速SD ) % > % dplyr公司 :: 重新安置( 转子速度和俯仰_Sim选项 , .之后 = 1 ) % > % 第三年 :: 枢轴加长器( 转子半径 : 风速SD ) % > % dplyr公司 :: 突变( 标准 = dplyr公司 :: if_else(grepl( " SD|SD|SD " , 名称 ), " 标准偏差 " , " 意思是 " ), 特征 = gsub公司( " (SD|SD|SD)|(平均值|Mean|Mean) " , " " , 名称 ) ) % > % dplyr公司 :: 选择( - 名称 ) % > % 第三年 :: 枢轴(_wider)( 名称来自 = 标准 , 值_来自 = 价值 ) % > % 第三年 :: 巢( 部分 = c(c)( 意思是 , 标准偏差 )) % > % 第三年 :: 枢轴(_wider)( 名称来自 = 特征 , 值_来自 = 部分 ) % > % dplyr公司 :: left_join(左_连接)( . , trb_opr_pars ) # >通过=“TurbineModel”连接
# ---风速、转速和叶片节距关系
wndspd_rtn_ptch示例
# >风速rtn速度bld变桨
# > 1 0 0.0 90
# > 2 1 0.0 90
# > 3 2 0.0 90
# > 4 3 6.8 0
# > 5 4 6.8 0
# > 6 5 6.8 0
# > 7 6 6.8 0
# > 8 7 6.8 0
# > 9 8 8.1 0
# > 10 9 9.1 0
# > 11 10 9.3 0
# > 12 11 9.4 4
# > 13 12 9.5 7
# > 14 13 9.7 9
# > 15 14 9.7 11
# > 16 15 9.9 13
# > 17 16 10.2 15
# > 18 17 10.2 16
# > 19 18 10.2 18
# > 20 19 10.2 19
# > 21 20 10.2 20
# > 22 21 10.2 22
# > 23 22 10.2 23
# > 24 23 10.2 24
# > 25 24 10.2 25
# > 26 25 10.2 26
# > 27 26 10.2 27
# > 28 27 10.2 28
# > 29 28 10.2 29
# > 30 29 10.2 30
# ---风电场参数
wf-par(_P) <- 数据帧 ( 工作流程id = c(c)( " wf_1型 " , " 工作流2 " ), n_涡轮机 = c(c)( 200 , 400 ), wf_宽度 = c(c)( 4 , 10 ), wf_标签 = c(c)( 55.8 , 55 ), td_关闭 = c(c)( 2.5 , 2 ), 大数组corr = c(c)( 错误的 , 真的 ) ) # -------------------------------------------------------------- #
# ----针对多个场景运行stoch_crm()----
# -------------------------------------------------------------- #
# ---设置场景组合
场景_规格 <- 第三年 :: 展开网格(_grid)( 特殊目的公司 = 鸟帕斯 $ 物种 , turb_id(土耳其_id) = trb_pars公司 $ 涡轮机模型 , 工作流程id = wf-par(_P) $ 工作流程id ) % > % 易怒的 :: add_column(添加列)( 场景_id = 粘贴0( " 场景_ " , 1 : nrow公司( . )), .之前 = 1 ) # ---为即将到来的迭代映射步骤设置进度条
铅 <- 进步 :: 进度栏 $ 新( 格式 = " 运行场景::what[:bar]:百分比eta::eta " , 宽度 = 100 , 全部的 = nrow公司( 场景_规格 )) # ---通过purrr::pmap将stoch_crm()映射到每个场景规范
输出 <- 场景_规格 % > % 呜呜声 :: pmap程序( 功能 ( 场景_id , 特殊目的公司 , turb_id(土耳其_id) , 工作流程id , ... ){ 铅 $ 滴答声( 代币 = 列表 ( 什么 = 场景_id )) # 当前物种的参数
c_spec(规范) <- 小鸟球 % > % dplyr公司 :: 过滤器( 物种 == {{ 特殊目的公司 }}) # 当前物种密度
c_dens(十进制) <- 密集路径 % > % dplyr公司 :: 过滤器( 物种 == {{ 特殊目的公司 }}) # 当前涡轮机场景的参数
c_turb(_turb) <- trb_pars公司 % > % dplyr公司 :: 过滤器( 涡轮机模型 == {{ turb_id(土耳其_id) }}) # 当前风电场场景的参数
c_wf公司 <- wf-par(_P) % > % dplyr公司 :: 过滤器( 工作流程id == {{ 工作流程id }}) # 列表列中的输入需要未列出,可以通过“unlist()”或
# 索引`[1]`
# 关闭“verbose”,否则控制台将
# 充斥着日志消息 stoch_crm(斯托克_克莱姆)( 模型选项(_O) = c(c)( 1 , 2 , 三 ), n比特 = 1000 , flt速度pars = 规格(_spec) $ 飞行_速度 [[ 1 ]], 车身_零件 = c_spec(规范) $ 车身_长度 [[ 1 ]], 机翼展开路径 = c_spec(规范) $ 翼展 [[ 1 ]], 避免使用bsc_pars = c_spec(规范) $ 回避基本 [[ 1 ]], 避免发送条形码 = c_spec(规范) $ 回避延长 [[ 1 ]], 无动作路径 = c_spec(规范) $ 夜间活动 [[ 1 ]], prop_crh_pars = 规格(_spec) $ 属性CRH_Obs [[ 1 ]], 鸟巢角 = c_dens(十进制) $ 密度(_O) , 鸟巢 = c_dens(十进制) $ mth密度 [[ 1 ]], 航班类型 = c_spec(规范) $ 航班 , prop_upwind(上风) = c_spec(规范) $ prop_upwind(上风) , 发电机_fhd_boots = 发电机_fhd_boots [[ 特殊目的公司 ]], n个叶片 = c_turb(_turb) $ 叶片 , rtr半径pars = c_turb(_turb) $ 转子半径 [[ 1 ]], 空气间隙棒 = c_turb(_turb) $ 轮毂高度添加 [[ 1 ]], bld宽度条 = c_turb(_turb) $ 叶片宽度 [[ 1 ]], rtn俯仰角 = c_turb(_turb) $ 转子速度和俯仰_Sim选项 , bld变桨杆 = c_turb(_turb) $ 螺距 [[ 1 ]], rtn速度pars = c_turb(_turb) $ 旋转速度 [[ 1 ]], windspd_pars(挡风玻璃面板) = c_turb(_turb) $ windSpeed(风速) [[ 1 ]], rtn俯仰风pd dt = wndspd_rtn_ptch示例 , trb_wind_avbl = c_turb(_turb) $ 风_矢量 [[ 1 ]], trb_downtime_pars = c_turb(_turb) $ trb_dwntm公司 [[ 1 ]], wf_n_标签 = c_wf公司 $ n_涡轮机 , wf_width(宽度) = c_wf公司 $ wf_width(宽度) , wf_纬度 = c_wf公司 $ wf_标签 , 潮汐_偏移 = c_wf公司 $ td_关闭 , lrg_错误_错误 = c_wf公司 $ 大数组corr , 冗长的 = 错误的 , 种子 = 1234 , 超出格式(_F) = " 摘要 " , 输出采样条 = 错误的 , 非周期 = " 季节 " , 季节性规格 = 季节_dt [[ 特殊目的公司 ]], 日志文件(_F) = 无效的 ) }) # ---关闭进度条
铅 $ 终止() # ---识别输出列表的元素 姓名( 输出 ) <- 场景_规格 $ 场景_id
输出
# >$场景1
# >$场景_1$opt1
# >#A台:3×10
# >季节id周期平均sd中位数pctl_2.5 pctl_25 pctl_75 pctl_9
# ><chr><chr><dbl><dbl><dbl><dbl
# >1繁殖5月_Aug 984。 75.9 982. 847. 933. 1031. 1130. 1241
# >2号饲料,552年9月至12月。 41.6 551. 478. 523. 579. 634. 683
# >3于626年1月迁移。 47.7 623. 541. 591. 658. 722. 780
# >#…带缩写变量名¹pctl_97.5
# >
# >$scenario_1$opt2
# >#A台:3×10
# >季节周期平均值sd中值pctl_2.5 pctl_25 pctl_75 pctl_9…¹pctl_99
# ><chr><chr><dbl><dbl><dbl><dbl
# >1育种May_Aug 60.7 67.3 41.9 2.41 25.3 65.7 293。 389
# >2喂食Sep_Dec 34.0 37.7 23.5 1.35 14.3 37.0 162。 217
# >3迁移1月至6月38.6 42.8 26.8 1.55 16.2 42.4 186。 249
# >#…带缩写变量名¹pctl_97.5
# >
# >$scenario_1$opt3
# >#A台:3×10
# >季节周期平均值sd中值pctl_2.5 pctl_25 pctl_75 pctl_9…¹pctl_99
# ><chr><chr><dbl><dbl><dbl><dbl
# >1育种5月_Aug 17.9 24.8 10.5 0.394 5.92 17.9 106。 151
# >2进料Sep_Dec 10.1 13.9 5.86 0.221 3.35 10.0 59.0 84.5
# >3迁移Jan_Apr 11.4 15.8 6.77 0.251 3.79 11.5 66.9 93.8
# >#…带有缩写的变量名“”“”“”“”“”“pctl_97.5”
# >
# >
# >$场景2
# >$场景2$opt1
# >#A台:3×10
# >季节id周期平均sd中位数pctl_2.5 pctl_25 pctl_75 pctl_9
# ><chr><chr><dbl><dbl><dbl><dbl
# >1836年5月至8月繁殖。 137. 1831. 1586. 1743. 1921. 2098. 2301
# >2进料Sep_Dec 1045。 76.4 1044. 909. 991. 1096. 1195. 1288
# >1181年1月3日迁移。 87.6 1177. 1025. 1118. 1241. 1356. 1463
# >#…带缩写变量名¹pctl_97.5
# >
# >$场景2$opt2
# >#A台:3×10
# >季节周期平均值sd中值pctl_2.5 pctl_25 pctl_75 pctl_9…¹pctl_99
# ><chr><chr><dbl><dbl><dbl><dbl
# >1繁殖5月_Aug 122。 132. 85.2 5.20 52.1 133. 574. 759.
# >2喂食9月_Dec 69.2 75.0 48.4 2.95 29.9 75.9 323。 430.
# >3迁移Jan_Apr 78.3 84.8 55.1 3.37 33.7 86.5 369。 491
# >#…带缩写变量名¹pctl_97.5
# >
# >$场景2$opt3
# >#A台:3×10
# >季节周期平均值sd中值pctl_2.5 pctl_25 pctl_75 pctl_9…¹pctl_99
# ><chr><chr><dbl><dbl><dbl><dbl
# >1育种五月—八月35.9 48.7 21.4 0.848 12.2 36.1 208。 294
# >2喂食9月_Dec 20.4 27.7 12.1 0.482 7.00 20.6 118。 167
# >3迁移日期:1月23日3月31日9月13日0.547 7.91 23.5 133。 185
# >#…带缩写变量名¹pctl_97.5
# >
# >
# >$场景_3
# >$scenario_3$opt1
# >#A台:3×10
# >季节周期平均值sd中值pctl_2.5 pctl_25 pctl_75 pctl_9…¹pctl_99
# ><chr><chr><dbl><dbl><dbl><dbl
# >1育种5月_Aug 2131。 143. 2126. 1875. 2034. 2223. 2415. 2690
# >2喂食9月~1149年12月。 76.8 1143. 1011. 1095. 1197. 1302. 1446
# >1300年1月3日迁移。 89.6 1292. 1131. 1239. 1357. 1491. 1638
# >#…带缩写变量名¹pctl_97.5
# >
# >$scenario_3$opt2
# >#A台:3×10
# >季节周期平均值sd中值pctl_2.5 pctl_25 pctl_75 pctl_9…¹pctl_99
# ><chr><chr><dbl><dbl><dbl><dbl
# >1育种5月_Aug 53.7 77.6 30.1 0.894 17.0 53.9 347。 477.
# >2喂食9月_ 12月28日41.8 16.1 0.496 9.16 29.0 185。 257
# >3迁移Jan_Apr 32.7 47.1 18.3 0.546 10.3 32.5 205。 291
# >#…带缩写变量名¹pctl_97.5
# >
# >$scenario_3$opt3场景
# >#A台:3×10
# >季节周期平均值sd中值pctl_2.5 pctl_25 pctl_75 pctl_9…¹pctl_99
# ><chr><chr><dbl><dbl><dbl><dbl
# >1育种五月—八月14.4 25.7 6.47 0.124 3.37 12.7 109。 162
# >2进料Sep_Dec 7.74 13.8 3.48 0.0685 1.81 6.84 59.8 86.3
# >3迁移Jan_Apr 8.74 15.6 3.93 0.0756 2.05 7.65 66.9 96.0
# >#…带缩写变量名¹pctl_97.5
# >
# >
# >$场景_4
# >$scenario_4$opt1
# >#A台:3×10
# >季节周期平均值sd中值pctl_2.5 pctl_25 pctl_75 pctl_9…¹pctl_99
# ><chr><chr><dbl><dbl><dbl><dbl
# >1号繁殖5月_Aug 3720。 240. 3713. 3290. 3560. 3872. 4188. 4654
# >2喂养日期:2037年9月至12月。 131. 2025. 1801. 1944. 2122. 2297. 2540
# >3于2297年1月迁移。 152. 2284. 2012. 2193. 2393. 2626. 2854
# >#…带缩写变量名¹pctl_97.5
# >
# >$scenario_4$opt2
# >#A台:3×10
# >季节周期平均值sd中值pctl_2.5 pctl_25 pctl_75 pctl_9…¹pctl_99
# ><chr><chr><dbl><dbl><dbl><dbl
# >1号繁殖5月_Aug 100。 142. 57.0 1.79 32.6 101. 633. 867.
# >2喂食9月_Dec 54.8 77.7 31.0 1.01 17.8 55.5 342。 474
# >3迁移1月至5月61.7 87.3 35.1 1.11 20.1 62.2 382。 531
# >#…带缩写变量名¹pctl_97.5
# >
# >$scenario_4$opt3场景
# >#A兼容:3×10
# >季节周期平均值sd中值pctl_2.5 pctl_25 pctl_75 pctl_9…¹pctl_99
# ><chr><chr><dbl><dbl><dbl><dbl
# >1育种May_Aug 26.8 47.1 12.3 0.249 6.47 23.9 199。 295
# >2进料Sep_Dec 14.7 25.7 6.69 0.140 3.53 13.1 111。 160
# >3迁移1月5日16.5 28.9 7.57 0.153 3.98 14.7 124。 177
# >#…带缩写变量名¹pctl_97.5
# >
# >
# >$场景5
# >$场景5$opt1
# >#A台:3×10
# >季节周期平均值sd中值pctl_2.5 pctl_25 pctl_75 pctl_9…¹pctl_99
# ><chr><chr><dbl><dbl><dbl><dbl
# >1176年5月至11月繁殖。 341. 1053. 825. 985. 1219. 2220. 2634
# >2日喂食,9月至12月615日。 176. 551. 430. 516. 640. 1136. 1338
# >3迁移Jan_Apr 703。 202. 630. 489. 591. 731. 1323. 1551
# >#…带缩写变量名¹pctl_97.5
# >
# >$场景5$op2
# >#A台:3×10
# >季节id周期平均sd中位数pctl_2.5 pctl_25 pctl_75 pctl_9
# ><chr><chr><dbl><dbl><dbl><dbl
# >1育种5月_Aug 0.837 2.75 0.131 0.000231 0.0341 0.390 8.68 31.9
# >2进料Sep_Dec 0.438 1.45 0.0687 0.000118 0.0180 0.202 4.66 17.0
# >3迁移Jan_Apr 0.500 1.63 0.0782 0.000136 0.0209 0.230 5.20 18.4
# >#…带缩写变量名¹pctl_97.5
# >
# >$场景5$opt3
# >#A台:3×10
# >季节id周期平均sd中位数pctl_2.5 pctl_25 pctl_75 pctl_
# ><chr><chr><dbl><dbl><dbl><dbl
# >1育种5月_Aug 0.221 0.795 0.0245 0.0000261 0.00602 0.0786 2.41 8.88
# >2进料Sep_Dec 0.116 0.418 0.0129 0.0000129 0.00319 0.0408 1.29 4.73
# >3迁移Jan_Apr 0.132 0.472 0.0147 0.0000151 0.00365 0.0462 1.44 5.12
# >#…带缩写变量名¹pctl_97.5
# >
# >
# >$场景_6
# >$场景_6$opt1
# >#A台:3×10
# >季节周期平均值sd中值pctl_2.5 pctl_25 pctl_75 pctl_9…¹pctl_99
# ><chr><chr><dbl><dbl><dbl><dbl
# >1育种5月_Aug 2164。 568. 1961. 1557. 1840. 2263. 3874. 4500
# >1148年9月至12月2日喂食。 298. 1043. 822. 977. 1210. 2021. 2323
# >3于1310年1月迁移。 341. 1188. 934. 1117. 1370. 2325. 2688
# >#…带缩写变量名¹pctl_97.5
# >
# >$场景_6$op2
# >#A台:3×10
# >季节周期平均值sd中值pctl_2.5 pctl_25 pctl_75 pctl_9…¹pctl_99
# ><chr><chr><dbl><dbl><dbl><dbl><dbl><dbl><dbl><dbl><dbl>
# >1育种5月_Aug 1.64 5.28 0.266 0.000499 0.0711 0.779 17.1 59.5
# >2进料Sep_Dec 0.872 2.82 0.142 0.000254 0.0379 0.415 9.33 32.2
# >3迁移Jan_Apr 0.991 3.18 0.162 0.000295 0.0437 0.472 10.4 34.7
# >#…带缩写变量名¹pctl_97.5
# >
# >$场景_6$opt3
# >#A台:3×10
# >季节周期平均值sd中值pctl_2.5 pctl_25 pctl_75 pctl_…¹pctl_99
# ><chr><chr><dbl><dbl><dbl><dbl><dbl><dbl><dbl><dbl><dbl>
# >1育种5月_Aug 0.433 1.53 0.0503 0.0000579 0.0126 0.157 4.77 16.6
# >2进料Sep_Dec 0.230 0.817 0.0271 0.0000291 0.00672 0.0832 2.59 8.97
# >3迁移Jan_Apr 0.261 0.921 0.0304 0.0000337 0.00768 0.0947 2.88 9.68
# >#…带有缩写的变量名“”“”“”“”“”“pctl_97.5”
# >
# >
# >$场景_7
# >$场景_7$opt1
# >#A台:3×10
# >季节周期平均值sd中值pctl_2.5 pctl_25 pctl_75 pctl_9…¹pctl_99
# ><chr><chr><dbl><dbl><dbl><dbl
# >429年1月1日繁殖。 93.1 423. 270. 360. 490. 626. 767
# >5月2日喂食,第697步。 151. 685. 442. 588. 791. 1006. 1252
# >3迁移10月至次年12月266日。 59.2 262. 164. 223. 303. 390. 498
# >#…带缩写变量名¹pctl_97.5
# >
# >$场景_7$op2
# >#A台:3×10
# >季节周期平均值sd中值pctl_2.5 pctl_25 pctl_75 pctl_9…¹pctl_99
# ><chr><chr><dbl><dbl><dbl><dbl
# >1繁殖Jan_Apr 257。 150. 223. 62.9 153. 320. 603. 1071
# >5月2日喂食_第417步。 244. 366. 100. 251. 517. 985. 1797
# >3迁移10月至次年12月159日。 93.0 137. 37.8 94.0 199. 377. 649.
# >#…带缩写变量名¹pctl_97.5
# >
# >$场景_7$opt3
# >#A台:3×10
# >季节周期平均值sd中值pctl_2.5 pctl_25 pctl_75 pctl_9…¹pctl_99
# ><chr><chr><dbl><dbl><dbl><dbl
# >1育种Jan_Apr 28.2 25.9 20.9 4.02 12.5 34.1 86.4 239。
# >2喂食5月_步骤45.9 42.4 34.3 6.51 20.3 55.1 145。 409
# >3迁移10月_Dec 17.5 16.0 12.8 2.48 7.78 21.8 54.6 144。
# >#…带缩写变量名¹pctl_97.5
# >
# >
# >$场景_8
# >$场景_8$opt1
# >#A台:3×10
# >季节周期平均值sd中值pctl_2.5 pctl_25 pctl_75 pctl_9…¹pctl_99
# ><chr><chr><dbl><dbl><dbl><dbl
# >855年1月1日繁殖。 186. 844. 539. 719. 977. 1248. 1531
# >2号喂料5月_第1376步。 299. 1352. 872. 1160. 1560. 1985. 2472.
# >3迁移10月_Dec 534。 119. 526. 329. 448. 608. 782. 999
# >#…带缩写变量名¹pctl_97.5
# >
# >$场景_8$op2
# >#A兼容:3×10
# >季节周期平均值sd中值pctl_2.5 pctl_25 pctl_75 pctl_9…¹pctl_99
# ><chr><chr><dbl><dbl><dbl><dbl
# >1号繁殖,1月至536年。 307. 468. 135. 324. 668. 1246. 2183
# >5月2日喂食第862步。 494. 760. 213. 524. 1067. 2013. 3624
# >3迁移10月_Dec 334。 192. 290. 81.5 200. 418. 783. 1330
# >#…带缩写变量名¹pctl_97.5
# >
# >$场景_8$opt3
# >#A台:3×10
# >季节周期平均值sd中值pctl_2.5 pctl_25 pctl_75 pctl_9…¹pctl_99
# ><chr><chr><dbl><dbl><dbl><dbl
# >1育种Jan_Apr 58.6 52.8 43.8 8.62 26.4 71.2 178。 485
# >2喂食5月_步骤94.4 85.7 71.2 14.0 42.5 114。 295. 820.
# >3迁移10月_ 12月36.6 32.8 27.1 5.36 16.5 45.6 113。 293
# >#…带缩写变量名¹pctl_97.5
# >
# >
# >$场景_9
# >$scenario_9$opt1
# >#A台:3×10
# >季节周期平均值sd中值pctl_2.5 pctl_25 pctl_75 pctl_9…¹pctl_99
# ><chr><chr><dbl><dbl><dbl><dbl
# >1个繁殖期,1月至571年1月。 123. 563. 360. 484. 643. 818. 1015.
# >2喂食5月_第965步。 210. 953. 615. 815. 1089. 1428. 1738
# >3迁移10月_Dec 354。 78.1 349. 223. 299. 401. 511. 646
# >#…带缩写变量名¹pctl_97.5
# >
# >$场景_9$op2
# >#A台:3×10
# >季节周期平均值sd中值pctl_2.5 pctl_25 pctl_75 pctl_9…¹pctl_99
# ><chr><chr><dbl><dbl><dbl><dbl
# >194年1月1日繁殖。 140. 157. 32.9 100. 245. 539. 1050
# >5月2日喂食_第329步。 238. 265. 55.1 169. 414. 905. 1788
# >3迁移10月至12月121日。 87.9 97.6 20.1 62.0 153. 332. 682.
# >#…带缩写变量名¹pctl_97.5
# >
# >$场景_9$opt3
# >#A兼容:3×10
# >季节周期平均值sd中值pctl_2.5 pctl_25 pctl_75 pctl_9…¹pctl_99
# ><chr><chr><dbl><dbl><dbl><dbl
# >1育种Jan_Apr 23.0 25.6 15.7 2.36 8.83 27.4 79.6 265。
# >2喂食5月_步骤38.9 43.5 26.3 3.96 15.0 47.4 137。 472
# >3迁移10月_ 12月14日16.2 9.75 1.46 5.51 17.2 48.9 183。
# >#…带缩写变量名¹pctl_97.5
# >
# >
# >$场景_10
# >$场景_10$opt1
# >#A台:3×10
# >季节周期平均值sd中值pctl_2.5 pctl_25 pctl_75 pctl_9…¹pctl_99
# ><chr><chr><dbl><dbl><dbl><dbl><dbl><dbl><dbl><dbl><dbl>
# >1136年1月1日繁殖。 245. 1120. 717. 963. 1279. 1628. 2019
# >2喂料May_Sep 1901。 413. 1875. 1212. 1605. 2145. 2811. 3422
# >3于709年10月至12月迁移。 156. 698. 446. 598. 802. 1024. 1294
# >#…带缩写变量名¹pctl_97.5
# >
# >$场景_10$选项2
# >#A台:3×10
# >季节周期平均值sd中值pctl_2.5 pctl_25 pctl_75 pctl_9…¹pctl_99
# ><chr><chr><dbl><dbl><dbl><dbl
# >1号繁殖站,Jan_Apr 404。 286. 328. 70.2 211. 509. 1101. 2131
# >5月2日喂食第676步。 481. 547. 116. 351. 852. 1833. 3592
# >3迁移10月至次年12月252日。 181. 205. 43.2 131. 320. 686. 1393.
# >#…带缩写变量名¹pctl_97.5
# >
# >$场景_10$opt3
# >#A台:3×10
# >季节id周期平均sd中位数pctl_2.5 pctl_25 pctl_75 pctl_9
# ><chr><chr><dbl><dbl><dbl><dbl
# >1育种Jan_Apr 47.6 52.0 32.8 5.05 18.6 56.9 164。 535
# >2喂食5月_第79.8 87.6 54.3 8.36 31.2 95.8 279步。 944
# >3迁移10月_ 12月29日33.2 20.5 3.15 11.7 35.6 101。 373
# >#…带缩写变量名¹pctl_97.5
# >
# >
# >$场景_11
# >$场景_11$opt1
# >#A台:3×10
# >季节周期平均值sd中值pctl_2.5 pctl_25 pctl_75 pctl_9…¹pctl_99
# ><chr><chr><dbl><dbl><dbl><dbl
# >1个繁殖期,1月至570年1月。 121. 558. 366. 484. 647. 826. 979
# >2喂食5月_第981步。 211. 957. 624. 830. 1115. 1431. 1682
# >3迁移10月至次年12月347日。 76.0 344. 221. 293. 395. 515. 606.
# >#…带缩写变量名¹pctl_97.5
# >
# >$scenario_11$opt2
# >#A台:3×10
# >季节周期平均值sd中值pctl_2.5 pctl_25 pctl_75 pctl_9…¹pctl_99
# ><chr><chr><dbl><dbl><dbl><dbl
# >1育种Jan_Apr 24.7 40.3 11.5 0.787 4.98 29.1 111。 440
# >2个喂食May_Sep 42.5 69.2 19.8 1.35 8.56 48.3 189。 768
# >3迁移10月_ 12月15日24.2 7.02 0.486 3.02 17.4 67.2 263。
# >#…带缩写变量名¹pctl_97.5
# >
# >$场景_11$选项3
# >#A台:3×10
# >季节周期平均值sd中值pctl_2.5 pctl_25 pctl_75 pctl_9…¹pctl_99
# ><chr><chr><dbl><dbl><dbl><dbl
# >1育种Jan_Apr 3.52 11.1 0.933 0.0455 0.337 2.66 21.6 133。
# >2喂食5月_步骤6.04 19.0 1.57 0.0814 0.584 4.48 36.6 233。
# >3迁移10月_ 12月2.13 6.60 0.568 0.0283 0.207 1.57 13.3 75.9
# >#…带缩写变量名¹pctl_97.5
# >
# >
# >$场景_12
# >$场景_12$opt1
# >#A台:3×10
# >季节周期平均值sd中值pctl_2.5 pctl_25 pctl_75 pctl_9…¹pctl_99
# ><chr><chr><dbl><dbl><dbl><dbl
# >1135年1月1日繁殖。 241. 1111. 728. 964. 1287. 1645. 1949
# >2喂食1932年5月至12月。 415. 1885. 1229. 1635. 2197. 2819. 3315
# >3迁移10月至次年12月696日。 152. 689. 442. 587. 791. 1032. 1212
# >#…带缩写变量名¹pctl_97.5
# >
# >$场景_12$op2
# >#A台:3×10
# >季节周期平均值sd中值pctl_2.5 pctl_25 pctl_75 pctl_9…¹pctl_99
# ><chr><chr><dbl><dbl><dbl><dbl
# >1育种Jan_Apr 50.9 81.9 24.0 1.68 10.5 60.3 229。 887.
# >2喂食5月_步骤86.6 139。 41.0 2.85 17.8 99.1 384. 1534
# >3迁移Oct_Dec 31.2 49.5 14.8 1.04 6.37 36.3 139。 535
# >#…带缩写变量名¹pctl_97.5
# >
# >$场景_12$opt3
# >#A台:3×10
# >季节周期平均值sd中值pctl_2.5 pctl_25 pctl_75 pctl_9…¹pctl_99
# ><chr><chr><dbl><dbl><dbl><dbl
# >1育种Jan_Apr 7.17 22.2 1.95 0.0973 0.707 5.54 43.2 268。
# >2喂食5月_步骤12.2 37.7 3.25 0.172 1.22 9.21 73.2 464。
# >3迁移10月_Dec 4.38 13.4 1.19 0.0606 0.439 3.28 27.1 154。
# >#…带缩写变量名¹pctl_97.5
# >
# >
# >$场景_13
# >$场景_13$opt1
# >#A台:3×10
# >季节周期平均值sd中值pctl_2.5 pctl_25 pctl_75 pctl_9…¹pctl_99
# ><chr><chr><dbl><dbl><dbl><dbl
# >1育种Dec_Feb 21.5 4.95 21.0 13.1 18.0 24.4 32.0 39.7
# >2喂食Mar_Aug 60.9 13.9 59.8 37.2 50.8 69.7 91.1 112。
# >3迁移9月24日5.70 24.0 14.9 20.3 27.8 37.0 45.4
# >#…带缩写变量名¹pctl_97.5
# >
# >$场景_13$op2
# >#A台:3×10
# >季节周期平均值sd中值pctl_2.5 pctl_25 pctl_75 pctl_9…¹pctl_99
# ><chr><chr><dbl><dbl><dbl><dbl
# >1育种Dec_Feb 21.6 4.72 21.3 13.2 18.4 24.6 31.5 36.8
# >2喂食Mar_Aug 61.3 13.3 60.9 37.6 52.2 69.8 89.0 106。
# >3迁移Sep_Nov 24.6 5.48 24.3 14.8 20.9 27.9 36.1 42.8
# >#…带缩写变量名¹pctl_97.5
# >
# >$场景_13$opt3
# >#A台:3×10
# >季节周期平均值sd中值pctl_2.5 pctl_25 pctl_75 pctl_9…¹pctl_99
# ><chr><chr><dbl><dbl><dbl><dbl
# >1育种Dec_Feb 18.6 3.99 18.7 11.4 15.9 21.1 27.2 35.6
# >2喂食Mar_Aug 52.9 11.3 52.9 32.7 45.1 59.8 77.3 100。
# >3迁移Sep_Nov 21.2 4.68 21.1 13.1 18.0 24.0 31.4 41.3
# >#…带缩写变量名¹pctl_97.5
# >
# >
# >$场景_14
# >$场景_14$opt1
# >#A兼容:3×10
# >季节周期平均值sd中值pctl_2.5 pctl_25 pctl_75 pctl_9…¹pctl_99
# ><chr><chr><dbl><dbl><dbl><dbl
# >1育种Dec_Feb 43.5 10.0 42.6 26.6 36.5 49.5 64.8 80.5
# >2喂食Mar_Aug 121。 27.6 119. 73.7 101. 138. 181. 223
# >3迁移9月11日49.0 11.4 48.1 30.0 40.8 55.7 74.3 91.2
# >#…带缩写变量名¹pctl_97.5
# >
# >$场景_14$op2
# >#A台:3×10
# >季节周期平均值sd中值pctl_2.5 pctl_25 pctl_75 pctl_9…¹pctl_99
# ><chr><chr><dbl><dbl><dbl><dbl
# >1育种Dec_Feb 46.0 9.96 45.4 28.2 39.2 52.4 66.7 78.3
# >2喂食Mar_Aug 128。 27.6 127. 78.7 109. 145. 185. 219
# >3迁移9月8日11.5 51.3 31.3 44.2 58.8 75.8 89.4
# >#…带缩写变量名¹pctl_97.5
# >
# >$场景_14$选项3
# >#A台:3×10
# >季节周期平均值sd中值pctl_2.5 pctl_25 pctl_75 pctl_9…¹pctl_99
# ><chr><chr><dbl><dbl><dbl><dbl
# >1育种Dec_Feb 39.6 8.40 39.7 24.3 33.9 45.0 57.7 75.1
# >2喂食Mar_Aug 110。 23.3 110. 68.7 94.2 125. 160. 206.
# >3迁移9月7日9.76 44.6 27.7 38.1 50.5 65.7 85.9
# >#…带缩写变量名¹pctl_97.5
# >
# >
# >$场景_15
# >$场景_15$opt1
# >#A兼容:3×10
# >季节周期平均值sd中值pctl_2.5 pctl_25 pctl_75 pctl_9…¹pctl_99
# ><chr><chr><dbl><dbl><dbl><dbl
# >1育种Dec_Feb 30.8 6.55 30.3 20.0 26.2 34.7 45.5 53.8
# >2喂食Mar_Aug 88.3 18.8 86.9 57.3 75.0 99.5 130。 157
# >3迁移9月11日34.5 7.50 33.9 22.1 29.1 38.9 51.6 60.5
# >#…带缩写变量名¹pctl_97.5
# >
# >$场景_15$op2
# >#A台:3×10
# >季节周期平均值sd中值pctl_2.5 pctl_25 pctl_75 pctl_9…¹pctl_99
# ><chr><chr><dbl><dbl><dbl><dbl
# >1育种Dec_Feb 16.7 3.64 16.5 10.1 14.0 19.1 23.8 29.4
# >2喂食Mar_Aug 47.8 10.4 47.6 29.1 40.2 54.7 68.6 83.1
# >3迁移9月18日-6月4日-6月11日11.2 15.7 21.3 27.0 33.1
# >#…带缩写变量名¹pctl_97.5
# >
# >$场景_15$opt3
# >#A台:3×10
# >季节周期平均值sd中值pctl_2.5 pctl_25 pctl_75 pctl_9…¹pctl_99
# ><chr><chr><dbl><dbl><dbl><dbl
# >1育种Dec_Feb 15.4 3.59 15.4 9.16 12.9 17.5 22.1 31.7
# >2喂食Mar_Aug 44.1 10.3 44.1 26.2 37.0 50.2 62.9 91.6
# >3迁移9月11日17.2 4.09 17.3 10.2 14.5 19.5 25.1 37.2
# >#…带缩写变量名¹pctl_97.5
# >
# >
# >$场景_16
# >$场景_16$opt1
# >#A台:3×10
# >季节周期平均值sd中值pctl_2.5 pctl_25 pctl_75 pctl_9…¹pctl_99
# ><chr><chr><dbl><dbl><dbl><dbl
# >1育种Dec_Feb 62.3 13.3 61.3 40.4 53.0 70.3 92.2 109。
# >2喂食Mar_Aug 175。 37.3 172. 114. 149. 197. 258. 311
# >3迁移9月11日69.1 15.0 68.0 44.2 58.3 78.1 104。 121
# >#…带缩写变量名¹pctl_97.5
# >
# >$场景_16$op2
# >#A台:3×10
# >季节周期平均值sd中值pctl_2.5 pctl_25 pctl_75 pctl_9…¹pctl_99
# ><chr><chr><dbl><dbl><dbl><dbl><dbl><dbl><dbl><dbl><dbl>
# >1育种Dec_Feb 35.4 7.67 35.0 21.5 29.9 40.6 50.5 62.3
# >2喂食Mar_Aug 99.4 21.5 98.9 60.8 83.8 114。 143. 173.
# >3迁移Sep_Nov 39.2 8.67 39.1 23.6 33.1 44.8 56.8 69.4
# >#…带缩写变量名¹pctl_97.5
# >
# >$场景_16$opt3
# >#A台:3×10
# >季节周期平均值sd中值pctl_2.5 pctl_25 pctl_75 pctl_9…¹pctl_99
# ><chr><chr><dbl><dbl><dbl><dbl
# >1育种Dec_Feb 32.7 7.54 32.7 19.6 27.6 37.1 46.9 66.3
# >2喂食Mar_Aug 91.9 21.1 91.8 54.9 77.1 104。 131. 188.
# >3迁移9月11日36.3 8.51 36.3 21.5 30.6 40.9 52.7 77.3
# >#…带缩写变量名¹pctl_97.5
# >
# >
# >$场景_17
# >$场景_17$opt1
# >#A台:3×10
# >季节周期平均值sd中值pctl_2.5 pctl_25 pctl_75 pctl_9…¹pctl_99
# ><chr><chr><dbl><dbl><dbl><dbl
# >1育种Dec_Feb 27.1 6.46 26.8 16.4 22.5 31.3 40.4 51.6
# >2喂食Mar_Aug 81.3 19.3 80.0 49.3 67.2 93.3 122。 157
# >3迁移9月10日7.53 30.4 18.5 35.7 47.4 58.1
# >#…带缩写变量名¹pctl_97.5
# >
# >$场景_17$op2
# >#A台:3×10
# >季节周期平均值sd中值pctl_2.5 pctl_25 pctl_75 pctl_9…¹pctl_99
# ><chr><chr><dbl><dbl><dbl><dbl
# >1繁殖Dec_Feb 1.10 0.470 1.04 0.417 0.761 1.34 2.18 3.68
# >2喂食Mar_Aug 3.28 1.40 3.12 1.24 2.27 4.02 6.51 11.0
# >3迁移Sep_Nov 1.25 0.541 1.19 0.471 0.869 1.54 2.51 4.42
# >#…带缩写变量名¹pctl_97.5
# >
# >$场景_17$opt3
# >#A台:3×10
# >季节周期平均值sd中值pctl_2.5 pctl_25 pctl_75 pctl_9…¹pctl_99
# ><chr><chr><dbl><dbl><dbl><dbl
# >1育种Dec_Feb 0.774 0.408 0.725 0.283 0.520 0.957 1.56 4.85
# >2喂食Mar_Aug 2.32 1.21 2.18 0.855 1.56 2.88 4.64 14.1
# >3迁移Sep_Nov 0.884 0.467 0.840 0.333 0.593 1.10 1.78 5.34
# >#…带缩写变量名¹pctl_97.5
# >
# >
# >$场景_18
# >$场景_18$opt1
# >#A兼容:3×10
# >季节周期平均值sd中值pctl_2.5 pctl_25 pctl_75 pctl_9…¹pctl_99
# ><chr><chr><dbl><dbl><dbl><dbl
# >1繁殖Dec_Feb 54.9 13.1 54.2 33.2 45.5 63.3 81.7 104。
# >2喂食Mar_Aug 161。 38.2 159. 97.8 133. 185. 242. 310
# >3迁移9月22日15.1 61.0 37.2 51.0 71.7 95.0 117。
# >#…带缩写变量名¹pctl_97.5
# >
# >$场景_18$选项2
# >#A台:3×10
# >季节周期平均值sd中值pctl_2.5 pctl_25 pctl_75 pctl_9…¹pctl_99
# ><chr><chr><dbl><dbl><dbl><dbl
# >1育种Dec_Feb 2.33 0.992 2.22 0.889 1.62 2.84 4.63 7.73
# >2喂食Mar_Aug 6.84 2.90 6.50 2.61 4.73 8.38 13.6 22.7
# >3迁移9月11日2.64 1.13 2.52 0.998 1.83 3.24 5.29 9.20
# >#…带缩写变量名¹pctl_97.5
# >
# >$场景_18$opt3
# >#A台:3×10
# >季节周期平均值sd中值pctl_2.5 pctl_25 pctl_75 pctl_9…¹pctl_99
# ><chr><chr><dbl><dbl><dbl><dbl
# >1育种Dec_Feb 1.65 0.855 1.54 0.606 1.11 2.03 3.31 10.1
# >2喂食Mar_Aug 4.83 2.49 4.55 1.79 3.26 6.00 9.54 28.7
# >3迁移9月11日1.86 0.969 1.77 0.707 1.25 2.31 3.68 11.0
# >#…带缩写变量名¹pctl_97.5
图书馆( 斯托克实验室 ) # ------------------------------------------------------
# 使用任意参数值运行,以进行说明
# ------------------------------------------------------
# 设置要从中绘制的参数的数据帧
参数 <- 数据帧 ( 飞行_速度 = 13.1 , # 飞行速度(m/s)
车身_ lt = 0.85 , # 主体长度(m)
翼展 = 1.01 , # 翼展(m)
航班类型 = " 拍打 " , # 扑翼或滑翔飞行
避免_基本 = 0.989 , # 选项1和选项2的回避率
避免文本 = 0.981 , # 方案3和4的延长回避率
noct_活性 = 0.5 , # 白昼鸟不活动的比例
prop _ crh _ surv(比例_相对曲线) = 0.13 , # 碰撞风险高度的鸟类比例(仅选项1)
prop_upwind(上风) = 0.5 , # 逆风航班比例
转子转速 = 15 , # 转子速度(m/s)
转子半径 = 120 , # 水轮机半径(m)
叶片_宽度 = 5 , # 涡轮叶片最厚点的宽度(m)
叶片缝 = 15 , # 平均半径节距(弧度)
n个叶片 = 三 , # 每个涡轮机的叶片总数
轮毂高度 = 150 , # HAT上方轮毂高度(m)
n台涡轮机 = 100 , # 风电场中的涡轮机数量
wf_宽度 = 52 , # 风电场最长部分的宽度
wf_纬度 = 56 , # 风电场质心纬度
潮汐_偏移 = 2.5 , # 风电场HAT的平均潮汐偏移
lrg_错误_错误 = 真的 # 应用大数组校正? ) # 每月鸟类密度
b_dens(密度) <- 数据帧 ( 月 = abb月 , 洞穴 = 运行( 12 , 0.8 , 1.5 ) ) # Johnston等人的飞行高度分布
发电机_fhd_dat <- 约翰斯顿_航班_高度_SOSS % > % dplyr公司 :: 过滤器( 变量 == " Gannet.est公司 " ) % > % dplyr公司 :: 选择( 高度 , 支柱 ) # 风电场的月运行时间
turb_oper(涡轮机操作) <- 数据帧 ( 月 = abb月 , 比例(_O) = 运行( 12 , 0.5 , 0.8 ) ) 斯托克实验室 :: 波段_厘米( 模型选项(_O) = c(c)( 1 , 2 , 三 ), 飞行速度 = 参数 $ 飞行_速度 , 车身_ lt = 参数 $ 车身_ lt , 翼展 = 参数 $ 翼展 , 航班类型 = 参数 $ 航班类型 , 避免_rt_basic = 参数 $ 避免_基本 , 避免文本 = 参数 $ 避免文本 , 无放射性 = 参数 $ noct_活性 , prop _ crh _ surv(比例_相对曲线) = 参数 $ prop _ crh _ surv(比例_相对曲线) , 密度_月 = b_dens(密度) , prop_upwind(上风) = 参数 $ prop_upwind(上风) , 发电机(_fhd) = 发电机_fhd_dat , 站点_fhd = 无效的 , # 仅选项4
转子转速 = 参数 $ 转子转速 , 转子半径 = 参数 $ 转子半径 , 叶片_宽度 = 参数 $ 叶片_宽度 , 叶片缝 = 参数 $ 叶片间距 , n个叶片 = 参数 $ n个叶片 , 轮毂_高度 = 参数 $ 轮毂_高度 , 弦乐_专业 = 弦_prof_5MW , n_涡轮机 = 参数 $ n台涡轮机 , turb_oper_月 = turb_oper(涡轮机操作) , wf_width(宽度) = 参数 $ wf_width(宽度) , wf_纬度 = 参数 $ wf_纬度 , 潮汐偏移量 = 参数 $ 潮汐_偏移 , lrg_错误_错误 = 参数 $ lrg_错误_错误 ) # >#A tibble:12×4
# >月opt1 opt2 opt3
# ><chr><dbl><dbl><dbl
# >98.6年1月1日25.4 10.4
# >2月2日53.2 13.7 5.61
# >83.7年3月3日21.6 8.83
# >4月4日61.4 15.8 6.48
# >117年5月5日。 30.2 12.4
# >91.7年6月6日23.6 9.67
# >114年7月7日。 29.3 12.0
# >8月8日91.1 23.5 9.61
# >2006年9月9日。 27.3 11.2
# >64.1年10月10日16.5 6.77
# >11月11日60.7 15.7 6.40
# >12月12日73.4 18.9 7.75