all,
I'm trying to simulate data for different sample size (let say n= 10,20, and 30). At the same time, i need the parameter (in my case are skewness and kurtosis) for the data for each n is the same. How can i do it in R.
I already tried do like this:
set.seed(100)
n<-10
true.skew<-0
true.kur<-3
library(PearsonDS)
e<-rpearson(n,moments=c(mean=0,variance=1, skewness=true.skew, kurtosis=true.kur))
library(moments)
skewness(e)
kurtosis(e)
and the result are:
>skewness(e)
[1] 0.2110712
> kurtosis(e)
[1] 1.950325
>
then i increase the sample size to n = 20,
set.seed(100)
n1<-20
true.skew<-0
true.kur<-3
library(PearsonDS)
e1<-rpearson(n1,moments=c(mean=0,variance=1, skewness=true.skew, kurtosis=true.kur))
library(moments)
skewness(e1)
kurtosis(e1)
output:
skewness(e1)
[1] 1.311147
> kurtosis(e1)
[1] 5.517818
the problems is, i got the different value of skewness and kurtosis.
Anybody can help me to do this?
Thank you.