2020-04-20

今天的内容

除说明外均引自：微信公众号生信星球
R语言基础

1.R与R Studio的安装
之前已安装并且运行无误，目前使用的是腾讯云的R Studio
 1.1 norm

> getwd()
[1] "/home/xzz"
> plot(rnorm(50))

plot 1

关于norm函数：
原文链接：https://blog.csdn.net/woodcorpse/article/details/84899814
norm是正态分布，前面加r表示生成随机正态分布的序列，其中rnorm(10)表示产生10个数；给定正太分布的均值和方差
Density(d), distribution function§, quantile function(q) and random® generation for the normal distribution with mean equal to mean and standard deviation equal to sd.
rnorm生成随机正态分布序列
pnorm可以输出正态分布的分布函数
dnorm可以输出正态分布的概率密度
qnorm给定分位数的正太分布
使用格式如下：

dnorm(x, mean = 0, sd = 1, log = FALSE)
pnorm(q, mean = 0, sd = 1, lower.tail = TRUE, log.p = FALSE)
qnorm(p, mean = 0, sd = 1, lower.tail = TRUE, log.p = FALSE)
rnorm(n, mean = 0, sd = 1)

 1.2 runif:生成-5和5之间100个随机数

> runif(100,-5,5)
  [1] -0.90212709 -4.83991612 -0.59004167  0.22565724 -4.70543667  2.47651273  1.08077521 -4.80853433 -1.10450940 -0.78815349
 [11] -0.38911068  2.15302920 -4.31591128 -3.64316848 -3.05244463  3.56707889 -1.22903412  4.50489305  1.76448619  3.69032276
 [21]  2.10145737  3.72847105  3.50437475  4.86666939  3.89960206  0.57480837 -4.10042485  4.30752521 -4.95337065 -2.28639249
 [31]  1.93132570  0.11171650 -1.61382420  4.75899099  0.56689665  2.15499275 -4.63539016  2.24101639  2.50918228 -4.28951762
 [41]  2.05307272 -0.22224267 -3.18690800  2.77842994  4.41588427 -3.92606352  2.55665017 -4.77578803  4.86153937 -4.00379830
 [51]  4.69304943  4.18273176  3.85914183  3.19773508 -1.31220351  1.46354505 -3.18789374  0.53698885  1.30677625 -0.95125769
 [61]  3.96410968 -0.84755272  0.63093391  1.08813009  4.46305885  1.79618433 -4.49605696 -1.42414226  0.17539350 -1.64502513
 [71]  3.60966034  4.88541333  3.34717601 -3.19869660 -4.99721821  0.90725231 -1.05652742  1.42876852 -1.54397673  2.45603270
 [81] -2.17845598  3.32584984  0.40630895 -1.67446191  0.28164664 -4.21120821  2.90648045  3.86621135  1.69262810  1.93347451
 [91] -4.55397531  4.50332641 -1.29783949 -1.45704195  0.37273407  4.20706636 -4.10025473 -1.62398655  0.02572065  4.68453673

 1.3 iris:鸢尾花三个亚种花萼、花瓣的数据集，以下仅展示部分

> iris
    Sepal.Length Sepal.Width Petal.Length Petal.Width    Species
1            5.1         3.5          1.4         0.2     setosa
2            4.9         3.0          1.4         0.2     setosa
3            4.7         3.2          1.3         0.2     setosa
4            4.6         3.1          1.5         0.2     setosa
5            5.0         3.6          1.4         0.2     setosa

51           7.0         3.2          4.7         1.4 versicolor
52           6.4         3.2          4.5         1.5 versicolor
53           6.9         3.1          4.9         1.5 versicolor
54           5.5         2.3          4.0         1.3 versicolor
55           6.5         2.8          4.6         1.5 versicolor

101          6.3         3.3          6.0         2.5  virginica
102          5.8         2.7          5.1         1.9  virginica
103          7.1         3.0          5.9         2.1  virginica
104          6.3         2.9          5.6         1.8  virginica

boxplot(iris$Sepal.Length~iris$Species,col = c("lightblue","lightyellow","lightpink"))

2. R基本操作
   2.1 用Rproject管理工作目录
   设置工作目录：setwd()
   查看工作目录：getwd()
   2.2 新建Rproject

点击Create Project，就打开了一个新的session，不干扰目前正在进行的项目。

 2.3 打开Rproject

专业人士Jimmy推荐只通过Rproject打开Rstudio
从此管理多个项目不成问题。
打开这个项目后，用getwd()查看一下，就会发现工作目录就是Rproject所在的目录。
用到的同一个项目的脚本/数据全都默认放在这里。
R-project管理多个R工作目录

 2.4 dir()显示文件列表
打开新建的Rproject：newbee.Rproj

> dir()
[1] "newbee.Rproj"
> list.files()
[1] "newbee.Rproj"

 2.5 基本算法

> 1+6
[1] 7
> a<-2*8
> a
[1] 16
> b=6
> c=3
> rm(a)     #移除对象a
> ls()
[1] "b" "c"
> rm(list = ls())     #移除所有对象
> ls()
character(0)

history()
ctrl+l:清空控制台