- 博客(19)
- 资源 (1)
- 问答 (1)
- 收藏
- 关注
RSS订阅原创 R语言如何最简单的写脚本并添加帮助文档
1. 默认的参数commandArgs超简单1.1 脚本示例head.Rargs = commandArgs(TRUE)if(length(args) != 2){ cat("运行命令方式：Rscript head.R dat.csv 5\n\thead.R 为脚本\n\tdat.csv 为数据\n\t5 为行数\n") quit("no")}else{ dd ...
2019-11-26 22:50:07
2761
原创 数据分析中转化数据对结果的影响
1. 问题提出在模型假定中，数据需要符合正态分布，在动物模型分析中，如果对数据进行转化，对育种值和方差组分变化的情况，进行测试。2. 示例正常数据library(asreml)library(learnasreml)data("animalmodel.dat")data("animalmodel.ped")dat = animalmodel.datped = animalmod...
2019-11-26 07:58:43
941
原创 学习R语言中什么时候感觉要突飞猛进了？
学习生涯回想我自己的学习R语言的经历，总结如下：1. 敲门第一级：跑别人的代码2. 入门第二级：修改别人的代码，安装包3. 晋级第三级：查看包中帮助文档，自己写代码4. 突飞猛进第四级：根据自己的需求，将代码写为函数或者脚本，方便调用5. 修成正果第五级：编写自己的R包，将函数进行打包6. 蓦然回首第六级：把自己的R包，放到CRAN，github，博...
2019-11-23 17:36:05
508
原创 cp mv 使用通配符时 对一些文件进行排除
问题描述现有几个文件和一个文件夹，想把这几个文件，放到文件夹中。常规操作(base) [dengfei@ny01 a]$ mv * testmv: 无法将目录"test" 移动至自身的子目录"test/test" 下结果也成功了，但是出现报错：mv: 无法将目录"test" 移动至自身的子目录"test/test" 下将test排除mv `ls|grep -v "test"` ...
2019-11-21 08:16:22
2226
原创 月亮与六便士--未选择的路
最近看了毛姆的《月亮与六便士》，为毛姆的洞察力和幽默所折服。看完这本书，让我想起了这首诗歌：《未选择的路》作者 美国 罗伯特·弗罗斯特黄色的林子里有两条路，很遗憾我无法同时选择两者身在旅途的我久久站立对着其中一条极目眺望直到它蜿蜒拐进远处的树丛。我选择了另外的一条，天经地义，也许更为诱人因为它充满荆棘，需要开拓；然而这样的路过并未引起太大的改变。那天清晨这两条小路一起...
2019-11-18 23:42:21
363
原创 R语言中多条件排序
R语言中，如何对数据框的数据根据某个条件进行排序呢？如何根据多条件进行排序呢，类似Excel中的排序效果：1. 示例数据R语言中鸢尾花的数据，数据有五列：> names(iris)[1] "Sepal.Length" "Sepal.Width" "Petal.Length" "Petal.Width" "Species" 前五行数据预览：> head(iris) ...
2019-11-13 19:53:28
18917
2
原创 pandas 数据框 将一列数据切割变为两列
1. 原始数据想把Sex一列，根据下划线且分为两列数据In [193]: user_info Out[193]: age city sex birthname ...
2019-11-13 19:52:37
1305
原创 基因组选择的几个概念
1. 基因组选择:(Genomic selection)中文：基因组选择利用覆盖全基因组的高密度SNP标记, 结合表型记录或系谱记录对个体育种值进行估计, 其假定这些标记中至少有一个标记与所有控制性状的QTL处于连锁不平衡状态.英文：Genomic selection is a form of marker-assisted selection in which genetic marke...
2019-11-13 19:48:54
5499
2
原创 多个性状 将多个性状画一张箱线图上R实战
示例数据数据描述：第一列是品种其它列为观测值，y1,y2…目的：相对每一个性状做箱线图，同时显示在一个图上面，横坐标为性状。> head(dd) Cul y1 y2 y3 y4 y5 y6 y7 y81 80001 0.334 0.405 0.358 29 130 239 420 6302 80002 0.348 0.393 0.365 24...
2019-11-13 12:01:25
3909
4
原创 top c命令查看内存
命令top cPID, 任务的ID号，可以用于killUSER，用户名，这里dengfei是我的用户名PR ：优先级NI ：负值表示高优先级，正值表示低优先级。VIRT ：虚拟内存RES ： 真实内存SHR ：共享内存S ：进程状态 D=不可中断的睡眠状态； R=运行； S=睡眠 ；T=跟踪/停止； Z=僵尸进程虚拟内存，真实内存，共享内存的关系VIRT=RES+SH...
2019-11-12 08:28:18
2890
原创 numpy将数据转化为缺失值
介绍numpy如何将某个条件的数据，转化为缺失值。import numpy as npdd = np.random.randn(3,5)dd将dd中小于0的数字，转化为缺失值：# 将dd小于0的赋值为缺失值dd[dd<0] = Nonedd...
2019-11-12 08:27:58
1805
原创 R语言中 判断对象是否存在
R语言中如何判断对象是否存在描述：比如对象a，在没有赋值之前，它是不存在的，想要通过一个函数，判断它是否存在。解决方案：比如tt，这个对象不存在> ttError: object 'tt' not found使用函数exists判断> exists("tt")[1] FALSE如果给tt赋值，再进行判断：> tt = 1> exists("tt"...
2019-11-12 08:27:29
10879
原创 一道数量遗传学题：如何计算育种值
题目设动物个体效应为随机遗传效应（a），日粮、性别和畜舍为固定环境效应（b），背膘厚的遗传力为0.4，请完成以下工作：1，建立背膘厚的线性模型2，写出模型的一般形式和矩阵形式3，写出混合线性模型方程组的各组分成分4，获得的估计值具有哪些特点5，不同日粮和性别的效应值是多少6，个体育种值是多少，是否和表型值排序一致？说明理由处理思路线性模型已经很清楚：固定因子：日粮，性别...
2019-11-08 08:29:47
4768
原创 提升工作效率的神器：电脑版番茄工作软件
番茄工作法很好，收集app居多，电脑版很少。这里介绍一个非常简洁，非常好用的电脑版：番茄钟10，介绍如下：1. 下载番茄钟10打开微软应用商店，搜索番茄钟10，点击安装。2. 使用番茄钟10软件点击三角符号，开始运行，进入专注模式：3. 番茄钟10 设置时间点击设置，可以设置专注时间和休息时间4. 提醒专注时间结束后，会弹窗提醒时间到，进入休息周期。5. 快速启动将番...
2019-11-08 08:16:58
1328
原创 python list 转化为矩阵报错：AttributeError: 'list' object has no attribute 'reshape'
报错类型>>> b['1', '0.1', '1', '0.1']>>> b.reshape(2,2)Traceback (most recent call last): File "<stdin>", line 1, in <module>AttributeError: 'list' object has no attri...
2019-11-07 21:55:38
31317
1
原创 Fortran 搭配vscode安装及配置
1. Fortran安装及配置1.1 安装MinGWMinGW-w64 - for 32 and 64 bit Windows下载安装包。双击打开安装程序，点击Next：点击Next:一直默认安装：1.2 配置环境变量我的路径为：C:\Program Files (x86)\mingw-w64\i686-8.1.0-posix-dwarf-rt_v6-rev0\mingw...
2019-11-07 08:22:19
6309
1
原创 Python中输出井号#
代码中的#是包含在字符串中的，字符串直到遇到下一个"为止，字符串里的#只是当做一个字母而不是注释处理。代码：print("hello # world")结果：
2019-11-07 08:19:11
3494
原创 根据日期 提取场年季效应
1， 提取准者根据年-月-日， 提取年+季节根据日期的年，提取年根据日期的月，提取季节将两者合并2. 示例> head(test) ID Ceding_dat1 1 2015-07-102 2 2015-07-103 3 2015-08-124 4 2016-02-285 5 2016-04-266 6 2016-04-26对数据Ceding_...
2019-11-06 07:55:56
546
TA创建的收藏夹 TA关注的收藏夹
TA关注的人