Fix The Fortran Code Below:

program resampking

implicit none

integer, parameter :: ns $=10000$

real, dimension$(4)$ :: WC$1=[0\mathrm{.}8,1\mathrm{.}2,1\mathrm{.}6,5\mathrm{.}2]$

real, dimension$(5)$ :: WCex $=[3\mathrm{.}1,4\mathrm{.}9,6\mathrm{.}8,8\mathrm{.}2,15\mathrm{.}9]$

integer :: i$,$ nge$\_$random, nge$\_$bootstrap

real :: median$\_$group$1,$ median$\_$group$2,$ test$\_$statistic, p$\_$value$\_$random, p$\_$value$\_$bootstrap

real :: combined$(9),$ temp$\_$median$1,$ temp$\_$median$2,$ pseudo$\_$stat

real, dimension$($ns$)$ :: pseudo

real :: actualstat, pvalue, corr

integer :: nge, nsize

real, dimension$($:$),$ allocatable :: x$,$ y$,$ xactual

real, dimension$(3)$ :: alpha $=[0\mathrm{.}01,0\mathrm{.}05,0\mathrm{.}10]$

character$($len$=3)$ :: ans

combined $=[$WC$1,$ WCex$]!$ Correct concatenation with $//$

call calculate$\_$statistics$($WC$1,$ WCex, median$\_$group$1,$ median$\_$group$2,$ test$\_$statistic$)$

nge$\_$random $=0$

nge$\_$bootstrap $=0$

do i $=1,$ ns

call random$\_$permutation$($combined$)$

call calculate$\_$statistics$($combined$(1$:$4),$ combined$(5$:$9),$ temp$\_$median$1,$ temp$\_$median$2,$ pseudo$\_$stat$)$

if $($abs$($pseudo$\_$stat$)>=$ abs$($test$\_$statistic$))$ then

nge$\_$random $=$ nge$\_$random $+1$

end if

temp$\_$median$1=$ bootstrap$\_$resample$($WC$1)$

temp$\_$median$2=$ bootstrap$\_$resample$($WCex$)$

pseudo$\_$stat $=$ temp$\_$median$2-$ temp$\_$median$1$

if $($abs$($pseudo$\_$stat$)>=$ abs$($test$\_$statistic$))$ then

nge$\_$bootstrap $=$ nge$\_$bootstrap $+1$

end if

end do

p$\_$value$\_$random $=$ real$($nge$\_$random $+1)/$ real$($ns $+1)$

p$\_$value$\_$bootstrap $=$ real$($nge$\_$bootstrap $+1)/$ real$($ns $+1)$

nsize $=$ size$($WC$1)$

call AllocateArrays$($nsize$)$

x $=$ xactual

actualstat $=-$corr$($nsize$,$ WC$1,$ WC$1)!$ Pass WC$1$ as both x and y

nge $=0$

do i $=1,$ ns

call simpleshuffle$($nsize$,$ WC$1)$

pseudo$($i$)=-$corr$($nsize$,$ WC$1,$ WC$1)!$ Pass WC$1$ as both x and y

if $($pseudo$($i$)>=$ actualstat$)$ nge $=$ nge $+1$

end do

pvalue $=$ real$($nge $+1)/$ real$($ns $+1)$

write $(*,21)$ nge, pvalue

$21$ format $(\text{'}$ Number of resampled test statistics $>-$r $=\text{'},$ I$10,/,$ &

$\text{'}$ p$-$value $=\text{'},$ F$10\mathrm{.}4)$

write $(*,22)$ 'Level of Significance', 'Null hypothesis rejected? $($Yes or No$)\text{'}$

$22$ format $($A$30,4$X$,$ A$30)$

do i $=1,3$

if $($pvalue $<$ alpha$($i$))$ then

ans $=$ 'Yes'

else

ans $=\text{'}$No$\text{'}$

end if

write $(*,23)$ alpha$($i$),$ ans

$23$ format $($F$4\mathrm{.}2,4$X$,$ A$3)$

end do

$!$ Deallocate arrays

call DeallocateArrays

contains

subroutine calculate$\_$statistics$($group$1,$ group$2,$ median$1,$ median$2,$ test$\_$stat$)$

real, dimension$($:$),$ intent$($in$)$ :: group$1,$ group$2$

real, intent$($out$)$ :: median$1,$ median$2,$ test$\_$stat

median$1=$ median$($group$1)$

median$2=$ median$($group$2)$

test$\_$stat $=$ median$2-$ median$1$

end subroutine calculate$\_$statistics

subroutine random$\_$permutation$($arr$)$

real, dimension$($:$),$ intent$($inout$)$ :: arr

integer :: i$,$ j$,$ n

real :: temp

n $=$ size$($arr$)$

do i $=1,$ n

call random$\_$number$($temp$)$

j $=$ i $+$ int$($temp $*($n $-$ i $+1))$

temp $=$ arr$($i$)$

arr$($i$)=$ arr$($j$)$

arr$($j$)=$ temp

end do

end subroutine random$\_$permutation

subroutine simpleshuffle$($m$,$ y$)$

implicit none

integer, intent$($in$)$ :: m

real, intent$($inout$)$ :: y$($m$)$

real :: temp, r

integer :: j$,$ k

do j $=1,$ m $-1$

call random$\_$number$($r$)$

k $=$ j $+$ int$($r $*($m $-$ j $+1))$

temp $=$ y$($k$)$

y$($k$)=$ y$($j$)$

y$($j$)=$ temp

end do

end subroutine simpleshuffle

real function bootstrap$\_$resample$($data$)$

real, dimension$($:$),$ intent$($in$)$ :: data

real, dimension$($size$($data$))$ :: sample

integer :: i

real :: temp

do i $=1,$ size$($data$)$

call random$\_$number$($temp$)$

sample$($i$)=$ data$($int$($temp $*$ size$($data$))+1)$

end do

bootstrap$\_$resample $=$ median$($sample$)$

end function bootstrap$\_$resample

real function median$($arr$)$

real, dimension$($:$),$ intent$($in$)$ :: arr

real, dimension$($size$($arr$))$ :: sorted$\_$arr

integer :: n

sorted$\_$arr $=$ arr

call quicksort$($sorted$\_$arr, size$($arr$))!$ Assuming quicksort is a sorting subroutine or function

n $=$ size$($arr$)$

if $($mod$($n$,2)==0)$ then

median $=0\mathrm{.}5*($sorted$\_$arr$($n $/2)+$ sorted$\_$arr$($n $/2+1))$

else

median $=$ sorted$\_$arr$(($n $+1)/2)$

end if

end function median

subroutine AllocateArrays$($n$)$

integer, intent$($in$)$ :: n

allocate$($x$($n$),$ y$($n$),$ xactual$($n$))$

end subroutine AllocateArrays

subroutine DeallocateArrays

deallocate$($x$,$ y$,$ xactual$)$

end subroutine DeallocateArrays

end program resampking