Publications by USER
Animal model
This is a page for my R publications R Pubs. You can checkout some of my few code. Animal model example problem with R Consider the mixed linear equation below: \[\begin{equation}y=Xb+Zu+e\end{equation}\] Variance components sigma_e = 40 sigma_a = 20 alpha = sigma_e/sigma_a alpha ## [1] 2 Set working directory setwd("E:/R-MME") # print worki...
1056 sym R (7586 sym/40 pcs) 3 img
R_Intro_3_Dealing with Missing Values
library(tidyverse) ## Warning: 패키지 'ggplot2'는 R 버전 4.1.2에서 작성되었습니다 ## Warning: 패키지 'tibble'는 R 버전 4.1.2에서 작성되었습니다 ## Warning: 패키지 'tidyr'는 R 버전 4.1.2에서 작성되었습니다 ## Warning: 패키지 'readr'는 R 버전 4.1.2에서 작성되었습니다 ## Warning: 패키지 ...
2557 sym R (16533 sym/112 pcs) 25 img
R_Intro_8_ggplot2 (1)
오늘 사용할 데이터는 palmerpenguin::penguins library(palmerpenguins) library(tidyverse) ## -- Attaching packages --------------------------------------- tidyverse 1.3.1 -- ## v ggplot2 3.3.5 v purrr 0.3.4 ## v tibble 3.1.6 v dplyr 1.0.7 ## v tidyr 1.2.0 v stringr 1.4.0 ## v readr 2.1.2 v forcats 0.5.1 ## -- C...
375 sym R (3054 sym/20 pcs) 8 img
R_Intro_7_variable summary (2)
이번에도 MASS 패키지 내 Cars93 데이터셋을 사용해보자. library(tidyverse) ## -- Attaching packages --------------------------------------- tidyverse 1.3.1 -- ## v ggplot2 3.3.5 v purrr 0.3.4 ## v tibble 3.1.6 v dplyr 1.0.7 ## v tidyr 1.2.0 v stringr 1.4.0 ## v readr 2.1.2 v forcats 0.5.1 ## -- Conflict...
1086 sym R (22084 sym/49 pcs)
Validation of genetic evaluations for categorical traits
Important measures in genetic evaluations Accuracy Bias dispersion Accuracy This is the correlation between the estimated breeding value \((EBV)\) and the true breeding value \((TBV)\). values near 1 indicate strong associations while values near zero indicate weak associations. Weak associations imply that the animals are poorly ranked and s...
8987 sym
Genomic selection in small populations
Genomic Selection in Small populations Julius Mugambe 2022-02-04 The accuracy of Genomic Predictions (GP) in small populations especially in developing countries and breeding companies is constrained by; terms of the structure of the reference and validation populations, response variables, genomic prediction models, validation methods. For...
3705 sym 2 img
R_Intro_10_Linear Regression Visualization
library(tidyverse) library(palmerpenguins) 1. Linear Regression without interaction terms geom_smooth()는 기본적으로 회귀선을 그려준다. 그런데 multiple regression의 경우는 어떻게 해야할까? 즉 covariate가 2개 이상일 경우에는? step 1. covariates 포함된 multi linear regression lm_fit <- lm(bill_depth_mm ~...
909 sym R (4176 sym/20 pcs) 4 img
Network Figure_220225
getwd() ## [1] "/Users/idlhy/Library/CloudStorage/OneDrive-개인/R FILE/network" rm(list=ls()) load("./network_220224.RData") library(tidygraph) library(ggraph) library(igraph) library(tidyverse) \(~\) 1990년대 linkdata_1990 <- read.csv("./network_1990_onemode.csv" , header = T #첫 행 변수명으로 ...
177 sym R (18872 sym/29 pcs) 6 img
R_Intro_9_ggplot (2)
지난번과 동일하게 오늘 사용할 데이터는 palmerpenguin::penguins library(palmerpenguins) library(tidyverse) 시작 전 데이터 훑어보기 str(penguins) %>% head(n=8) ## tibble [344 x 8] (S3: tbl_df/tbl/data.frame) ## $ species : Factor w/ 3 levels "Adelie","Chinstrap",..: 1 1 1 1 1 1 1 1 1 1 ... ## $ island ...
614 sym R (2516 sym/17 pcs) 7 img
R_Intro_5_dplyr
\(~\) 이번 예제에서는 gss가 아니라 MASS 패키지 내 Cars93 데이터를 활용 참고: (https://rfriend.tistory.com/234) library(dplyr) ## Warning: 패키지 'dplyr'는 R 버전 4.1.2에서 작성되었습니다 library(MASS) \(~\) Cars93 <- Cars93[, 1:10] #변수는 그냥 10개만 쓰겠음 str(Cars93) ## 'data.frame': 93 obs...
3558 sym R (120078 sym/123 pcs)