attach(mtcars)
plot(wt,mpg)
abline(lm(mpg~wt))
title("Regression of MPG on Weight")
detach(mtcars)
函数colors()可以返回所有可用颜色的名称
R中也有多种用于创建连续型颜色向量的函数，包括rainbow()、heat.colors()、terrain.colors()、topo.colors()以及cm.colors()。举例来说，rainbow(10)可以生成10种连续的“彩虹型”颜色。多阶灰度色可使用gray()函数生成。这时要通过一个元素值为0和1之间的向量来指定各颜色的灰度。gray(0:10/10)将生成10阶灰度色。
n<-10
mycolors<-rainbow(n)
pie(rep(1,n),labels=mycolors,col=mycolors)
mygrays<-gray(0:n/n)
pie(rep(1,n),labels=mygrays,col=mygrays)
title("rainbow")
pie(rep(1,n),labels=mygrays,col=mygrays)
文本大小的参数【cex/cex.axis:缩放倍数；cex.lab坐标轴标签的缩放倍数；cex.main/cex.sub标题和副标题的缩放倍数
par(font.lab=3,cex.lab=1.5,font.main=4,cex.main=2)
图形外观实例1
dose<-c(20,30,40,45,60)
drugA<-c(16,20,27,40,60)
drugB<-c(15,18,25,31,40)
opar<-par(no.readonly=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,drugB,type="b",pch=23,lty=6,col="blue",bg="green")
par(opar)
图形外观实例2
plot(dose,drugA,type="b",col="red",lty=2,pch=2,lwd=2,
main="Clinical Trials for DrugA",
sub="This is hypothetical data",
xlab="Dosage",ylab="Drug Response",
xlim=c(0,60),ylim=c(0,70))
图形外观实例3
plot(dose,drugA,type="b",col="red",lty=2,pch=2,lwd=2,xlim=c(0,60),ylim=c(0,70))
title(main="main title",sub="sub-title",xlab="x-axis label",ylab="y-axis label")
数学标注 查阅help(plotmath)
图形的组合1
attach(mtcars)
opar<-par(no.readonly=TRUE)
par(mfrow=c(2,2))
plot(wt,mpg,main="Scatterplot of wt vs.mpg")
plot(wt,disp,main="Scatterplot of wt vs disp")
hist(wt,main="Histogram of wt")
boxplot(wt,main="Boxplot of wt")
par(opar)
detach(mtcars)
图形的组合2
attach(mtcars)
opar<-par(no.readonly=TRUE)
par(mfrow=c(1,3))
boxplot(wt)
boxplot(mpg)
boxplot(disp)
par(opar)
detach(mtcars)
图形的组合3
attach(mtcars)
layout(matrix(c(1,1,2,3),2,2,byrow=TRUE))
boxplot(wt)
boxplot(mpg)
boxplot(disp)
detach(mtcars)
图形的组合4
attach(mtcars)
layout(matrix(c(1,1,2,3),2,2,
byrow=TRUE),width=c(3,1),heights=c(1,2))
boxplot(wt)
boxplot(mpg)
boxplot(disp)
detach(mtcars)
创建leadership数据框
manager<-c(1,2,3,4,5)
data<-c("10/24/08","10/28/08","10/1/08","10/12/08","5/1/09")
country<-c("US","US","UK","UK","UK")
gender<-c("M","F","F","M","F")
age<-c(32,45,25,39,99)
q1<-c(5,3,3,3,2)
q2<-c(4,5,5,3,2)
q3<-c(5,2,5,4,1)
q4<-c(5,5,5,NA,2)
q5<-c(5,5,2,NA,1)
leadership<-data.frame(manager,data,country,gender,age,q1,q2,q3,q4,q5,stringsAsFactors=FALSE)
变量的重命名
1. fix(leadership)
leadership
缺失值
is.na(leadership)
在分析中排除缺失值,na.rm=TRUE选项，可以在计算之前移除缺失值并使用剩余值进行计算
函数na.omit()移除所有含有缺失值的观测
newdata<-na.omit(leadership)
getwd()
apply(x,2,mean)
apply函数的使用格式：
apply(x,MARGIN,FUN) x为数据对象，MARGIN=1表示行，MARGIN=2表示列
mydata<-matrix(rnorm(30),nrow=6)
apply(mydata,1,mean)
apply(mydata,2,mean)
函数t()即可对一个矩阵或数据框进行转置
整合数据
aggregate(x,by,FUN) x待折叠的数据对象，by是一个变量名组成的列表，FUN函数
attach(mtcars)
aggdata<-aggregate(mtcars,by=list(cyl,gear),FUN=mean,na.rm=TRUE)
attach(x)
y<-aggregate(x,by=list(Length,Number),FUN=mean,na.rm=TRUE)