R-test

#iris data
print ("this is first project script test")
str (iris);

#body fat data
library("mboost");
data("bodyfat", package = "TH.data");
str (bodyfat)

#save object in file
a = 1:10;
print (a);
save (a, file="data/a.Rdata")
rm (a);
load("data/a.Rdata")
print (a)

#read and save data in csv
var1 = 1:5;
var2 = (1:5)/10
var3 = c("hello", "world", "china", "usa", "uk")
df1 = data.frame(var1, var2, var3)
names(df1) = c("column1", "column2", "column3");
write.csv (df1, "data/b.csv", row.names = FALSE)
df2 = read.csv("data/b.csv");
print (df2)

#cat datas
print (dim (iris));
print (names(iris))
print (str(iris))
print (attributes(iris))
iris[1:5,]
head (iris)
tail (iris)
summary (iris)

#var
var (iris$Sepal.Length)
#hist zhifangtu
hist (iris$Sepal.Length)
#midu
plot (density(iris$Sepal.Length))
#table
table (iris$Species)
#pie
pie (table(iris$Species))
#bar
barplot (table(iris$Species))

#xiefangcha
cov (iris$Sepal.Length, iris$Petal.Length)
#iris 
iris[, 1:4]
#xiefnagcha matrix
cov (iris[, 1:4])
#xiefangcha xishu
cor (iris$Sepal.Length, iris$Petal.Length)
#boxplot
boxplot (iris$Sepal.Length~Species, data=iris)
#pairs
pairs (iris)
#plot jitty
plot (jitter(iris$Sepal.Length), jitter(iris$Sepal.Width))

#3d
library ("scatterplot3d")
scatterplot3d (iris$Sepal.Length, iris$Sepal.Width, iris$Petal.Length)