原创 育種值 表型值 迴歸係數 相關係數 遺傳力之間的關係
假定表型值由均值+育種值+殘差 yi=μ+ai+ϵi y_i = \mu + a_i + \epsilon_iyi=μ+ai+ϵi 表型值 VS 育種值 他們之間的相關係數 cor(y,a)=Cov(y,a)Var(y)∗Var
2019-04-17 23:25:44
11
原创 R語言中怎麼講矩陣變爲三元組(即三列: row, col, value)
bing上搜索: 矩陣怎麼變爲三元組形式, 在bing上搜索結果: 各種C語言, 完全找不到R語言的實現方法. R語言寫循環很慢, 就用現成的包吧 Matrix Matrix是一個很強大的包, 它可以將矩陣轉化爲稀疏矩陣, 將0變爲.
2019-04-17 23:25:43
1
原创 R 語言如何提取年份-月份-季節-天
如何提取季節，月份，使用lubridate包 年份： year(datae) 月份： month(datae) 日期： day(datae) 季節： quarter(datae) 代碼 library(lubridate) d<-c(
2019-04-17 23:25:43
102
原创 R語言中對數據框的列名重命名
報錯類型 Error: All arguments must be named plyr中的rename和dplyr中的rename用法是不同的. plyr::rename rename(data, c(old=new)) dplyr
2019-04-17 23:25:43
24
原创 Julia數據框轉化爲矩陣
Julia構建數據框 using DataFrames t1 = DataFrame(y1=collect(1:3),y2=collect(4:6),y3=collect(7:9)) 結果: Main> t1 = DataFrame
2019-04-17 23:25:43
8
原创 基因組選擇中如何清洗基因組數據
打油詩 題目：plink666 不用不知道 一用忘不掉 朝思暮也想 真的好奇妙 1. 背景 一個朋友之前問過這個問題，問題可以分爲： 基因組選擇中如何將ATCG的數據轉化爲012的形式 如何進行基因組數據的過濾篩選 有沒有示例代碼可
2019-04-17 23:25:43
43
原创 如何高效的在服務器和本地進行上傳和下載文件
昨天, 師弟告訴我可以在xshell中使用sz進行下載, 想要上傳的話用rz就行了. 然後我竟然沒有聽過. 學習最好的方法就是寫一篇博客, 比如這篇. 1. 從putty到xshell 最開始， 是使用putty進行的服務器連接，痛點在
2019-04-17 00:04:25
2
原创 如何高效的在服務器和本地進行上傳和下載文件
昨天, 師弟告訴我可以在xshell中使用sz進行下載, 想要上傳的話用rz就行了. 然後我竟然沒有聽過. 學習最好的方法就是寫一篇博客, 比如這篇. 1. 從putty到xshell 最開始， 是使用putty進行的服務器連接，痛點
2019-04-17 00:04:25
原创 反對996爲何是一場鬧劇
概念定義 PSD是Poor, Smart and Deep desire to become rich 的縮寫，是對一種心理狀態的形容，所謂PSD是形容那些貧窮，但是很聰明，有遠大理想抱負，有成爲富人和上流社會的人的強烈慾望的人.
2019-04-15 23:57:32
原创 tidyverse安裝報錯/haven/libs/haven.so: undefined symbol: libiconv
library(tidyverse) Error: package or namespace load failed for ‘tidyverse’ in dyn.load(file, DLLpath = DLLpath, ...):
2019-04-15 23:57:22
原创 R語言的對象以json的形式儲存
參考 Getting started with JSON and jsonlite 0. 安裝包 install.packages("rjson") 1. 向量 library(rjson) vector = c(1:10) vec
2019-04-15 23:57:11
2
原创 如何解決-bash: updatedb: command not found
報錯如下: a]# updatedb -bash: updatedb: command not found 解決方法: yum install mlocate 然後使用updatedb更新 使用locate查找文件 具體操作如下
2019-04-12 23:51:15
原创 一步法中混合線性模型方程組構建和控制--blupf90
參考文獻 http://nce.ads.uga.edu/wiki/lib/exe/fetch.php?media=singlestepblupf90.pdf 1，ABLUP VS SSGBLUP 傳統ABLUP與SSGBLUP的區別
2019-04-12 23:51:15
1
原创 R語言安裝units包報錯 解決方法
解決方案 返回命令行, 我的系統是centos, 安裝 udunits2-devel, 安裝成功後, 再安裝units即可. yum install udunits2-devel 進入R語言: install.packages("
2019-04-09 23:28:30
1
原创 snakemake 學習筆記4
snakemake如何連接不同的rule 我在stackoverflow中問了一個問題, 獲得了答案, 對snakemake的理解也加深了一成. 經驗所得 每一個snakemake的rule都要有input,output, 裏面的內
2019-04-09 23:28:30
1