分段回归的代码

赖江山 2019-3-22 19:26

#分段回归的代码 x - c(1:10, 13:22) y - numeric(20) ## Create first segment y - 20:11 + rnorm(10, 0, 1.5) ## Create second segment y - seq(11, 15, len=10) + rnorm(10, 0, 1.5) ## Plot it library(segmented) lin.mod - lm(y~x) #psi starting values for the breakpoints to b ...

个人分类: RDA|5697 次阅读|没有评论

关于批量读入数据文件并组合为一个文件的代码

赖江山 2019-3-8 09:23

在数据处理过程，我们经常要将很多独立的文件拼接起来，然后在组合到一个文件里面！其实干这个事情简单一个循环就可以搞定！ #批量读入数据并合并为一个文件 #在当前工作目录是下新建一个文件夹叫myname #把要读入的数据存储在myname文件夹里面，这里只能是1层不能再有下一层文件夹 #先用dir()函数获取目录中的 ...

个人分类: RDA|2407 次阅读|没有评论

在R里面演奏生日歌

热度 1 赖江山 2019-3-5 19:35

今天为了想知道在R里面是否能作曲演奏？在网上找了代码，果然可以，大家可以尝试一下，把这段代码拷贝到R里面运行。 在谱写生日歌.r

个人分类: RDA|3669 次阅读|1 个评论 热度 1

R语言在生态学研究当中的应用分析

赖江山 2019-1-12 09:38

R 语言在生态学研究当中的应用分析 编写代码处理数据已经成为很多生态学者日常的工作之一。能够用于编程的语言很多，近年来， R 语言以其灵活、开源、易于掌握、免费等诸多优点，在生态学研究各领域迅速传播并赢得广大研究者的青睐。 为了研究 R 语言在生态学研究论文里面使用情况，我们 ...

个人分类: RDA|9226 次阅读|没有评论

NMDS的是否有解释率？

热度 1 赖江山 2018-11-13 22:21

最近老有人问NMDS是否有排序轴是否有解释率？好像有些文章，特别是SBB的文章经常会有文章提供NMDS前几轴解释量（%)，根据我的理解 ，NMDS是非参数的方法，不应该有解释量，不知道他们这些NMDS的解释率是怎么算出来的？我估计是用主坐标分析（PCoA)的解释量来作为NMDS的解释量，但这种张冠李戴的做法可能是不 ...

个人分类: RDA|10890 次阅读|2 个评论 热度 1

关于排序图标识的虚化代码

赖江山 2018-10-13 21:55

#有时候排序图的中物种或样方过多，造成标识重叠的情况，可以用vegan包中的orditorp将标识虚化 ## A cluttered ordination plot : library(vegan) data(BCI) mod - cca(BCI) plot(mod, dis=sp, type=t) # Now with orditorp and abbreviated species names cnam - make.cepnames(names( ...

个人分类: RDA|3370 次阅读|没有评论

本页有 4 篇博文因作者的隐私设置或未通过审核而隐藏