Using OWN-structure
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

Using OWN-structure



Hi ASREMLers!
I am trying to specify own variance structures using facility OWN. Using the
'clasical' definition of the variance structure, I get the residual
component AR=0.231273.
File 'varstr.as' was:

********************
Using own

 ident !P                       # Animal
 x                              # Testday
 anim 2476 !I                   # Animal to permanent
 mfi                            # Mean daily feed intake
 iniwt                          # Initial weight-mean iniwt 46
 birthfa 5 !I                   # Birth family
 season 21 !I                   # Season
 x2                             # Testday*testday
 apparat 26 !I                  # Apparat

/home/suatt/suatt141/random4/woche.ped
/home/suatt/suatt141/random4/woche.dat

mfi ~ mu iniwt birthfa season apparat x x2,
      in.bi in.se in.app bi.se bi.app se.app

1 2 0

anim                    # R=IxRo
5 0 OWN1 0.01 !TR
********************

and generated output file 'varstr.asr'

*******************
   5  OWNxx       0.01
   12380 records assumed sorted    5 within 2476
   Forming  889  equations:  343  dense
   Initial updates will be shrunk by factor    0.316
   Fault     4  READING OWN structures varstr
*******************

I am working under Linux. I have tried with original MYOWNGDG.f90 in
directory $HOME/asreml/ (together with asreml) and in directory
$HOME/random/ (together with the asreml-program). File 'varstrB.own' was
also generated.

_________________________________________________________
    5    1    1
0.1000000D-01

 This file was written by ASReml for reading by your
 program MYOWNGDG
 ASReml writes this file, runs your program and then reads
 varstrB.own
 which it presumes has the following format:
 The first lines should agree with the top of this file
 specifying the order of the matrices         (    5)
            the number of variance parameters (    1)
      and a control  parameter you can specify (    1).
 These are written in (3I5 ) format. They are followed by
 the list of variance parameters  written in (6D13.7) format.
 Follow this with     2 matrices  written in (6D13.7) format.
 These are to be each of    15 elements being lower triangle
 row-wise of the G matrix and  its derivatives with respect to
 the parameters in turn.

___________________________________________________________

Any suggestions, perhaps changes of the f90-program?
Thanks in advance and have a nice weekend.

Justo Lorenzo Bermejo


 			                                                                     
  Justo Lorenzo Bermejo                                    
  Institut fr Tierzucht und Tierhaltung                   
  Olshausenstr. 40                                           
  24118 Kiel  	                                             
  Germany                                    &_____         
                                             /     ^..^     
  T.n. 0431  880 7315                        \_____(oo)    
  Fax  04340 402 761                          WW WW        
  jlorenzo@tierzucht.uni-kiel.de                           
                                                           


          La vida est hecha solo de momentos

--
Asreml mailinglist archive: http://www.chiswick.anprod.csiro.au/lists/asreml