# If you donft install sgt package in your R software, firstly please run this R code .
install.packages("sgt", dependencies = TRUE)


#Condition 1-------------------------------------------------------------------
datasetD1 <- c(0.0483*1.05, 0.0483*1.025, 0.0483, 0.0483*0.975, 0.0483*0.95)
datasetD1
logD1 <- log(datasetD1)
logD1
meanD1 <- mean(logD1)
meanD1
varD1 <- var(logD1)
varD1
library(sgt)#package optimx and numDeriv are need to use sgt package 
log(0.0334)
x <- seq(-8, 0, 0.01)
plot(x, dsgt(x, mu=meanD1, sigma = sqrt((1+1/5)*varD1^2), lambda = 0, p=2, q=5/2), type="n")
curve(dsgt(x, mu=meanD1, sigma = sqrt((1+1/5)*varD1^2), lambda = 0, p=2, q=5/2), type="l", add="T")
par(new=T)
px <- c(0, 0, 0, 0, 0)
points(logD1, px)
lower_limit <- log(0.0334)
abline(v=lower_limit)
psgt(lower_limit, mu=meanD1, sigma = sqrt((1+1/5)*varD1^2), lambda = 0, p=2, q=5/2)

#Condition 2-------------------------------------------------------------------
datasetD2 <- c(0.0483*1.5, 0.0483*1.25, 0.0483, 0.0483*0.75, 0.0483*0.5)
datasetD2
logD2 <- log(datasetD2)
logD2
meanD2 <- mean(logD2)
meanD2
varD2 <- var(logD2)
varD2
library(sgt)#package optimx and numDeriv are need to use sgt package 
log(0.0334)
x <- seq(-8, 0, 0.01)
plot(x, dsgt(x, mu=meanD2, sigma = sqrt((1+1/5)*varD2^2), lambda = 0, p=2, q=5/2), type="n")
curve(dsgt(x, mu=meanD2, sigma = sqrt((1+1/5)*varD2^2), lambda = 0, p=2, q=5/2), type="l", add="T")
par(new=T)
px <- c(0, 0, 0, 0, 0)
points(logD2, px)
lower_limit <- log(0.0334)
abline(v=lower_limit)
psgt(lower_limit, mu=meanD2, sigma = sqrt((1+1/5)*varD2^2), lambda = 0, p=2, q=5/2)

#Condition 3-------------------------------------------------------------------
datasetD3 <- c(0.0483*1.64, 0.0483*1.32, 0.0483, 0.0483*0.68, 0.0483*0.36)
datasetD3
logD3 <- log(datasetD3)
logD3
meanD3 <- mean(logD3)
meanD3
varD3 <- var(logD3)
varD3
library(sgt)#package optimx and numDeriv are need to use sgt package 
log(0.0334)
x <- seq(-8, 0, 0.01)
plot(x, dsgt(x, mu=meanD3, sigma = sqrt((1+1/5)*varD3^2), lambda = 0, p=2, q=5/2), type="n")
curve(dsgt(x, mu=meanD3, sigma = sqrt((1+1/5)*varD3^2), lambda = 0, p=2, q=5/2), type="l", add="T")
par(new=T)
px <- c(0, 0, 0, 0, 0)
points(logD3, px)
lower_limit <- log(0.0334)
abline(v=lower_limit)
psgt(lower_limit, mu=meanD3, sigma = sqrt((1+1/5)*varD3^2), lambda = 0, p=2, q=5/2)

#Condition 4-------------------------------------------------------------------
datasetD4 <- c(0.0483*1.82, 0.0483*1.41, 0.0483, 0.0483*0.59, 0.0483*0.18)
datasetD4
logD4 <- log(datasetD4)
logD4
meanD4 <- mean(logD4)
meanD4
varD4 <- var(logD4)
varD4
library(sgt)#package optimx and numDeriv are need to use sgt package 
log(0.0334)
x <- seq(-8, 0, 0.01)
plot(x, dsgt(x, mu=meanD4, sigma = sqrt((1+1/5)*varD4^2), lambda = 0, p=2, q=5/2), type="n")
curve(dsgt(x, mu=meanD4, sigma = sqrt((1+1/5)*varD4^2), lambda = 0, p=2, q=5/2), type="l", add="T")
par(new=T)
px <- c(0, 0, 0, 0, 0)
points(logD4, px)
lower_limit <- log(0.0334)
abline(v=lower_limit)
psgt(lower_limit, mu=meanD4, sigma = sqrt((1+1/5)*varD4^2), lambda = 0, p=2, q=5/2)