科学网

 找回密码
  注册
搜索
热搜: 科学 论文
搜索
Use R: 画泰勒图 Taylor Diagram,模型比较一目了然
庄伟 2020-10-14 22:54
IPCC推荐的模型评价方法: 在做模型相关的工作时,我们经常会需要比较模拟效果。当使用多个模型时,看那个模型的模拟效果更好些,模拟误差更小,仅仅画散点图进行比较是不直观的。Taylor Diagram 泰勒图是一种可以直观进行模型比较的方法。它巧妙利用三角函数几何原理,把三种评价模型的指标在一张图上展现。IPCC的报告 ...
个人分类: R|753 次阅读|没有评论
USE R 画出散点密度图
热度 1 庄伟 2019-11-6 16:24
有三种方法 一种如下: 可以用第一个文件中的代码生成图例。这样方法用R里面基本的画图功能,有些格式要自己调整。图例不能自己生成,需要出图后,自己修改。 r_scatter_density_plot.R 这是第二种方法 借助ggplot 不过这个示例中,如何画出图例,需要思考一下,有待解决 r_scat ...
个人分类: R|1499 次阅读|2 个评论 热度 1
Penman-Monteith模型模拟蒸散R代码示例
庄伟 2019-10-28 09:34
运行代码在下方zip附件中,请下载:) PM_HowardSprings_20190918.zip
个人分类: R|1158 次阅读|没有评论
Century 模型 R代码分享
庄伟 2019-8-30 13:45
Century 是传统的土壤碳过程模型,这里分享本人最近找到的三个相关版本: 1. Fortran 版本(linux) 该版本的下载地址: https://github.com/bsulman/INTERFACE-model-experiment-synthesis DAYCENT.zip Sulman et al (2018)发表了多个土壤碳模型比较的文章,并分享了部分模型的代码,century的源代 ...
个人分类: R|1980 次阅读|没有评论
USE R: SCEUA优化方法进行Shuttleworth-Wallace 模型的参数优化
庄伟 2019-7-9 17:02
sceua_SW.zip 分享一个SCEUA进行参数优化的例子,模型是Shuttleworth-Wallace 模型,其中的冠层导度子模型采用自Leuning et al. (2008) WRR. 用到了Ozflux一个站点的数据,数据为白天的日均值,优化的参数是最大气孔导度gmax,和土壤表面阻力有关的参数b3。感兴趣的老师同学可以看看。这里写的R代码是在j ...
个人分类: R|1464 次阅读|没有评论
ggplot画图时间坐标轴显示英文格式的方法
庄伟 2019-4-17 15:14
发表在国际期刊的论文,一般要求图中的坐标都是英文。用R画图时,输出时间坐标轴有时会用中文时间来显示,用下面的命令,可以把输出的语言变成英语: Sys.setlocale(LC_TIME, English)
个人分类: R|2060 次阅读|没有评论
jupyter lab:交互式的编程环境,它可以帮助我们完成一个完成的实验报告
庄伟 2019-3-9 16:45
Jupyter lab 是Jupter Notebook的升级版。 它继承了Jupyter Notebook的优点:将可执行的代码、描述性文本(Markdown)、方程(LaTex)、图片、交互式可视化和运行输出结果等整合显示到一个笔记本文件中。 Jupyter lab可以配置R、Python等运行环境。 以下是官网对Jupyter lab的介绍视频 https://jupyterlab.read ...
个人分类: R|2093 次阅读|没有评论
[转载]如何在jupyterlab 中运行R (windows)效果可肩比RSTUDIO
庄伟 2019-3-7 13:18
如何在jupyterlab 中运行R Installing RKernel on Windows JupyterLab 本文转载自 https://qiita.com/crowdy/items/3738ebe1edb95839753f 准备工作: 安装anaconda, 下载地址: https://www.anaconda.com/distribution/ 安装完毕后,打开anaconda Prompt 输入: conda install -c conda-f ...
个人分类: R|3433 次阅读|2 个评论
R中multisensi和sensitivity 全局敏感性分析的运行代码示例
庄伟 2019-3-4 10:17
全局敏感性分析(Global Sensitivity Analysis)的例子见: https://cran.r-project.org/web/packages/multisensi/vignettes/multisensi-vignette.pdf 附上一个来自multisensi包的全局敏感性分析的例子: 使用的数据和模型函数是来自multisensi,使用的敏感性分析的fast99来自sensitivity. 以下为code: #&n ...
个人分类: R|2189 次阅读|没有评论
USE R: 使用sensitivity包中的fast99命令进行局部敏感性分析
热度 1 庄伟 2018-9-11 17:14
R中进行敏感性分析的包sensitivity中的fast的全称是Extended Fourier Amplitude Sensitivity Test。 按照R自带的HELP的描述:Description fast99 implements the so-called extended-FAST method (Saltelli et al. 1999). This method allows the estimation of first order and total Sobol' indices for all ...
个人分类: R|3185 次阅读|3 个评论 热度 1

Archiver|手机版|科学网 ( 京ICP备07017567号-12 )

GMT+8, 2020-11-29 10:37

Powered by ScienceNet.cn

Copyright © 2007- 中国科学报社

返回顶部