gamma1 <- skewness.cp(10) # the skewness when the shape is 10
alpha1 <- shape.dp(gamma1) # the shape when the skewenss is gamma1; should be 10
ret <- D.gamma(10) # the derivatives of the skewness as a function of the shape
Run the code above in your browser using DataLab