LinearRegression/phenotype_generation.sh at master · Electric-Jam/LinearRegression · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
#For type1 error

basedir=$1
outbase=$2
regpath=$3

mkdir -p $outbase/CSE284_phenotype_1/100_samples
mkdir -p $outbase/CSE284_phenotype_1/1k_samples
mkdir -p $outbase/CSE284_phenotype_1/10k_samples

Rscript $regpath/phenotype_generation_1.R \
--fam $basedir/simulate.100.fam \
--n_phenotypes 100 \
--out $outbase/CSE284_phenotype_1/100_samples


Rscript $regpath/phenotype_generation_1.R \
--fam $basedir/simulate.1000.fam \
--n_phenotypes 100 \
--out $outbase/CSE284_phenotype_1/1k_samples


Rscript $regpath/phenotype_generation_1.R \
--fam $basedir/simulate.10000.fam \
--n_phenotypes 100 \
--out $outbase/CSE284_phenotype_1/10k_samples


#For type2 error

mkdir -p $outbase/CSE284_phenotype_2/100_samples
mkdir -p $outbase/CSE284_phenotype_2/1k_samples
mkdir -p $outbase/CSE284_phenotype_2/10k_samples

Rscript $regpath/phenotype_generation_2.R \
--bfile $basedir/simulate.100 \
--n_phenotypes 100 \
--n_negative_causal_variants 100 \
--n_positive_causal_variants 100 \
--negative_beta_min -1 \
--negative_beta_max -0.5 \
--positive_beta_min 1 \
--positive_beta_max 2 \
--out $outbase/CSE284_phenotype_2/100_samples

Rscript $regpath/phenotype_generation_2.R \
--bfile $basedir/simulate.1000 \
--n_phenotypes 100 \
--n_negative_causal_variants 100 \
--n_positive_causal_variants 100 \
--negative_beta_min -1 \
--negative_beta_max -0.5 \
--positive_beta_min 1 \
--positive_beta_max 2 \
--out $outbase/CSE284_phenotype_2/1k_samples

Rscript $regpath/phenotype_generation_2.R \
--bfile $basedir/simulate.10000 \
--n_phenotypes 100 \
--n_negative_causal_variants 100 \
--n_positive_causal_variants 100 \
--negative_beta_min -1 \
--negative_beta_max -0.5 \
--positive_beta_min 1 \
--positive_beta_max 2 \
--out $outbase/CSE284_phenotype_2/10k_samples