Generating series of distributions with the same skewness and different kurtosis or with same kurtosis and different skewness?
And here is one way of using JohnsonFit to display how a distribution changes with changes in the moments:
library(TeachingDemos)
library(SuppDists)
tmpfun <- function(m1=0,m2=1,m3=-0.2,m4=3.7, xmin=-3, xmax=3){
parms <- JohnsonFit( c(m1,m2,m3,m4), 'use')
plot(function(x) dJohnson(x,parms), xmin, xmax)
}
tmplst <- list( Parameters=list(
m1=list('slider', init=0, from=-1, to=1, resolution=0.1),
m2=list('slider', init=1, from=0 , to=3, resolution=0.1),
m3=list('slider', init=-0.2, from=-1, to=1, resolution=0.1),
m4=list('slider', init=3.7, from=0 , to=10, resolution=0.5)),
xmin=list('numentry', init=-3, width=5),
xmax=list('numentry', init=3 , width=5) )
tkexamp( tmpfun, tmplst, plotloc='left' )
of course assuming that you have SuppDists and TeachingDemos packages installed, play around with this (and possibly change parameters and/or limits in tmplst) before demoing it as there are many combinations of the parameters that cause errors.
--
Gregory (Greg) L. Snow Ph.D.
Statistical Data Center
Intermountain Healthcare
greg.snow at imail.org
801.408.8111
-----Original Message----- From: r-help-bounces at r-project.org [mailto:r-help-bounces at r- project.org] On Behalf Of Gabor Grothendieck Sent: Wednesday, September 24, 2008 6:29 AM To: zhijie zhang Cc: R-help at stat.math.ethz.ch Subject: Re: [R] Generating series of distributions with the same skewness and different kurtosis or with same kurtosis and different skewness? Check out JohnsonFit in SuppDists. On Tue, Sep 23, 2008 at 10:59 AM, zhijie zhang <epistat at gmail.com> wrote:
Dear R users, I hope to explain the concepts of skewness and kurtosis by
generating
series of distributions with same skewness and different kurtosis or
with
same kurtosis and different skewness, but it seems that i cannot find
the
right functions. I have searched the mailing list, but no answers were found. Is it possible to do that in R? Which function could be used? Thanks a lot. -- With Kind Regards, oooO::::::::: (..)::::::::: :\.(:::Oooo:: ::\_)::(..):: :::::::)./::: ::::::(_/:::: :::::::::::::
[********************************************************************** *]
Zhi Jie,Zhang ,PHD Tel:+86-21-54237149 Dept. of Epidemiology,School of Public Health,Fudan University Address:No. 138 Yi Xue Yuan Road,Shanghai,China Postcode:200032 Email:epistat at gmail.com <Email%3Aepistat at gmail.com> Website: www.statABC.com
[********************************************************************** *]
oooO:::::::::
(..):::::::::
:\.(:::Oooo::
::\_)::(..)::
:::::::)./:::
::::::(_/::::
:::::::::::::
[[alternative HTML version deleted]]
______________________________________________ R-help at r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-
guide.html
and provide commented, minimal, self-contained, reproducible code.
______________________________________________ R-help at r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting- guide.html and provide commented, minimal, self-contained, reproducible code.