科学网

 找回密码
  注册
置顶 · 提供数据分析服务-本帖长期有效(5月更新)
2019-5-30 16:01
本人从事林业统计分析,主编R与ASReml软件的教材2部,主译学术专著1部。熟悉植物遗传评估,擅长R软件与ASReml软件。熟悉R编程,包括常见的基础统计(方差分析、协方差分析、t检验等)和高级统计(主成分分析、聚类分析、判别分析等);同时,也善于使用R软件进行图形的绘制,包括高级程序包ggplot2和lattice;此外,本 ...
个人分类: R|4498 次阅读|没有评论
AAfun4包的所有函数及其功能
林元震 2019-7-26 14:15
--------------------- 作者:yzhlinscau 来源:CSDN 原文: https://blog.csdn.net/yzhlinscau/article/details/97385205 版权声明:本文为博主原创文章,转载请附上博文链接! AAfun4包是本人根据多年的遗传评估经验,基于asreml包开发的辅助程序包,用于增加一些辅助函数,实现一些特定的功能 ...
个人分类: ASReml|2114 次阅读|没有评论
AAfun4包的新函数
林元震 2019-7-21 11:36
https://blog.csdn.net/yzhlinscau/article/details/96703282 主要针对原有的**asreml.batch()**,增加了两个函数:**repin.batch**和**update.batch**,用于计算新的遗传参数,或者重新运行新的批量分析。最新版的asreml.batch(),基本上能涵盖常见所有的可用模型,比如家系模型,个体模型,空间模型,基因 ...
个人分类: ASReml|1814 次阅读|没有评论
如何批量载入R程序包
林元震 2019-7-19 11:54
有时,我们需要载入一系列R程序包,如果仅使用library()来载入,需要写多次。下述代码,基本2行搞定。 pkgs-c('asreml','AAfun4','ggplot2') lapply(pkgs,library,character.only=TRUE)
个人分类: R|3337 次阅读|没有评论
AAfun4包的部分更新
林元震 2019-7-19 11:40
因科学网对程序代码支持差,具体博文请见下述网页: https://blog.csdn.net/yzhlinscau/article/details/96477198 注:最新的AAfun4包 不再免费提供。
个人分类: ASReml|1446 次阅读|没有评论
用R语言玩玩茅台股票
林元震 2019-7-2 15:22
本人不炒股,只是对R语言的分析一直感兴趣,加上偶尔会看看茅台酒的一些情况。最近,刚好看了一本用R语言量化投资的书,于是借用此书程序代码,玩玩茅台股票。仅供娱乐使用!切勿当真! 因科学网对代码支持得差,感兴趣的,全文请见下述网页: https://blog.csdn.net/yzhlinscau/article/details/94444992 &nb ...
个人分类: R|2521 次阅读|没有评论
R语言的使用体会
林元震 2019-7-2 10:10
从2012年接触R语言至今,已过去了7年。R语言一路在发展,版本也在不断地更新,程序包也源源不断地涌向。但最近2年,我发现,R语言使用起来有些麻烦,尤其当你更新R语言版本,或者买了新电脑外,重新安装R语言及其R程序包,有时就会发现程序包安装不了!之前能够使用的R包,不再工作了!此外,不同版本R包,有时与R语言不 ...
个人分类: R|2757 次阅读|没有评论
一步法GBLUP(SS-GBLUP)的实施
林元震 2019-7-1 21:04
一步法GBLUP(SS-GBLUP)可综合利用基因分型数据和谱系文件的全部信息,进行基因组BLUP,获得更多可靠结果。AAfun4包最新版里添加了函数AGH.inv(),用于计算H矩阵的逆矩阵,可用于SS-GBLUP。 简单示范如下: library(AAfun4) data(ugped)#未基因分型的谱系文件 data(gped)&nbs ...
个人分类: ASReml|5285 次阅读|没有评论
林木基因型与环境互作的研究方法及其应用--综述
林元震 2019-6-9 08:37
林业科学 2019 , Vol. 55 Issue (5) : 142-151 DOI : 10.11707/j.1001-7488.20190516 摘要:我国是全球第一大木材进口国和第二大木材消费国,木材对外依存度已连续多年超过50%。然而,我国每公顷森林年均生长量约为林业发达国家水平的一半,这说明 ...
个人分类: 数量遗传学|2773 次阅读|没有评论
生态R包spatstat的部分用法(三)--地形分析
林元震 2019-1-7 11:20
继续之前的帖子,下述代码主要是针对地形分析。 1. 地形数据读取 esdf-readxl::read_excel(data2.xlsx,sheet=1,skip=1) 2.地形因子生成 im 形式 ##选出海拔的数据 elev=esdf$meanelev ##将海拔数据转换为30行50列的矩阵 e=matrix(data=elev,nrow=30,ncol=50) ...
个人分类: R|5070 次阅读|没有评论

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

GMT+8, 2024-9-26 02:59

Powered by ScienceNet.cn

Copyright © 2007- 中国科学报社

返回顶部