#change directory
datafile <- read.table("datanodates.csv", header = TRUE, sep = ",", na.string=".") 

SPX<-datafile[1]
NAS<-datafile[2]
INDU<-datafile[3]
T<-datafile[4]
MSFT<-datafile[5]
PTR<-datafile[6]
Y<-data.frame(SPX,NAS,INDU,T,MSFT,PTR)

##SPX
par(mfrow=c(1,1))
xname<-"S&P Daily Returns (01/2005-11/2009)"
xlab<-"S&P Returns"
hist(Y$SPX, breaks = "Sturges")
hist(Y$SPX, breaks = "FD", main = paste("Histogram of" , xname),xlab)
#hist(Y$SPX, breaks = 20)

par(mfrow=c(1,1))
plot(density(Y$SPX, kernel = "gaussian",bw="ucv"),main = paste("Density of" , xname))
lines(density(Y$SPX, kernel = "epanechnikov",bw="ucv"),col="red")


##PTR
par(mfrow=c(1,1))
xname<-"PetroChina Daily Returns (01/2005-11/2009)"
xlab<-"Petrochina Returns"
#hist(Y$PTR, breaks = "Sturges")
hist(Y$PTR, breaks = "FD", main = paste("Histogram of" , xname),xlab)
#hist(Y$PTR, breaks = 20)

par(mfrow=c(1,1))
plot(density(Y$PTR, kernel = "gaussian",bw="ucv"),main = paste("Density of" , xname))
lines(density(Y$PTR, kernel = "epanechnikov",bw="ucv"),col="red")