library(BSDA)
## 载入需要的程辑包:lattice
## 
## 载入程辑包:'BSDA'
## The following object is masked from 'package:datasets':
## 
##     Orange
#Engagement
zvalue <- zsum.test(mean.x=.48, sigma.x = .117,  n.x = 3,
                    mean.y = .486,  sigma.y = .134,  n.y = 11,
                    alternative = "two.sided",
                    mu = 0,
                    conf.level = 0.95
)
zvalue
## 
##  Two-sample z-Test
## 
## data:  Summarized x and y
## z = -0.076229, p-value = 0.9392
## alternative hypothesis: true difference in means is not equal to 0
## 95 percent confidence interval:
##  -0.16027  0.14827
## sample estimates:
## mean of x mean of y 
##     0.480     0.486

#Transformatiol leadership

zvalue <- zsum.test(mean.x=.423, sigma.x = .222,  n.x = 5,
                    mean.y = .469,  sigma.y = .163,  n.y = 18,
                    alternative = "two.sided",
                    mu = 0,
                    conf.level = 0.95
)
zvalue
## 
##  Two-sample z-Test
## 
## data:  Summarized x and y
## z = -0.4321, p-value = 0.6657
## alternative hypothesis: true difference in means is not equal to 0
## 95 percent confidence interval:
##  -0.2546497  0.1626497
## sample estimates:
## mean of x mean of y 
##     0.423     0.469

#Proactive behavior

zvalue <- zsum.test(mean.x=.351, sigma.x = .110,  n.x = 8,
                    mean.y = .412,  sigma.y = .190,  n.y = 26,
                    alternative = "two.sided",
                    mu = 0,
                    conf.level = 0.95
)
zvalue
## 
##  Two-sample z-Test
## 
## data:  Summarized x and y
## z = -1.1326, p-value = 0.2574
## alternative hypothesis: true difference in means is not equal to 0
## 95 percent confidence interval:
##  -0.16656479  0.04456479
## sample estimates:
## mean of x mean of y 
##     0.351     0.412