[R] how to create skewness curve

quantum quantum at live.dk
Wed Jun 27 11:29:54 CEST 2012


I can make kurtosis but now skewness curve. How can I do that? I have used
the following code.



x<-(1:601-301)/50
x
dt3<-dt(x,3)
pt3<-pt(x,3)

# t distribution has mean 0; variance f/(f-2); skewness 0; excess kurtosis
6/(f-4)

dn01<-dnorm(x,0,1)
pn01<-pnorm(x,0,1)
dlogis<-dlogis(x,0,1)
plogis<-plogis(x,0,1)
# Logistic distribution is a long-tailed distribution with
# mean m and variance (pi^2 /3) s^2; skewness 0 and excess kurtosis 5/6. .
s<-sqrt(3)/pi # To get mean 0 and variance 1
dlogis<-dlogis(x,0,s)
plogis<-plogis(x,0,s)

plot(x,dlogis,col="red",type="l")
lines(x,dt3,col="blue")
lines(x,dn01)


plot(x,pn01,type="l")
lines(x,pt3)

plot(x,qnorm(pn01,0,1),type="l")
lines(x,qnorm(pt3,0,1),col="blue")
lines(x,qnorm(plogis,0,1),col="red")

--
View this message in context: http://r.789695.n4.nabble.com/how-to-create-skewness-curve-tp4634601.html
Sent from the R help mailing list archive at Nabble.com.



More information about the R-help mailing list