-
Notifications
You must be signed in to change notification settings - Fork 15
/
.Rhistory
43 lines (43 loc) · 1.01 KB
/
.Rhistory
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
rmarkdown::render('index.Rmd')
setwd(:..)
setwd("..")
rmarkdown::render('index.Rmd')
load("C:/Users/Andrew/Documents/Lieber/Research/GitHub/summerR_2015/data/charmcirc.rda")
circ[155,]
head(circ2)
dim
table(complete.cases(circ2)
)
?c
year=c(1999,1999,2000,2001)
head(year)
unique(year)
?mean
x = c(1,4,7,6,4,12,9,3)
sort(x)
x[order(x)]
order(x)
x = cut(circ$daily, c(0,10000,20000, 50000))
head(x)
table(x)
cc = factor(c("case","case","case",
"control","control","control"))
cc
relevel(cc, ref="control")
x = 1:100
cx = cut(x, breaks=c(10,25,50), labels=FALSE)
table(cx, useNA = "always")
table(cx[15:50], useNA = "always")
table(cx[15:50], useNA = "ifany")
table(cx[15:50])
table(cx)
?quantile
x
quantile(x)
cut(x, quantile(x))
cut(x, quantile(x), include.lowest = TRUE)
table(cut(x, quantile(x), include.lowest = TRUE))
table(cut(x, quantile(x,prob=c(0,33,66,100)), include.lowest = TRUE))
table(cut(x, quantile(x,prob=c(0,.33,.66,1)), include.lowest = TRUE))
source("http://bioconductor.org/biocLite.R")
biocLite("Biostrings")