-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathlocal.system.R
93 lines (73 loc) · 2.9 KB
/
local.system.R
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
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
################################################################################
# LOCAL SYSTEM #
################################################################################
### Operating System Version
### MacBook Pro
ROOT_loc = "/Users/slaan3"
STORAGE_loc = "/Users/slaan3"
CLOUD_loc = "/Users/slaan3/Library/Mobile Documents/com~apple~CloudDocs/"
ONEDRIVE_loc = "/Users/slaan3/Library/CloudStorage/OneDrive-UMCUtrecht/Genomics"
# STORAGE_loc = "/Volumes/LaCie"
### Generic
LAB_loc = paste0(CLOUD_loc, "/Genomics/LabBusiness")
AEDB_loc = paste0(ONEDRIVE_loc, "/Athero-Express/AE-AAA_GS_DBs")
### Genetic and genomic data
PLINK_loc=paste0(STORAGE_loc,"/PLINK")
GWAS_loc=paste0(PLINK_loc,"/_GWAS_Datasets")
# genetic
AEGSQC_loc = paste0(PLINK_loc, "/_AE_ORIGINALS/AEGS_COMBINED_QC2018")
MICHIMP_loc=paste0(PLINK_loc,"/_AE_ORIGINALS/AEGS_COMBINED_EAGLE2_1000Gp3v5HRCr11")
# sc rna
AESCRNA_loc = paste0(PLINK_loc, "/_AE_ORIGINALS/AESCRNA/prepped_data")
# bulk rna
AERNA_loc = paste0(PLINK_loc, "/_AE_ORIGINALS/AERNA")
### Project
PROJECT_loc = paste0(STORAGE_loc, "/git/CirculatoryHealth/EndoMT_in_AE")
# use this if there is relevant information here.
ifelse(!dir.exists(file.path(PROJECT_loc, "/targets")),
dir.create(file.path(PROJECT_loc, "/targets")),
FALSE)
TARGET_loc = paste0(PROJECT_loc,"/targets")
### SOME VARIABLES WE NEED DOWN THE LINE
TRAIT_OF_INTEREST = "EndoMT_in_AE" # Phenotype
PROJECTNAME = "EndoMT_in_AE"
cat("\nCreate a new analysis directories.\n")
cat("\n- general directory\n")
ifelse(!dir.exists(file.path(PROJECT_loc)),
dir.create(file.path(PROJECT_loc)),
FALSE)
ANALYSIS_loc = paste0(PROJECT_loc)
cat("\n- for plots\n")
ifelse(!dir.exists(file.path(ANALYSIS_loc, "/PLOTS")),
dir.create(file.path(ANALYSIS_loc, "/PLOTS")),
FALSE)
PLOT_loc = paste0(ANALYSIS_loc,"/PLOTS")
ifelse(!dir.exists(file.path(PLOT_loc, "/QC")),
dir.create(file.path(PLOT_loc, "/QC")),
FALSE)
QC_loc = paste0(PLOT_loc,"/QC")
cat("\n- for output of summary results\n")
ifelse(!dir.exists(file.path(ANALYSIS_loc, "/OUTPUT")),
dir.create(file.path(ANALYSIS_loc, "/OUTPUT")),
FALSE)
OUT_loc = paste0(ANALYSIS_loc, "/OUTPUT")
# cat("\n- for baseline tables\n")
# ifelse(!dir.exists(file.path(ANALYSIS_loc, "/BASELINE")),
# dir.create(file.path(ANALYSIS_loc, "/BASELINE")),
# FALSE)
# BASELINE_loc = paste0(ANALYSIS_loc, "/BASELINE")
#
# cat("\n- for genetic analyses\n")
# ifelse(!dir.exists(file.path(ANALYSIS_loc, "/SNP")),
# dir.create(file.path(ANALYSIS_loc, "/SNP")),
# FALSE)
# SNP_loc = paste0(ANALYSIS_loc, "/SNP")
#
# cat("\n- for Cox regression results\n")
# ifelse(!dir.exists(file.path(PLOT_loc, "/COX")),
# dir.create(file.path(PLOT_loc, "/COX")),
# FALSE)
# COX_loc = paste0(PLOT_loc, "/COX")
setwd(paste0(PROJECT_loc))
getwd()
list.files()