-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathstepsRUriWatch.r
62 lines (48 loc) · 1.15 KB
/
stepsRUriWatch.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
remove(mirs)
install.packages("rlang")
install.packages("readxl")
library(readxl)
df <- read_excel("C:\\Users\\dsa\\Desktop\\UriWatchDataSet001.xlsx", sheet=1 , col_types = "numeric")
p24=df[4:613,11]
p25=df[4:613,20]
p26=df[4:613,29]
c7=df[4:613,38]
c8=df[4:613,47]
c9=df[4:613,56]
df <- read_excel("C:\\Users\\dsa\\Desktop\\UriWatchDataSet001.xlsx", sheet=1 )
mirs <- matrix(0,610, 7)
mirs[1:610,1]=t(mirList)
mirs[1:610,2]= as.numeric(
t(c7))
mirs[1:610,3]= as.numeric(
t(c8))
mirs[1:610,4]=t(c9)
mirs[1:610,5]=t(p24)
mirs[1:610,6]=t(p25)
typeof(mirs)
#[1] "double"
dim(mirs)
#[1] 610 6
install.packages("FactoMineR")
library(FactoMineR)
res <- PCA(mirs)
df <- read_excel("C:\\Users\\dsa\\Desktop\\UriWatchDataSet001.xlsx", sheet=1)
mirList=df[4:613,1]
mirList
# A tibble: 610 x 1
X__1
#<chr>
# 1 hsa-miR-10b-5p
#2 hsa-miR-10a-5p
##3 hsa-miR-486-5p
#4 hsa-let-7f-5p
#5 hsa-let-7a-5p
#6 hsa-miR-99b-5p
#7 hsa-miR-191-5p
#8 hsa-miR-100-5p
#9 hsa-miR-30a-5p
#10 hsa-let-7b-5p
# ... with 600 more rows
mirs[1:610,2:7]=as.numeric(mirs[1:610,2:7])
PCA(df[,2:7])
PCA(as.numeric(df3[,1:610]))