作者:马丁摩根
日期:2019年7月22日
1 + 2
## [1] 3
x = c(1,2,3)1:3#整数的序列从1到3
## [1] 1 2 3
x + c(4,5,6)#矢量化
## [1] 5 7 9
x + 4#回收
## [1] 5 6 7
vectors.
数字()
那特点()
那逻辑()
那整数()
那复杂的()
,......NA.
: '无法使用'因素()
:受限制的“级别”集的价值观。运营
==.
那<
那<=
那>
那> =
,......|
(或者),&
(和),!!
(不是)[
,例如,X [C(2,3)]
[< -
,例如,x [c(1,3)] = x [c(1,3)]
is.na()
职能
x = rnorm(100)y = x + rnorm(100)绘图(x,y)
data.frame.
df < - data.frame(独立= x,dependent = y)头(df)
##独立依赖性## 1 -0.4338047 -0.5779168 ## 2 -0.2769985 -1.0665115 ## 3 -1.6966211 -1.8769578 ## 4 -2.648116 -1.95408116 -1.9540841 ## 5 -2.1015776 -1.1166887#6 0.7109163-0.3363154 -1.11.116688 -1.1166887#6 0.7109163 -1.1.1166876 -1.1166887#
DF [1:5,1:2]
##独立依赖性## 1 -0.4338047 -0.5779168 ## 2 -0.2769985 -1.0665115 ## 3 -1.6966211 -1.8769578 ## 4 -2.6481076 -1.954081076-0.9540841 ## 5 -2.1015776 -1.1166887 -1.1166887 -1.1166888
DF [1:5,]
##独立依赖性## 1 -0.4338047 -0.5779168 ## 2 -0.2769985 -1.0665115 ## 3 -1.6966211 -1.8769578 ## 4 -2.6481076 -1.954081076-0.9540841 ## 5 -2.1015776 -1.1166887 -1.1166887 -1.1166888
绘图(依赖〜独立,df)#'公式'界面
df [,1]
那DF [,“INDEP”]
那df [[1]]
那DF [[“INDEP”]]
那DF $ INDEP.
练习:仅绘制价值依赖> 0.
那独立> 0.
选择行
RIDX < - (DF $依赖> 0)&(DF $独立> 0)
绘制子集
绘图(依赖于〜独立,DF [RIDX,])
皮肤猫另一种方式
绘图(依赖于〜独立,df,子集=(依赖> 0)&(独立> 0))
Fit < - LM(依赖于〜独立,DF)#线性模型 - 回归ANOVA(FIT)#摘要表
##方差的分析表## ##响应:依赖于## DF SUM SQ平均SQ F值Pr(> F)##独立1 92.664 92.664 70.32 3.787E-13 *** ## Residuals 98 129.139 1.318 ## -- ## signif。代码:0'***'0.001'**'0.01'*'0.05'。'0.1''1
绘图(依赖〜独立,df)abline(适合)
lm()
:普通功能合身
: 一个目的课堂“lm”Anova()
: 一种通用的有一个特定的方法课堂“LM”班级(适合)
## [1]“lm”
方法(类=“LM”)
##Kappa ##Vcov ##查看“?方法”用于访问帮助和源代码
?“绘图”#纯旧功能或通用?“plot.formula”#方法?“plot.lm”#coptor'lm'对象的#方法,plot(fit)
库(GGPLOT2)GGPLOT(DF,AES(x =独立,y =依赖项))+ geom_point()+ geom_smooth(方法=“lm”)
图书馆(GGPLOT2)
,每次会话一次)2002年作为理解微阵列数据的平台
1,750包。专业领域:
重要的主题
资源
一个独特的特征生物体- 用于对象用于代表数据
图书馆(Biostrings)DNA < - dnastringset(C(“Aactcc”,“CTGCA”))DNA
##长度2 ##宽度SEQ ## [1] 6 AACTCC ## [2] 5 CTGCA
逆转概述(DNA)
## Dnastringset长度2 ##宽度SEQ ## [1] 6 GGAGTT ## [2] 5 TGCAG
1750'软件'包,https://biocumon.org/packages.
发现和使用,例如,deseq2.
还:
sessioninfo()
## R版本3.6.1修补(2019-07-16 R76845)##平台:X86_64-Apple-Darwin17.7.0(64位)##正在运行:Macos High Sierra 10.13.6 ## ##矩阵产品:默认## blas:/users/ma38727/bin/r-3-6-branch/lib/liblblas.dylib ## lapack:/users/ma38727/bin/r-3-6-branch/lib/librlapack.dylib### locale:## [1] en_US.UTF-8 / en_US.UTF-8 / EN_US.UTF-8 / C / EN_US.UTF-8 / EN_US.UTF-8 ## ##附加基础包:##[1] stats4并行统计图形grdevices utils datasets ## [8]方法基础## ##其他附加包:## [1] BioStrings_2.53.2 Xvector_0.25.0 iranges_2.19.10 ## [4] S4Vectors_0.23.17 Biocgenerics_0。31.5 ggplot2_3.2.0 ## [7] biocstyle_2.13.2 ## ##通过命名空间加载(且未附加):## [1] RCPP_1.0.1 Pillar_1.4.2 Compiler_3.6.1 ## [4] Biocmanager_1.30.4 zlibbioc_1。31.0 tools_3.6.1 ## [7] Digest_0.6.20评估_0.14 Tibble_2.1.3 ## [10] gtable_0.3.0 pkgconfig_2.0.2 rlang_0.4.0 ## [13] yaml_2.2.0 xfun_0.8 withr_2.1.2 ## [16]] stringr_1.4.0 dplyr_0.8.3 KNITR_1.23 ## [19] Grid_3.6.1 tidyselect_0.2.5 glue_1.3.1 ## [22] R6_2.4.0 RmarkDown_1.14 Bookdown_0.12 ## [25] purrr_0.3.2 magrittr_1.5 scales_1.0.0 ## [28] CodeTools_0.2-16 Htmltools_0.3.6Assertthat_0.2.1 ## [31] ColorSpace_1.4-1 Labeling_0.3 Stringi_1.4.3 ## [34] lazyeval_0.2.2 munsell_0.5.0 crayon_1.3.4
本教程中报告的研究得到了国家人类基因组研究所和国家癌症研究所的奖励号码U41HG004059和U24CA180996的国家癌症研究所。
该项目已收到欧洲联盟Horizo n 2020研究和创新计划下的欧洲研究委员会(ERC)的资金(Grant协议号633974)