Program g07cafe
! G07CAF Example Program Text
! Mark 24 Release. NAG Copyright 2012.
! .. Use Statements ..
Use nag_library, Only: g07caf, nag_wp
! .. Implicit None Statement ..
Implicit None
! .. Parameters ..
Integer, Parameter :: nin = 5, nout = 6
! .. Local Scalars ..
Real (Kind=nag_wp) :: clevel, df, dl, du, prob, t, xmean, &
xstd, ymean, ystd
Integer :: ifail, nx, ny
Character (1) :: equal, tail
! .. Executable Statements ..
Write (nout,*) 'G07CAF Example Program Results'
Write (nout,*)
! Skip heading in data file
Read (nin,*)
! Read in the sample sizes, means and standard deviations
Read (nin,*) nx, xmean, xstd
Read (nin,*) ny, ymean, ystd
! Display data
Write (nout,*) 'Sample X'
Write (nout,99996) ' Sample size = ', nx
Write (nout,99995) ' Mean =', xmean
Write (nout,99995) ' Standard deviation =', xstd
Write (nout,*)
Write (nout,*) 'Sample Y'
Write (nout,99996) ' Sample size = ', ny
Write (nout,99995) ' Mean =', ymean
Write (nout,99995) ' Standard deviation =', ystd
d_lp: Do
! Read in the type of statistic and CI required
Read (nin,*,Iostat=ifail) clevel, tail, equal
If (ifail/=0) Then
Exit d_lp
End If
! Calculate statistic
ifail = 0
Call g07caf(tail,equal,nx,ny,xmean,ymean,xstd,ystd,clevel,t,df,prob, &
dl,du,ifail)
! Display results
Write (nout,*)
If (equal=='E' .Or. equal=='e') Then
Write (nout,*) 'Assuming population variances are equal.'
Else
Write (nout,*) 'No assumptions about population variances.'
End If
Write (nout,*)
Write (nout,99999) 't test statistic = ', t
Write (nout,99998) 'Degrees of freedom = ', df
Write (nout,99997) 'Significance level = ', prob
Write (nout,*)
Write (nout,99999) 'Difference in means'
Write (nout,99999) ' Value = ', xmean - ymean
Write (nout,99999) ' Lower confidence limit = ', dl
Write (nout,99999) ' Upper confidence limit = ', du
Write (nout,99999) ' Confidence level = ', clevel
End Do d_lp
99999 Format (1X,A,F10.4)
99998 Format (1X,A,F8.1)
99997 Format (1X,A,F8.4)
99996 Format (1X,A,I5)
99995 Format (1X,A,E11.4)
End Program g07cafe