R_L3

/*

link https://yadi.sk/d/mSNB0KgfG0kCGA
*/

plot(x)
plot(10,5)
data(Indometh)
names(Indometh)
attach(Indometh)
plot(time,conc)
means<-tapply(conc,time,mean)
means
names(means)
indo.times<-as.numeric(names(means))
plot(indo.times,means)

#Параметры xlab и ylab
plot(indo.times, means, xlab="Время", ylab="Концентрация")

#Параметры xlim и ylim
plot(indo.times,means,xlab="Время",ylab="Концентрация",ylim=c(0,15))
plot(indo.times,means,xlab="Время",ylab="Концентрация",ylim=c(0,5))

#Параметры axes и ann
plot(indo.times, means, xlab="Время",ylab="Концентрация",axes=TRUE,ann=TRUE)
plot(indo.times, means, xlab="Время",ylab="Концентрация",axes=FALSE,ann=TRUE)
plot(indo.times, means, xlab="Время",ylab="Концентрация",axes=TRUE,ann=FALSE)

#Параметр log
plot(indo.times, means, xlab="Время", ylab="Концентрация", log="x")
plot(indo.times, means, xlab="Время", ylab="Концентрация", log="y")
plot(indo.times, means, xlab="Время", ylab="Концентрация", log="xy")

#Параметр main

plot(indo.times, means, xlab="Время", ylab="Концентрация", main="Скорость выделения индометацина", type="o")
  # Пример из main
detach(Indometh)

plot(sin, -3*pi, 3*pi, col="red", lwd="3", main="График функции Sin(x)", xlab="x", ylab="Sin(x)", sub="Интервал изменения аргумента=(-3*pi, 3*pi)")
  # функция из примера
f<-function(x) {return(2*x^2-3*x+5)}
plot(f,-5,5)
  # ЗАДАНИЕ 1
f2<-function(x) {return(3*x^3-x^2+3*x-4)}
plot(f2,-5,5)

#Пример 1 R-код построение графических примитивов
X=10
Y=10
plot(x=c(0,X),y=c(0,Y),main="Графические примитивы",sub="Разноцветные прямоугольники",xlab="x",ylab="y",type="p",pch=".:")
rect(0,2,5,8,col="red")

x1=runif(1,0,X)
x1
y1=runif(1,0,Y)
y1

y2=runif(1,x1,X)
y2
y2=runif(1,y1,Y)
y2

rr=runif(1,0,1)
rg=runif(1,0,1)
rb=runif(1,0,1)
rect(x1,y1,x2,y2,col=rgb(rr,rg,rb))

xc1=4
yc1=2
xc2=8
yc2=5
r1=1
r2=2

symbols(x=c(xc1,xc2),y=c(yc1,yc2),circles=c(r1,r2),bg=c("yellow","green"),add=TRUE)
points(x=c(xc1,xc2),y=c(yc1,yc2),type="p",pch=20)
segments(xc1,yc1,xc2,yc2,lwd=2)

M=10
xx=runif(M+1,0,X)
yy=runif(M+1,0,Y)

for(i in 1:M)
{
  rr=runif(1,0,1)
  rg=runif(1,0,1)
  rb=runif(1,0,1)
  segments(xx[i],yy[i],xx[i+1],yy[i+1],lwd=3,col=rgb(rr,rg,rb))
}

#Круговые диаграммы
percent.voted<-c(60,40)
votes<-c(49.3,19.2,13.2,11.7,3.4,1.0,0.6)

names(percent.voted)<-c("Проголосоали","Не явились")
names(votes)<-c("Единая Россия","КПРФ","Справедливая Россия","ЛДПР","Яблоко","Патриоты России","Правое дело")
percent.voted
votes
par(mfrow=c(2,1))
pie(percent.voted,radius=0.9,cex=0.6,main="Явка",col=c("black","gray80"))
pie(votes,cex=0.9,radius=0.9,init.angle=-10,main="Распределение голосов",col=c(2:8))

#Точечные диаграммы
par(mfrow=c(1,1))
data(mtcars)
mtcars
dotchart(mtcars$mpg,labels=row.names(mtcars),main="Экономия топлива у 32 моделей автомобилей",xlab="Миль/галлон",cex=0.8)
x<-mtcars[order(mtcars$mpg),]
x$cyl<-factor(x$cyl)
x$color[x$cyl==4]<-1
x$color[x$cyl==6]<-2
x$color[x$cyl==8]<-3
dotchart(x$mpg,labels=row.names(x),groups=x$cyl,gcolor="blue",pch=16,color=x$color,main="Экономичность двигателя у 32 моделей автомобилей",xlab="Миль/галлон",cex=0.8)

#Гистограммы
X<-rnorm(n=50,mean=15,sd=5)
hist(X,breaks=20,freq=FALSE,col="lightblue")
plot(density(X))
plot(density(X,bw=0.8))
hist(X,breaks=20,freq=FALSE,col="lightblue",xlab="Переменная X", ylab="Плотность вероятности",main="Гистограмма, совмещенная с кривой плотности")
lines(density(X),col="red",lwd=2)

#Диаграммы размахов
boxplot(count~spray,data=InsectSprays)
boxplot(count~spray, xlab="Инсектициды",ylab="Количество выживших насекомых",main="Эффективность инсектицидов",col="coral",data=InsectSprays)
boxplot(count~spray, xlab="Инсектициды",ylab="Количество выживших насекомых",main="Эффективность инсектицидов",col="coral",horizontal=TRUE,data=InsectSprays)

#Столбиковые диаграммы

data(InsectSprays)
InsectSprays
attach(InsectSprays)
Means<-tapply(count,spray,mean)
Means

barplot(Means,col="steelblue",xlab="Инсектицид",ylab="Количество выживших насекомых")
barplot(Means,col="steelblue",xlab="Инсектицид",ylab="Количество выживших насекомых",border="red",width=sqrt(Means))
barplot(Means,density=20,col="red",horiz=T,las=1,ylab="Инсектицид",xlab="Количество выживших насекомых")

library(MASS)
data(genotype)
genotype
means=with(genotype,tapply(Wt,list(Litter,Mother),mean))
means
barplot(means,beside=TRUE,col=topo.colors(4),legend.text=rownames(means),xlab="Выводок",ylab="Вес, г",ylim=c(0,100))
barplot(means,beside=FALSE,col=topo.colors(4),xlab="Выводок",ylab="Вес, г")
sds=with(genotype,tapply(Wt,list(Litter,Mother),sd))
sds

b<-barplot(means,ylim=c(min(pretty(means-sds)),max(pretty(means+sds))),col=topo.colors(4),beside=TRUE,xpd=FALSE,ylab="Вес, г", xlab="Выводок", legend.text=rownames(means))
arrows(b, means+sds, b, means-sds,angle=90,code=3,length=0.05)

detach(InsectSprays)