• Welcome to LiuJason's Blog!
R语言中使用readOGR读取shp文件坐标转换|经纬度值很大

R语言中使用readOGR读取shp文件坐标转换|经纬度值很大

我在绘制比利时海岸线的研究地点地图时,下载了geo.be的shp文件,使用如下代码导入的时候出现了问题:经纬度值非常大be_map <- readOGR("SiteMap/resources/AD_3_District.shp")这是由于使用的坐标系不同导致的,因此需要做转换:be_map <- spTransform(be_map, ……

[已解决]libudunits2.so was not found安装sf时报错|CentOS

[已解决]libudunits2.so was not found安装sf时报错|CentOS

安装sf的时候出现如下错误:checking udunits2.h usability... nochecking udunits2.h presence... nochecking for udunits2.h... nochecking udunits2/udunits2.h usability... nochecking udunits2/u……

Supporting Materials (DO and DTA files) of Blossfeld's book Event History Analysis with Stata|Attachment

Supporting Materials (DO and DTA files) of Blossfeld's book Event History Analysis with Stata|Attachment

Authors:Hans-Peter Blossfeld, Katrin Golsch, and Götz RohwerPublisher: Lawrence ErlbaumCopyright: 2007ISBN-13: 978-0-8058-6047-4Download Here:EventHistoryAnalysisWithStata.z……

[已解决]

[已解决]"proj_api.h presence and usability... no"安装proj-devel——R语言中GIS包rgdal的依赖包|CentOS7为例

首先看这里:https://www.liujason.com/article/152.html然后,ssh里安装gdal-devel, proj-devel, proj-nad, proj-epsglibgdal.so问题:cp /usr/local/lib/libgdal.so.20* /usr/lib64/libkmldom.so问题:wget……

R语言中在表格中显示Tukey HSD事后检验的分组结果|和SPSS一样可视化ANOVA方差分析

R语言中在表格中显示Tukey HSD事后检验的分组结果|和SPSS一样可视化ANOVA方差分析

在多层次方差分析之后,我们可能需要看看到底是哪些处理组之间有显著差异,或者虽然ANOVA结果是显著的但是我们还要看看哪些处理组之间是不显著的,这种时候我们就需要TukeyHSD事后检验了。用过SPSS的人应该知道,SPSS的事后检验结果非常的人类友好——可读性很高,会把每个层次的处理告诉你,但是R自带的结果就非常不友好了...我拿同一组数据在R和SPSS里跑……

Centos7上的SF安装——GIS依赖包

Centos7上的SF安装——GIS依赖包

依赖列表gdal> = 2.1#用于数据结构geos> = 3.5#用于数据操作pro4j#用于投影udunits#为单位的PostgreSQL [可选]lwgeom [可选]#对于geohashCentos7默认仓库为您提供gdal1.1,geos 2.x. 但是,sf :: st_voronoi需要geos> =……

[已解决]

[已解决]"checking for gdal-config... no" 安装gdal2.x——R语言中GIS包rgdal的依赖包|CentOS7为例

在编译R环境的时候,最讨厌的额就是一些特殊的依赖包,只能通过ssh安装,而且还经常找不到这些包或者是版本不匹配...这次的gdal 2.x就是一个例子,浪费了我很长的时间去找相关的安装方法,这里记一下笔记,下次碰到了的话就方便了。这里还是强烈建议,安装各种包之前看一看包的介绍:packages/rgdal。 有很明确的说明需要什么样子的环境:for bu……

使用ggplot2中的分面参数facet对数据进行分类作图/分页/一页多图

使用ggplot2中的分面参数facet对数据进行分类作图/分页/一页多图

最近在整理文章用到的分析的代码,准备发表的时候一并上传给大家作参考,所以要仔细检查有没有错误或者是可以优化代码的地方。这里给大家介绍一下ggplot2里一个超级好用的参数:facet。用过R的应该都知道,普通的绘图我们可以用par(nfrow())来进行分页,但是这一方法在ggplot2中并不适用,事实上ggplot2的facet更加好用一些,可以自动帮你……

使用R语言插入缺失数据-平均值法(附代码)

使用R语言插入缺失数据-平均值法(附代码)

做实验嘛总会莫名其妙少掉一些样本,但是样本少了文章还是要写啊!那么缺失了数据对分析有什么影响呢?我觉得最重要的影响应该是权重的变化,这里用我这次的实验数据做例子吧:本实验是研究9个不同来源地(纬度梯度下)的植物凋落物分解速率的差异,同时还有3种不同的凋落物袋看分解者和来源地的交互作用,具体可以看这个poster:Poster (draft) of [R……