Re: Tranlate ASreml call from R to Fortran

From: <arthur.gilmour_at_DPI.NSW.GOV.AU>
Date: Wed, 23 Apr 2008 16:09:34 +1000

Dear Hongsheng

re:

Does anybody know how to translate the following ASreml function call from
R to Fortran.

#######################################################################################################################################

asreml(Pheno ~ Geno, random = ~ ped(PID, var = T)+ units, rcov =
~PID:ar1(EID), ginverse = list(PID = com.inv) , data = dataset.long))

#######################################################################################################################################

If more information needed, please let us know.

I'm not sure what you mean by this question.
Are you wanting to run a procedure within R that you have written in
Fortran which invokes a call to ASReml?
I would not know how to do that. The asreml-r R functions call a C
interface to the ASReml core routines (which are in fortran).

Are you wanting to write a Fortran program to run outside of R but which
runs R and asreml() within R?

Otherwise, you can write a fortran program fairly easily which can run
standalone ASReml. The usual command
(depending on platform) is something like

      CALL SYSTEM("ASReml myjob")

May Jesus Christ be gracious to you in 2008,

Arthur Gilmour, His servant .
 
Mixed model regression mapping for QTL detection in experimental crosses.
Computational Statistics and Data Analysis 51:3749-3764 at
http://dx.doi.org/10.1016/j.csda.2006.12.031

Profile:http://www.dpi.nsw.gov.au/research/staff/arthur-gilmour
      http://www.dpi.nsw.gov.au/reader/17263
Personal website: http://www.cargovale.com.au/

Skype: arthur.gilmour
mailto:Arthur.Gilmour_at_dpi.nsw.gov.au, arthur_at_cargovale.com.au
Principal Research Scientist (Biometrics)
NSW Department of Primary Industries
Orange Agricultural Institute, Forest Rd, ORANGE, 2800, AUSTRALIA
 
fax: 02 6391 3899; 02 6391 3922 Australia +61
telephone work: 02 6391 3815; home: 02 6364 3288; mobile: 0438 251 426

ASREML 2 is now available from http://www.VSNi.co.uk/products/asreml
The ASReml discussion group is at ASREML-L_at_dpi.nsw.gov.au
To join it, mailto:arthur.gilmour_at_dpi.nsw.gov.au
Archives are at
https://gatekeeper.dpi.nsw.gov.au/Listserv/archives/asreml-l.html
Cookbook: http://uncronopio.org/ASReml

The ASReml3 alpha download registration page is:
http://www.austatgen.org/asreml3alpha.php
There are asreml3alpha forums at: http://www.austatgen.org/asreml3forum

Proposed travel:
Argentina 9-19 September ASReml WS + Biometrics 
><><><><><><><><><><><><><><><><><><><><><><><>
Received on Fri Apr 23 2008 - 16:09:34 EST

This webpage is part of the ASReml-l discussion list archives 2004-2010. More information on ASReml can be found at the VSN website. This discussion list is now deprecated - please use the VSN forum for discussion on ASReml. (These online archives were generated using the hypermail package.)