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
[EMAIL PROTECTED]
801.408.8111
> -----Original Message-----
> From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED]
> project.org] On Behalf Of Gabor Grothendieck
> Sent: Wednesday, September 24, 2008 6:29 AM
> To: zhijie zhang
> Cc: [EMAIL PROTECTED]
> 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 <[EMAIL PROTECTED]>
> 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:[EMAIL PROTECTED] <[EMAIL PROTECTED]>
> > Website: www.statABC.com
> >
> [**********************************************************************
> *]
> > oooO:::::::::
> > (..):::::::::
> > :\.(:::Oooo::
> > ::\_)::(..)::
> > :::::::)./:::
> > ::::::(_/::::
> > :::::::::::::
> >
> > [[alternative HTML version deleted]]
> >
> > ______________________________________________
> > R-help@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@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@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.
