CosmoPhotoz (version 0.1)

CosmoPhotoz-package: Photometric redshift estimation based on generalized linear models

Description

This package provides an user-friendly interfaces to perform fast and reliable photometric redshift estimation. The code makes use of generalized linear models and can adopt gamma or inverse gaussian families, either from a frequentist or a Bayesian perspective. The code additionally provides a Shiny application providing a simple user interface.

Details

Package:

CosmoPhotoz

Type:

Package

Version:

0.1

Date:

2014-08-22

License:

GPL (>= 3)

The CosmoPhotoz package aims to provide an user-friendly interface to enable the estimation of photometric redshifts. The present version employs generalized linear models and the user can adopt either gamma or inverse gaussian families from either a frequentist or a Bayesian perspective at the fitting step.

The package includes a plotting function to enable the production of diagnostic plots. Four examples of the implemented visual tests can be seen in the figures bellow:

merged-Diag-Plots.jpgmerged-Diag-Plots

Additionally, the code is also accompanied by a Shiny application that can be hosted locally or deployed by the user at a webserver. This application allows the user to run the photometric redshift estimation, to configure many parameters of the code visually and experiment with the results. It also allows the user to either use the PHAT0 data, or to upload its own data files (the expected format can be found at the application's help tab). A screenshot of this application can be seen in the next figure.

CosmoPhotoz-Shiny.jpgCosmoPhotoz-Shiny-Plots

To run the graphical interface locally, it is sufficient to call:

runApp(paste(find.package("CosmoPhotoz"),"/glmPhotoZ-2/",sep=""))

Finally, a short tutorial on how to use the package can be found at:

--------------- http://rafaelsdesouza.github.io/CosmoPhotoz/ -----------------

See Also

CosmoPhotoZestimator, glmTrainPhotoZ, glmPredictPhotoZ, glm, bayesglm

Examples

## Not run: 
# # Load the data
# data(PHAT0train)
# data(PHAT0test)
# # Run the analysis
# photoZest <- CosmoPhotoZestimator(PHAT0train, PHAT0test, 6)
# # This is considerably faster, but the results are not so good
# # photoZest <- CosmoPhotoZestimator(PHAT0train, PHAT0test, 6, robust=FALSE) 
# # Create a boxplot showing the results
# plotDiagPhotoZ(photoz = photoZest, specz = PHAT0test$redshift, type = "box")
# ## End(Not run)