代码编辑器怎么下载R语言学习代码(1-6章节)
注释:参考代码来自《R语言实战 第二版》
说明:代码下载即可使用,作者使用的是R 3.6.1
#调用数据编辑器指数据框:适用于小数据集的使用
mydata<-data.frame(age=numeric(0),gender=character(0))
mydata<-edit(mydata)
fix(mydata)
getwd()
setwd("C:/Users/18226/Documents/rdatas")
N1 <- read.table("xx.txt",header = TRUE,row.names = "StudentID",sep = ",")
help(table)
#读取excel数据
install.packages("readxl")
library(readxl)
data = read_excel("xx.xlsx",sheet = 1)
help("readxl")
install.packages("openxlsx")#出现点问题,待处理
library(xlsx)
read.xlsx("ee.xlsx",sheet=1,colnames=TRUE)
#读取SPSS数据
library("foreign")#系统已经存在,可以使用read.spss()
install.packages("Hmisc")
library(Hmisc)
#use.value.labels表示让函数将带有标签的变量导入为R水平对应相同的因子
sps1 <- ("Brand.sav",use.value.labels = TRUE)
summary(sps1)
str(sps1)
#数据集的标注
mtcars
p <- data.frame(mtcars$mpg,mtcars$cyl)
names(p)[2] <- "test var labels"#变量标签
#创建值标签,leves代表变量的实际值,labels表示包含了理想值标签的字符型向量
p$`test var labels`<- factor(p$`test var labels`,levels = c(4,8),labels = c("male","female"))
#处理数据对象的实用函数
mtcars
length(mtcars)#显示对象中元素/成分的数量
dim(mtcars)#显示数据对象的维度
str(mtcars)#显示数据对象的结构
class(mtcars)#显示数据对象的类或类型
names(mtcars)#显示数据对象各成分的名称
head(mtcars)#显示数据对象的开始部分
tail(mtcars)#显示数据对象的最后部分
ls()#显示当前的对象列表
fix(mtcars)#直接编辑对象
rm(object)
rm(list = ls())#删除当前工作环境中的几乎所有对象
#图形初阶
#保存文件:png(),jpeg(),bmp(),tiff(),afile(),postcript()
pdf("mygraph.pdf")
attach(mtcars)#添加搜索路径
plot(wt,mpg)
abline(lm(mpg~wt))
title("Regression of MPG on weight")
detach(mtcars)#删除搜索路径
dev.off()
w()#打开一个新的图形窗口
()
hist(mpg)
dose<-c(20,30,40,45,60)
drugA<-c(16,20,27,40,60)
drugb<-c(15,18,25,31,40)
plot(dose,drugA,type = "b")
opar <- adonly = TRUE)#复制一份当前图形参数设置
par(lty=2,pch=17)
plot(dose,drugA,type = "b")
par(opar)
plot(dose,drugA,type = "b",lty=3,lwd=3,pch=15,cex=2)
#颜
colors()#返回所有可用的颜的名称
install.packages("RColorBrewer")
library(RColorBrewer)
n <- 7
mycolors <- brewer.pal(n,"Set1")
barplot(rep(1,n),col = mycolors)
#文本属性
plot(dose,drugA,type = "b",main = "thesfd",cex.main=3,font.main=1,cex.lab=1.5,cex.axis=1.5)
#图形尺寸和边界尺寸
dose<-c(20,30,40,45,60)
drugA<-c(16,20,27,40,60)
drugb<-c(15,18,25,31,40)
opar<- adonly = TRUE)#获取当前图形参数
par(pin=c(2,3))
par(lwd=2,cex=1.5)
par(cex.axis=0.75,font.axis=3)
plot(dose,drugA,type = "b",pch=19,lty=2,col="red")
plot(dose,drugA,type = "b",pch=23,lty=6,col="blue",bg="green")
par(opar)
help(par)
#添加文本、自定义坐标和图例
plot(dose,drugA,type = "b",col="red",lty=2,pch=2,lwd=2,
    main = "Cline Tris foe Gtd ",
    sub = "FG",xlab = "xlab",ylab = "ylab",xlim = c(0,60),ylim = c(0,70))
#标题
title()
#坐标轴
axis()
#代码清单3-2
x <-c(1:10)
y <-x
z<- 10/x#生成数据
opar <-adonly = TRUE)#获取当前图形参数
par(mar=c(5,4,4,8)+0.1)#增加边界大小
plot(x,y,type = "b",pch=21,col="red",yaxt="n",lty=3,ann = FALSE)#绘制x对y的图形
lines(x,z,type = "b",pch=22,col="blue",lty=2)
axis(2,at=x,labels = x,col.axis="red",las=2)
axis(4,at=z,labels = round(z,digits = 2),col = "blue",las=2,cex.axis=0.7,tck=-0.01)
mtext("y=1/x",side = 4,line = 3,cex.lab=1,las=2,col="blue")#添加文本
title("An Example of Creative Axes",xlab = "X Values",ylab = "Y=x")#添加标题
par(opar)
#次要参考线
library(Hmisc)
minor.tick(nx=n,ny=n,tick.ratio = n)
minor.tick(nx=2,ny=3,tick.ratio = 0.5)
#代码清单3-3
dose<-c(20,30,40,45,60)
drugA<-c(16,20,27,40,60)
drugb<-c(15,18,25,31,40)
opar<- adonly = TRUE)#获取当前图形参数
par(lwd=2,cex=1.5,font.lab=2)#增加线宽、文本大小,标签字体
plot(dose,drugA,type = "b",pch=15,lty=1,col="red",ylim = c(0,60),
    main = "Drug A VS Drug b",
    xlab = "Drug DoseG",ylab = "Dtug Response")
abline(h=c(30),lwd=1,lty=2,col="gray")
library(Hmisc)
minor.tick(nx=3,ny=3,tick.ratio = 0.5)#次要参考线
legend("topleft",inset = 0.05,title = "Drug Type",c("A","B"),
      lty = c(1,2),pch=c(15,17),col=c("red","blue"))#添加图例
#文本标注
attach(mtcars)
plot(wt,mpg,
    main = "Mile vc sd",xlab = "sdf",ylab = "sde",
    pch=18,col="red")
text(wt,mpg,
    row.names(mtcars),
    cex=0.6,pos=4,clo="blue")
detach(mtcars)
#数学标注
help("plotmath")
demo("plotmath")
#图形组合
attach(mtcars)
opar<- adonly = TRUE)#获取当前图形参数
par(mfrow=c(2,2))#组合四幅图形
plot(wt,mpg,main = "dgagg sdaf df")
plot(wt,disp,main = "dgagg sdaf df")
hist(wt,main = "dgagg sdaf df")
boxplot(wt,main = "dgagg sdaf df")
par(opar)
detach(mtcars)
attach(mtcars)
layout(matrix(c(1,1,2,3),2,2,byrow = TRUE))
plot(wt,mpg,main = "dgagg sdaf df")
plot(wt,disp,main = "dgagg sdaf df")

版权声明:本站内容均来自互联网,仅供演示用,请勿用于商业和其他非法用途。如果侵犯了您的权益请与我们联系QQ:729038198,我们将在24小时内删除。