# #——风格,回声= FALSE,结果= '飞机 '-------------------------------------------------------- BiocStyle:减价()选项(宽度= 100,max.print = 1000)选项(useFancyQuotes = FALSE) knitr:: opts_chunk设置(eval = as.logical (Sys美元。采用“KNITR_EVAL”,“真正的”)),缓存= as.logical (Sys。采用“KNITR_CACHE”,“真正的”)),错误= FALSE ) ## ------------------------------------------------------------------------------------------------- 意思是(1:10)rnorm(1:10)总结(rnorm (1:10 )) ## ------------------------------------------------------------------------------------------------- 数据(iris) #找到那些行,花瓣。iris$花瓣. width ==0.2，] #找到那些行萼片。长度小于4.5虹膜[虹膜$萼片。长度< 4.5,)#发现所有行属于setosa setosa_iris =虹膜(虹膜$物种= =“setosa]暗(setosa_iris)头(setosa_iris ) ## ------------------------------------------------------------------------------------------------- # 下降的列包含字符即物种虹膜< -虹膜[,!#简单:colMeans(iris) #简化结果sapply(iris, mean) #找出每一行的平均值。Apply (iris, 1, mean) #更简单:rowMeans(虹膜 ) ## ------------------------------------------------------------------------------------------------- # 定义一个向量x < - rnorm(1000) #矢量化计算y < - x + rnorm(1000年,sd = 8) #对象构造df < - data.frame (x = x, y = y) #线性模型符合< - lm (y ~ x,df ) ## ------------------------------------------------------------------------------------------------- 票面价值(mfrow = c(1、2)情节(y ~ x, df, cex.lab = 2) abline(健康,坳=“红色”,lwd = 2)图书馆(ggplot2) ggplot (df, aes (x,y)) + geom_point () + stat_smooth(=“lm方法 ") ## ---- sessionInfo ---------------------------------------------------------------------------------- sessionInfo ()