The R Project SVN R

Rev

Rev 23351 | Go to most recent revision | Show entire file | Ignore whitespace | Details | Blame | Last modification | View Log | RSS feed

Rev 23351 Rev 25034
Line 64... Line 64...
64 
c   This version dated Aug 24, 1996.
 64 
c   This version dated Aug 24, 1996.
65 
c   Ross Ihaka, University of Auckland.
 65 
c   Ross Ihaka, University of Auckland.
66 
c   `docoef' option added Feb 17, 2003;  Martin Maechler ETH Zurich.
 66 
c   `docoef' option added Feb 17, 2003;  Martin Maechler ETH Zurich.
67 
 
 67 
 
68 
      subroutine lminfl(x, ldx, n, k, docoef, qraux, resid,
 68 
      subroutine lminfl(x, ldx, n, k, docoef, qraux, resid,
69 
     +     hat, coef, sigma)
 69 
     +     hat, coef, sigma, tol)
70 
      integer ldx, n, k, docoef
 70 
      integer ldx, n, k, docoef
71 
      double precision x(ldx,k), qraux(k), resid(n),
 71 
      double precision x(ldx,k), qraux(k), resid(n),
72 
     +     hat(n), coef(n,k), sigma(n)
 72 
     +     hat(n), coef(n,k), sigma(n)
73 
c   coef(.,.) can be dummy(1) when docoef is 0(false)
 73 
c   coef(.,.) can be dummy(1) when docoef is 0(false)
74 
 
 74 
 
Line 88... Line 88...
88 
        sigma(j) = 1.0d0
 88 
        sigma(j) = 1.0d0
89 
        call dqrsl(x, ldx, n, k, qraux, sigma, sigma, dummy,
 89 
        call dqrsl(x, ldx, n, k, qraux, sigma, sigma, dummy,
90 
     .       dummy, dummy, dummy, 10000, info)
 90 
     .       dummy, dummy, dummy, 10000, info)
91 
        do 30 i = 1, n
 91 
        do 30 i = 1, n
92 
          hat(i) = hat(i)+sigma(i)*sigma(i)
 92 
          hat(i) = hat(i)+sigma(i)*sigma(i)
- 
 
 93 
          if(hat(i) .ge. 1.0d0 - tol) hat(i) = 1.0d0
93 
   30   continue
 94 
   30   continue
94 
   40 continue
 95 
   40 continue
95 
c
 96 
c
96 
c     changes in the estimated coefficients
 97 
c     changes in the estimated coefficients
97 
c
 98 
c