Wood's method by R program

--
0. Download 

Download a program file of R and data files:

  wood.R
  astro.txt
  mag.txt

The R package is available from
http://cran.r-project.org/

--
1. astronomical data

Start up R and type as follows.

  source("wood.R")
  astro <- list()
  astro$data <- read.table("astro.txt")
  astro$F <- F.optim(astro$data)

It takes about some seconds.
The result is checked by typing

  astro

The variable "astro$F$z" is the optimal point 
and "astro$F$likeinv" is the optimal value.

--
2. magnetism data

Change the file name "astro.txt" to "mag.txt" in the above.

The file "mag_make.R" was used to make the correct data "mag.txt"
(The data in Sec. 7 of Kent (1982) is wrong.)

--
3. Note

Monte Carlo method is also available

  astro$MC <- MC.optim(astro$data,tnum=1e6)

It takes about some minutes. The option "tnum" is the sampling number.

--
EOF