如果CRAN上找不到包，那怎么办？

热度 1 赖江山 2017-5-3 13:37

#如果CRAN包找不到的包，可以使用githubinstall函数来找包的资源 install.packages('githubinstall') #已发布至CRAN library(githubinstall) githubinstall('AnomalyDetection')#输入任何包的名称即可在web上找到这些包的位置

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

已知经纬度直接算地理距离的函数

热度 1 赖江山 2017-4-11 08:20

install.packages(geosphere) library(geosphere) xy - rbind(c(0,0),c(90,90),c(10,10),c(-120,-45)) distm(xy) xy2 - rbind(c(0,0),c(10,-10)) distm(xy, xy2)

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

用R直接从全球生物多样性信息网络(GBIF)获取物种地理分布数据

热度 6 赖江山 2017-4-10 14:23

#install.packages(c('raster','rgdal','dismo','rJava')) #loadsthedismo(distributionmodeling)library #usinganexamplefilethatisinstalledwiththedismopackag ...

个人分类: RDA|19152 次阅读|8 个评论 热度 6

R里面读取netcdf（.nc)格式数据的函数

赖江山 2017-4-10 13:48

# NetCDF 全称为network Common Data Format， 网络通用数据格式是由美国大学大气研究协#（University Corporation for Atmospheric Research，UCAR) 的Unidata项目科学家针对科学数据的特点开发的，是一种面向 #数组型并适于网络共享的数据的描述和编码标准。目前，N etCDF广泛应用于大气科学、水文、海洋学、 ...

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

将大样地划分为任意矩形的小样地的函数

赖江山 2017-4-9 06:58

gxgy.to.rectindex=function (gx, gy, gridsize = c(20,20), plotdim = c(1000, 500)) { badgxgy = (gx 0 | gy 0 | gx = plotdim | gy = plotdim | is.na(gx) | is.na(gy)) colno = 1 + floor(gx/gridsize ) rowno = 1 + flo ...

4085 次阅读|没有评论

只知道相关系数，如何算显著性p值

赖江山 2017-4-6 18:43

#从别的地方的相关系数，如果知道样本量，在R里面输入下面的代码可以获得显著性，r是相关系数，n是样本量 #如果 r是负值， pt(r*sqrt((n-2)/(1-r^2)),n-2)*2 #如果 r是正值 (1-pt(r*sqrt((n-2)/(1-r^2)),n-2))*2

13925 次阅读|没有评论

画线性回归线的置信区间的线

赖江山 2017-4-6 17:56

x - rnorm(15) y - x + rnorm(15) predict(lm(y ~ x)) new - data.frame(x = seq(-3, 3, 0.5)) predict(lm(y ~ x), new, se.fit = TRUE) pred.w.plim - predict(lm(y ~ x), new, interval = prediction)#拟合值置信区间 matplot(new$x, pred.w.plim,lty = c(1,2,2), type = l, ylab = pr ...

8455 次阅读|没有评论

设定回归线的长度

热度 1 赖江山 2017-4-6 10:50

#abline是画回归线的函数，但是划线的边际延伸到图形边界，有时候不需要这么长回归线，可以clip来控制 x - rnorm(15) y - rnorm(15) plot(y~x) clip(min(x), max(x), min(y), max(y)) abline(lm(y~x))

4071 次阅读|1 个评论 热度 1

如何去掉数据框带某一字符的行的R代码

赖江山 2017-3-30 21:32

#最近处理一份植物名录，其中要将带有“科”和“属”的行去掉，只保留一个种的行，用grep()函数解决 badong=read.csv(巴东2.csv,header=T) head(badong) #去掉科和属，因为最后有一个总的名录可以匹配 a=badong$中文名 n=grep(科,a) #选择带科的行的index m=grep(属,a) #选择带属的行 ...

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

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