sampled_ple_analysis: Run the profile likelihood estimation analysis.

In pdp10/sbpiper: Data Analysis Functions for 'SBpipe' Package

Description

Run the profile likelihood estimation analysis.

Usage

sampled_ple_analysis(model, filename, parameter, plots_dir,
  fileout_param_estim_summary, logspace = TRUE, scientific_notation = TRUE)

Arguments

model	the model name
filename	the filename containing the fits sequence
parameter	the parameter to compute the PLE analysis
plots_dir	the directory to save the generated plots
fileout_param_estim_summary	the name of the file containing the summary for the parameter estimation
logspace	true if parameters should be plotted in logspace. (default: TRUE)
scientific_notation	true if the axis labels should be plotted in scientific notation (default: TRUE)

Examples

dir.create(file.path("pe_datasets"))
dir.create(file.path("pe_plots"))
data(insulin_receptor_all_fits)
write.table(insulin_receptor_all_fits, 
            file=file.path("pe_datasets", "all_fits.csv"), 
            row.names=FALSE)
# generate the global statistics for the parameter estimation
pe_ds_preproc(filename=file.path("pe_datasets", "all_fits.csv"), 
              param.names=c('k1', 'k2', 'k3'), 
              logspace=TRUE, 
              all.fits=TRUE, 
              data_point_num=33, 
              fileout_param_estim_summary=file.path("pe_datasets", "param_estim_summary.csv"))
sampled_ple_analysis(model="ir_beta", 
                     filename=file.path("pe_datasets", "all_fits_log10.csv"), 
                     parameter="k1", 
                     plots_dir="pe_plots", 
                     fileout_param_estim_summary=file.path("pe_datasets", 
                                                           "param_estim_summary.csv"),
                     logspace=TRUE)

pdp10/sbpiper documentation built on May 17, 2019, 11:17 p.m.