Program g07befe

!     G07BEF Example Program Text

!     Mark 26.1 Release. NAG Copyright 2016.

!     .. Use Statements ..
      Use nag_library, Only: g07bef, nag_wp
!     .. Implicit None Statement ..
      Implicit None
!     .. Parameters ..
      Integer, Parameter               :: nin = 5, nout = 6
!     .. Local Scalars ..
      Real (Kind=nag_wp)               :: beta, corr, dev, gamma, sebeta,      &
                                          segam, tol
      Integer                          :: ifail, lic, maxit, n, nit
      Character (1)                    :: cens
!     .. Local Arrays ..
      Real (Kind=nag_wp), Allocatable  :: wk(:), x(:)
      Integer, Allocatable             :: ic(:)
!     .. Executable Statements ..
      Write (nout,*) 'G07BEF Example Program Results'
      Write (nout,*)

!     Skip heading in data file
      Read (nin,*)

!     Read in the problem size and control parameters
      Read (nin,*) cens, n, tol, maxit

      If (cens=='N' .Or. cens=='n') Then
        lic = 0
      Else
        lic = n
      End If
      Allocate (x(n),ic(lic),wk(n))

!     Read in data
      Read (nin,*) x(1:n)

!     Read in censor flag if required
      If (lic>0) Then
        Read (nin,*) ic(1:n)
      End If

!     Read in initial estimate of GAMMA
      Read (nin,*) gamma

!     Calculate estimates
      ifail = 0
      Call g07bef(cens,n,x,ic,beta,gamma,tol,maxit,sebeta,segam,corr,dev,nit,  &
        wk,ifail)

!     Display results
      Write (nout,99999) ' BETA  = ', beta, ' Standard error = ', sebeta
      Write (nout,99999) ' GAMMA = ', gamma, ' Standard error = ', segam

99999 Format (1X,2(A,F10.4))
    End Program g07befe