13-15.1
9-6.9
35-32.9
2.1^2 * (2/15.1 + 1/32.9 + 1/6.9)
pbinom(1,6,0.02)
pbinom(1,6,0.02,lower.tail = F)
pbinom(0,6,0.02,lower.tail = F)
8000*30
5000*30
0.055*9e6
0.4^4
0.4^3*0.6*4
0.4^3*0.6
0.4^3*0.6+0.4^4
227.69/40
15+1.96*12/sqrt(48)
qt(0.975,47)
15+qt(0.975,47)*12/sqrt(48)
35*8*22
0.05/20000
2000*0.5*0.5
sqrt(2000*0.5*0.5)
qnorm(0.99)
s = sqrt(2000*0.5*0.5)
s
qnorm(0.99)*s+1000
qnorm(0.999)*s+1000
qnorm(0.999)*22+1000
0.01*0.1
0.01*0.1/0.03
x = c(45,60,70,65,55)
y = c(45,30,22,27,41)
lm(y~x)
17.67*5
17.67*5-73
mean(y)
13/3
qnorm(-2/3)
qnorm(0.67)
pnorm(0.67)
pnorm(-2/3)
pbinom(2,7,pnorm(-2/3),lower.tail = F)
pbinom(2,7,0.25,lower.tail = F)
sqrt(2.8)
-3/sqrt(2.8)
qnorm(-3/sqrt(2.8))
pnorm(-3/sqrt(2.8))
choose(12,4)
choose(6,3)
choose(6,3)*6
choose(6,4)
13/495
135/495
choose(12,4)*choose(8,4)
600/34650
135/495*2
135/495*2-600/34650
0.273*2-0.017
pbinom(9,12,0.8,lower.tail = F)
33+0.97*59
cor(X,y)
cor(x,y)
cor(x,y)*sd(y)/sd(x)
33-59*cor(x,y)*sd(y)/sd(x)
20/128e3
60e3/9e6
pnorm(2)
pnorm(2.2)
(pnorm(2.2)-pnorm(2))/pnorm(2)
(pnorm(2.2)-pnorm(2))/(1-pnorm(2))
63/428000
50e3/5e6
20/1283e
20/128e3
20/128e3/0.01
10^2.55
10^2.59
10^2.64
10^2.67
10^2.70
10^4.14
10^2.55
sqrt(2*350/4)
qnorm(0.999)
sqrt(2*350/4)*qnorm(0.999)
sqrt(2*500/4)
sqrt(2*500/4)*qnorm(0.999)
sqrt(2*500/4)*qnorm(0.999)+500
qnorm(2)
pnorm(2)
pnorm(2.2)
pnorm(2+3/15)-pnorm(2)
(pnorm(2+3/15)-pnorm(2))/(1-pnorm(2))
(pnorm(2+2/15)-pnorm(2))/(1-pnorm(2))
(pnorm(2+4/15)-pnorm(2))/(1-pnorm(2))
pnorm(qnorm(0.98),0,sqrt(1/2))
pnorm(qnorm(0.98),0,sqrt(1/2),lower.tail = F)
qnorm(.98)
pnorm(qnorm(0.98),0,sqrt(1/2),lower.tail = F)^5
qnorm(0.9)
pnorm(qnorm(0.98)-qnorm(0.9),0,sqrt(1/2),lower.tail = F)
pnorm(qnorm(0.98)-qnorm(0.9),0,sqrt(1/2),lower.tail = F)^5
pnorm(0,0,sqrt(1/2),lower.tail = F)^5
pnorm(0,0,sqrt(1/2),lower.tail = F)
pnorm(qnorm(0.98)-qnorm(0.99),0,sqrt(1/2),lower.tail = F)^5
.2*.3/0.53
10000*0.66+20000*0.23+50000*0.11
10^(-1.35)
0.01*0.11/0.03
6000^2*0.6^2
6000^2*(1-0.6^2)
200^2*100*0.7*0.3
setwd("D:/Work/Projects/Prediction/Code/embryo_selection")
pnorm(-2)
2:20
source('D:/Work/Projects/Prediction/Code/embryo_selection/EmbryoSelectionAnalysis_Lowest.R')
source('D:/Work/Projects/Prediction/Code/embryo_selection/DichotomizedTraits.R')
source('D:/Work/Projects/Prediction/Code/embryo_selection/DichotomizedTraits.R')
source('D:/Work/Projects/Prediction/Code/embryo_selection/DichotomizedTraits.R')
source('D:/Work/Projects/Prediction/Code/embryo_selection/DichotomizedTraits.R')
source('D:/Work/Projects/Prediction/Code/embryo_selection/DichotomizedTraits.R')
source('D:/Work/Projects/Prediction/Code/embryo_selection/DichotomizedTraits.R')
source('D:/Work/Projects/Prediction/Code/embryo_selection/DichotomizedTraits.R')
source('D:/Work/Projects/Prediction/Code/embryo_selection/DichotomizedTraits.R')
source('D:/Work/Projects/Prediction/Code/embryo_selection/DichotomizedTraits.R')
source('D:/Work/Projects/Prediction/Code/embryo_selection/DichotomizedTraits.R')
source('D:/Work/Projects/Prediction/Code/embryo_selection/DichotomizedTraits.R')
source('D:/Work/Projects/Prediction/Code/embryo_selection/DichotomizedTraits.R')
source('D:/Work/Projects/Prediction/Code/embryo_selection/DichotomizedTraits.R')
source('D:/Work/Projects/Prediction/Code/embryo_selection/DichotomizedTraits.R')
text(15,10,'A',cex=2)
source('D:/Work/Projects/Prediction/Code/embryo_selection/DichotomizedTraits.R')
source('D:/Work/Projects/Prediction/Code/embryo_selection/DichotomizedTraits.R')
source('D:/Work/Projects/Prediction/Code/embryo_selection/DichotomizedTraits.R')
x = c(45,60,70,65,55)
mean(x)
sd(x)
(62-50)/sqrt(100*0.25)
pnorm((62-50)/sqrt(100*0.25),lower.tail = F)
270/495 - 60/3465
1-(270/495 - 60/3465)
1-(270/495)
1-(270/495 - 60/3465)
1-(270/495 - 60/3465)-0.4671
setwd("D:/Work/Projects/Prediction/Code/embryo_selection")
pwd
pwd()
readRDS("./Data/high_risk_exclude_th_cond_parents.rds")
readRDS("./Data/high_risk_exclude_sim_cond_parents.rds")
x = c(3,1,5,4,2)
which.min(x)
source('D:/Work/Projects/Prediction/Code/embryo_selection/RunEmbryoSelectionSimulations_MixedStrategy.R')
source('D:/Work/Projects/Prediction/Code/embryo_selection/SimulateEmbryoSelection_Exclude.R')
source('D:/Work/Projects/Prediction/Code/embryo_selection/RunEmbryoSelectionSimulations_MixedStrategy.R')
source('D:/Work/Projects/Prediction/Code/embryo_selection/RunEmbryoSelectionSimulations_MixedStrategy.R')
risk_red_highrisk_sim_all
dim(risk_red_highrisk_sim_all)
source('D:/Work/Projects/Prediction/Code/embryo_selection/GenerateFigures_Mixed.R')
source('D:/Work/Projects/Prediction/Code/embryo_selection/GenerateFigures_Mixed.R')
source('D:/Work/Projects/Prediction/Code/embryo_selection/GenerateFigures_Mixed.R')
source('D:/Work/Projects/Prediction/Code/embryo_selection/GenerateFigures_Mixed.R')
source('D:/Work/Projects/Prediction/Code/embryo_selection/GenerateFigures_Mixed.R')
source('D:/Work/Projects/Prediction/Code/embryo_selection/EmbryoSelectionAnalysis_Exclude.R')
source('D:/Work/Projects/Prediction/Code/embryo_selection/GenerateFigures_Mixed.R')
source('D:/Work/Projects/Prediction/Code/embryo_selection/GenerateFigures_Mixed.R')
source('D:/Work/Projects/Prediction/Code/embryo_selection/GenerateFigures_Mixed.R')
source('D:/Work/Projects/Prediction/Code/embryo_selection/GenerateFigures_Mixed.R')
source('D:/Work/Projects/Prediction/Code/embryo_selection/GenerateFigures_Mixed.R')
source('D:/Work/Projects/Prediction/Code/embryo_selection/GenerateFigures_Mixed.R')
source('D:/Work/Projects/Prediction/Code/embryo_selection/GenerateFigures_Mixed.R')
source('D:/Work/Projects/Prediction/Code/embryo_selection/GenerateFigures_Mixed.R')
source('D:/Work/Projects/Prediction/Code/embryo_selection/GenerateFigures_Mixed.R')
pnorm(0.98,0,sqrt(0.5))
pnorm(0.98,0,0.5)
pnorm(0.98,0,sqrt(0.5),lower.tail = F)
pnorm(qnorm(0.98),0,sqrt(0.5),lower.tail = F)
pnorm(qnorm(0.98),0,sqrt(0.5),lower.tail = F)^5
pnorm(qnorm(0.98),qnorm(0.9),sqrt(0.5),lower.tail = F)^5
pnorm(qnorm(0.98),qnorm(0.9),sqrt(0.5),lower.tail = F)
0.57*28
0.47+0.227
0.697*28
sqrt(0.29)
0.46^2
0.46*0.54*2
0.5/.71
0.7*0.5*0.54
0.71*0.81 + 0.29*0.54
sqrt(0.71)
0.84*0.16*2
0.71+0.27
0.27/0.98
0.71/0.98
0.28*0.42+0.72*0.84
0.71*0.84+0.29*0.72
z
x
y
x = c(45,60,70,65,55)
y = c(45,30,22,27,41)
mean(x)
mean(y)
mean(y^2)
mean(y^2)-mean(y)
mean(y^2)-5*mean(y)
.64/(1.32*9.62)
0.05*9.62/1.32
.36*58+53
37.5/3
x1 = c(45,60,70,56,55)
mean(x)
mean(x1)
cov(x1,y)
cor(x1,y)
lm(y~x1)
cov(x1,y)/(sd(x)*sd(y))
sd(x1)
.875*9.67/9.62
33-0.879*59
.879*58-18.89
-0.77*58+78.43
sum((x-1.65)^2)/4
sqrt(sum((x-1.65)^2)/4)
sqrt(sum((y-1.65)^2)/4)
sum((y-1.65)^2)/
sum((y-1.65)^2)
sum((y-1.65)^2)/100
sqrt(sum((y-1.65)^2)/100)
sum((y-1.65)*(x-mean(x)))
3.64/(9.62*7.31)
1.65+0.04*59
-0.04*58+4.01
7.55*5-73
-0.97*58+90.39
33-ץ97*59
33-.97*59
.97*58-24.35
1.5/sqrt(2)
pnorm(0.71,lower.tail = F)
-.98*58+91.34
.98*9.67/9.62
2*456
12*456
15+53/85*100
15+53/100*85
20+50/100*85
11*306
106-54+1
106-23+1
80/125
1629+897
sqrt(1.43)
-1.07/sqrt(1.43)
pnorm(-1.07/sqrt(1.43))
qnorm(0.8)
pnorm(0.84)
9.3+2.6*0.84
9.3-2.6*0.84
9.3+2.6*qnorm(0.8)
9.3-2.6*qnorm(0.8)
2*2.6*qnorm(0.8)
11.48-7.12
9.3+2.6*qnorm(0.8)+2*2.6*qnorm(0.8)
9.3+3*2.6*qnorm(0.8)
11.48+4.36
2.4^2
2.8^2
(249*5.76+149*7.84)/(250+150-2)
0.5/sqrt(6.54*(1/250+1/150))
pnorm(0.5/sqrt(6.54*(1/250+1/150)),lower.tail = F)
pnorm(1.89,lower.tail = F)
2.4^2/150+2.8^2/250
6.54*(1/250+1/150)
2.4^2/250+2.8^2/150
sqrt(2.4^2/250+2.8^2/150)
sqrt(6.54*(1/250+1/150))
96/150
1.64*sqrt(0.24)/0.04
(1.64*sqrt(0.24)/0.04)^2
96/150*400
0.6*400
16/sqrt(400*0.4*0.6)
90+1.64*6
90+qnorm(0.95*6
90+qnorm(0.95)*6
5*0.8^4*0.2+0.8^5
pbinom(3,5,0.8,lower.tail = F)
pbinom(3,5,0.3)
1-5*0.3^4*0.7-0.3^5
pbinom(0,5,0.3)
pbinom(-1,5,0.3)
pbinom(11,5,0.3)
pbinom(1,5,0.3)
0.7^5
0.74*0.1/(0.74*0.1+0.03*0.9)
0.74*0.2/(0.74*0.2+0.03*0.8)
-2.08/2.4*2.4/15.8
qnorm(0.71)
qnorm(0.59)
qnorm(0.8)
dnorm(0)
1/sqrt(2*pi)
K = 0.1
KR = 0.2
T = qnorm(1-K)
TR = qnorm(1-KR)
aR = 0.5 # For first degree relatives
z = dnorm(T)
i = z/K
h2 = (T-TR*sqrt(1-(1-T/i)*(T^2-TR^2)))/(aR*(i+(i-T)*TR^2))
h2
z
i
T
TR
h2 = (T-TR*sqrt(1-(1-T/i)*(T^2-TR^2)))/(aR*(i+(i-T)*TR^2))
h2
K = 0.2
KR = 0.3
T = qnorm(1-K)
TR = qnorm(1-KR)
aR = 0.5 # For first degree relatives
z = dnorm(T)
i = z/K
h2 = (T-TR*sqrt(1-(1-T/i)*(T^2-TR^2)))/(aR*(i+(i-T)*TR^2))
T
TR
z
i
h2
h2^2
sqrt(h2)
h2 = (T-TR*sqrt(1-(1-T/i)*(T^2-TR^2)))/(aR*(i+(i-T)*TR^2))
h2
h2 = (T-TR*sqrt(1-(1-T/i)*(T^2+TR^2)))/(aR*(i+(i-T)*TR^2))
h2
h2 = (T-TR*sqrt(1-(1+T/i)*(T^2-TR^2)))/(aR*(i+(i-T)*TR^2))
h2
h2 = (T-TR*sqrt(1+(1-T/i)*(T^2-TR^2)))/(aR*(i+(i-T)*TR^2))
h2
h2 = (T-TR*sqrt(1-(1-i/T)*(T^2-TR^2)))/(aR*(i+(i-T)*TR^2))
h2
h2 = (T-TR*sqrt(1-(1-T/i)*(T^2-TR^2)))/(aR*(i-(i-T)*TR^2))
h2
h2 = (T-TR*sqrt(1-(1-T/i)*(T^2-TR^2)))/(aR*(i+(i-TR)*T^2))
h2
h2 = (T-TR*sqrt(1-(1-T/i)*(T^2-TR^2)))/(aR*(i+(i-T)*TR))
h2
K = 0.1
KR = 0.2
T = qnorm(1-K)
TR = qnorm(1-KR)
aR = 0.5 # For first degree relatives
z = dnorm(T)
i = z/K
h2 = (T-TR*sqrt(1-(1-T/i)*(T^2-TR^2)))/(aR*(i+(i-T)*TR))
h2
h2 = (T-TR*sqrt(1-(1-T/i)*(T^2-TR^2)))/(aR*(i+(i-T)*TR^2))
h2
h2 = (T-TR*sqrt(1-(1-T/i)*(T^2-TR^2)))/(aR*(i+(i-T^2)*TR^2))
h2
h2 = (T-TR*sqrt(1-(1-T/i)*(T^2-TR^2)))/(aR*(i+(i-T)*TR^2))
h2 = (T-TR*sqrt(1-(1-T/i)*(T^2-TR^2)))/(aR*(i+(i-T)*TR^2)); cat(sprintf("h2=%g\n"))
h2 = (T-TR*sqrt(1-(1-T/i)*(T^2-TR^2)))/(aR*(i+(i-T)*TR^2)); cat(sprintf("h2=%g\n",h2))
h2 = (TR-TR*sqrt(1-(1-T/i)*(T^2-TR^2)))/(aR*(i+(i-T)*TR^2)); cat(sprintf("h2=%g\n",h2))
h2 = (T-T*sqrt(1-(1-T/i)*(T^2-TR^2)))/(aR*(i+(i-T)*TR^2)); cat(sprintf("h2=%g\n",h2))
h2 = (T-TR*sqrt(1-(1-TR/i)*(T^2-TR^2)))/(aR*(i+(i-T)*TR^2)); cat(sprintf("h2=%g\n",h2))
h2 = (T-TR*sqrt(1-(1-T/i)*(TR^2-TR^2)))/(aR*(i+(i-T)*TR^2)); cat(sprintf("h2=%g\n",h2))
h2 = (T-TR*sqrt(1-(1-T/i)*(T^2-TR^2)))/(aR*(i+(i-T)*T^2)); cat(sprintf("h2=%g\n",h2))
h2 = (T-TR*sqrt(1-(1-T/i)*(T^2-TR^2)))/(aR*(i+(i-TR)*TR^2)); cat(sprintf("h2=%g\n",h2))
K = 0.1
KR = 0.3
T = qnorm(1-K)
TR = qnorm(1-KR)
aR = 0.5 # For first degree relatives
z = dnorm(T)
i = z/K
h2 = (T-TR*sqrt(1-(1-T/i)*(T^2-TR^2)))/(aR*(i+(i-TR)*TR^2)); cat(sprintf("h2=%g\n",h2))
K = 0.1
KR = 0.3
T = qnorm(1-K)
TR = qnorm(1-KR)
aR = 0.5 # For first degree relatives
z = dnorm(T)
i = z/K
h2 = (T-TR*sqrt(1-(1-T/i)*(T^2-TR^2)))/(aR*(i+(i-T)*TR^2)); cat(sprintf("h2=%g\n",h2))
K = 0.1
KR = 0.2
T = qnorm(1-K)
TR = qnorm(1-KR)
aR = 0.5
z = dnorm(T)
i = z/K
h2 = (T-TR*sqrt(1-(1-T/i)*(T^2-TR^2)))/(aR*(i+(i-T)*TR^2)); cat(sprintf("h2=%g\n",h2))
T
TR
z
i
-0.13*25+9.36
pnorm(-1.89)
r2A = 0.1
r2B = 0.2
rA = sqrt(r2A)
rB = sqrt(r2B)
mu = c(0,0)
sigma = 1/2 * matrix(c(r2A,rho*rA*rB,rho*rA*rB,r2B),nrow=2)
mu
rho - 0.01
rho = 0.01
mu = c(0,0)
sigma = 1/2 * matrix(c(r2A,rho*rA*rB,rho*rA*rB,r2B),nrow=2)
sigma
sigma*2
rho*rA*rB
rA
rB
?mvrnorm
mvrnorm(n = 1,mu,sigma)
mvrnorm(n = 1,mu,sigma)
mvrnorm(n = 1,mu,sigma)
nfam = 5
mvrnorm(n = nfam,mu,sigma)
source('D:/Work/Projects/Prediction/Code/embryo_selection/SimulateEmbryoSelection_TwoTraits.R')
simulate_lowest_risk_two_traits(0.1,0.2,0.05,0.02,0.06,1:20)
source('D:/Work/Projects/Prediction/Code/embryo_selection/SimulateEmbryoSelection_TwoTraits.R')
source('D:/Work/Projects/Prediction/Code/embryo_selection/SimulateEmbryoSelection_TwoTraits.R')
results
source('D:/Work/Projects/Prediction/Code/embryo_selection/SimulateEmbryoSelection_TwoTraits.R')
res
source('D:/Work/Projects/Prediction/Code/embryo_selection/SimulateEmbryoSelection_TwoTraits.R')
res
source('D:/Work/Projects/Prediction/Code/embryo_selection/SimulateEmbryoSelection_TwoTraits.R')
res
res = simulate_lowest_risk_two_traits(0.1,0.1,-0.3,0.02,0.06,1:20)
res
res = simulate_lowest_risk_two_traits(0.1,0.1,-0.3,0.02,0.02,1:20)
res
rhos = c(-0.3,seq(-0.2,0,by=0.05))
rhose
rhos
source('D:/Work/Projects/Prediction/Code/embryo_selection/SimulateEmbryoSelection_TwoTraits.R')
source('D:/Work/Projects/Prediction/Code/embryo_selection/SimulateEmbryoSelection_TwoTraits.R')
source('D:/Work/Projects/Prediction/Code/embryo_selection/SimulateEmbryoSelection_TwoTraits.R')
risk_red_lowest_sim_all
risk_red_lowest_sim_all[1,,]
simulate_lowest_risk_two_traits(0.1,0.1,0,0.05,0.05,1:20,1e4)
risk_red_lowest_sim_all[1,,]
source('D:/Work/Projects/Prediction/Code/embryo_selection/SimulateEmbryoSelection_TwoTraits.R')
risk_red_lowest_sim_all[1,,]
install.packages("viridis")
library(viridis)
?viridis
?magma
magma(10)
saveRDS(risk_red_lowest_sim_all, file = "./Data/lowest_risk_two_taits.rds")
as.expression(bquote('rho'*'='*.(rhos[rhi])))
source('D:/Work/Projects/Prediction/Code/embryo_selection/GenerateFigures_TwoTraits.R')
source('D:/Work/Projects/Prediction/Code/embryo_selection/GenerateFigures_TwoTraits.R')
source('D:/Work/Projects/Prediction/Code/embryo_selection/GenerateFigures_TwoTraits.R')
source('D:/Work/Projects/Prediction/Code/embryo_selection/GenerateFigures_TwoTraits.R')
source('D:/Work/Projects/Prediction/Code/embryo_selection/GenerateFigures_TwoTraits.R')
risk_red_lowest_sim_all[1,,]
source('D:/Work/Projects/Prediction/Code/embryo_selection/GenerateFigures_TwoTraits.R')
source('D:/Work/Projects/Prediction/Code/embryo_selection/GenerateFigures_TwoTraits.R')
source('D:/Work/Projects/Prediction/Code/embryo_selection/GenerateFigures_TwoTraits.R')
source('D:/Work/Projects/Prediction/Code/embryo_selection/GenerateFigures_TwoTraits.R')
source('D:/Work/Projects/Prediction/Code/embryo_selection/GenerateFigures_TwoTraits.R')
source('D:/Work/Projects/Prediction/Code/embryo_selection/GenerateFigures_TwoTraits.R')
legs
source('D:/Work/Projects/Prediction/Code/embryo_selection/GenerateFigures_TwoTraits.R')
source('D:/Work/Projects/Prediction/Code/embryo_selection/GenerateFigures_TwoTraits.R')
source('D:/Work/Projects/Prediction/Code/embryo_selection/GenerateFigures_TwoTraits.R')
source('D:/Work/Projects/Prediction/Code/embryo_selection/GenerateFigures_TwoTraits.R')
source('D:/Work/Projects/Prediction/Code/embryo_selection/GenerateFigures_TwoTraits.R')
source('D:/Work/Projects/Prediction/Code/embryo_selection/GenerateFigures_TwoTraits.R')
source('D:/Work/Projects/Prediction/Code/embryo_selection/GenerateFigures_TwoTraits.R')
source('D:/Work/Projects/Prediction/Code/embryo_selection/GenerateFigures_TwoTraits.R')
source('D:/Work/Projects/Prediction/Code/embryo_selection/GenerateFigures_TwoTraits.R')
source('D:/Work/Projects/Prediction/Code/embryo_selection/GenerateFigures_TwoTraits.R')
source('D:/Work/Projects/Prediction/Code/embryo_selection/GenerateFigures_TwoTraits.R')
source('D:/Work/Projects/Prediction/Code/embryo_selection/GenerateFigures_TwoTraits.R')
source('D:/Work/Projects/Prediction/Code/embryo_selection/GenerateFigures_TwoTraits.R')
source('D:/Work/Projects/Prediction/Code/embryo_selection/GenerateFigures_TwoTraits.R')
source('D:/Work/Projects/Prediction/Code/embryo_selection/GenerateFigures_TwoTraits.R')
source('D:/Work/Projects/Prediction/Code/embryo_selection/GenerateFigures_TwoTraits.R')
source('D:/Work/Projects/Prediction/Code/embryo_selection/GenerateFigures_TwoTraits.R')
source('D:/Work/Projects/Prediction/Code/embryo_selection/GenerateFigures_TwoTraits.R')
source('D:/Work/Projects/Prediction/Code/embryo_selection/GenerateFigures_TwoTraits.R')
source('D:/Work/Projects/Prediction/Code/embryo_selection/GenerateFigures_TwoTraits.R')
source('D:/Work/Projects/Prediction/Code/embryo_selection/GenerateFigures_TwoTraits.R')
source('D:/Work/Projects/Prediction/Code/embryo_selection/GenerateFigures_TwoTraits.R')
source('D:/Work/Projects/Prediction/Code/embryo_selection/GenerateFigures_TwoTraits.R')
?plotmath
source('D:/Work/Projects/Prediction/Code/embryo_selection/GenerateFigures_TwoTraits.R')
source('D:/Work/Projects/Prediction/Code/embryo_selection/GenerateFigures_TwoTraits.R')
source('D:/Work/Projects/Prediction/Code/embryo_selection/SimulateEmbryoSelection_TwoTraits.R')
source('D:/Work/Projects/Prediction/Code/embryo_selection/SimulateEmbryoSelection_TwoTraits.R')
source('D:/Work/Projects/Prediction/Code/embryo_selection/GenerateFigures_TwoTraits.R')
?rev
source('D:/Work/Projects/Prediction/Code/embryo_selection/GenerateFigures_TwoTraits.R')
source('D:/Work/Projects/Prediction/Code/embryo_selection/GenerateFigures_TwoTraits.R')
source('D:/Work/Projects/Prediction/Code/embryo_selection/GenerateFigures_TwoTraits.R')
source('D:/Work/Projects/Prediction/Code/embryo_selection/GenerateFigures_TwoTraits.R')
source('D:/Work/Projects/Prediction/Code/embryo_selection/GenerateFigures_TwoTraits.R')
source('D:/Work/Projects/Prediction/Code/embryo_selection/GenerateFigures_TwoTraits.R')
dim(risk_red_lowest_sim_all)
risk_red_lowest_sim_all[,5,5]
source('D:/Work/Projects/Prediction/Code/embryo_selection/GenerateFigures_TwoTraits.R')
source('D:/Work/Projects/Prediction/Code/embryo_selection/GenerateFigures_TwoTraits.R')
source('D:/Work/Projects/Prediction/Code/embryo_selection/GenerateFigures_TwoTraits.R')
source('D:/Work/Projects/Prediction/Code/embryo_selection/GenerateFigures_TwoTraits.R')
risk_red_lowest_sim_all[,5,2]