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