LIRA Tutorial example

[Borlaff]

Dear LIRA team:

I am experiencing a possible bug following the tutorial described in page 13 from ira/tutorials/liratutorial_quick.pdf
I attach the required files from the tutorial example, plus the null_q1_c1.fits required file.
test_lira.zip

The error can be reproduced as follows:

source('liraOutput.R')
require(astro)

obsFile <- 'img_64x64_0.5.fits'
bkgFile <- 'null_q1_c1.fits'
psfFile <- 'psf_center_33x33.fits'
startFile <- 'null_q1_c1.fits'

testLira <- liraOutput(obsFile=obsFile, startFile=startFile,mapFile=FALSE,
bkgFile=bkgFile,
psfFile=psfFile,
maxIter=1000,
alpha.init <- c(3,4,5,6,7,8),
fit.bkg.scale=TRUE,
outDir='test_')

---

Output:
Loading required package: lira
Loading required package: FITSio
Loading required package: astro
Loading required package: MASS
Loading required package: plotrix

Code will run in posterior sampling mode.

A scale parameter will be fit to the bkg model.

The total number of Gibbs draws is 1000, every 1th draws will be saved.

The model will be fit using the Multi Scale Prior.

The data matrix is 64 by 64.

The data file should contain a 2^6 by 2^6 matrix of counts.

Starting Values for the smoothing parameter (alpha):
Aggregation Level: 0, alpha: 3 (Full Data)
Aggregation Level: 1, alpha: 4
Aggregation Level: 2, alpha: 5
Aggregation Level: 3, alpha: 6
Aggregation Level: 4, alpha: 7
Aggregation Level: 5, alpha: 8 (In the 2x2 table)

The prior distribution on the total count from the multiscale component is
Gamma(1.000000, 0.050000).

The hyper-prior smoothing parameter (kappa 2) is 1000.

Pixel: (3, 58)
Error in lira(obs.matrix = obsmat, start.matrix = strtmat, map.matrix = mapmat, :
The psf does not allow data in above pixel

---

Thanks a lot!