Pomic (version 0.9.1)

analyse.pomics: Analyse a data frame with parameter values and pomic scores

Description

This function calculate the POMIC weights values and do parameter averaging.

Usage

analyse.pomics(dataset,colnbparameters,colnbpomic,colnbvarparam=NULL,order=TRUE)

Arguments

dataset

Numerical, a dataset containing the parameter values and pomic scores

colnbparameters

Numerical, index of columns where the parameters are found

colnbpomic

Numerical, index of column where the pomic scores are found

colnbvarparam

Numerical, index of columns where the variances of parameters given each model can be found (should be of same length than colnbparameters)

order

Boolean, should the dataframe be ordered according to the pomic scores

Value

• Return a list containing:
• dataset
  an ordered data frame containing the entered data set and the columns of delta POMIC, Likelihood of model and POMIC weights
• params
  a data set of best parameter values, averaged parameter values and standard deviation of parameters

Details

The scores of POMIC should be raw values considering only one pattern. In case of multiple patterns fitting, for each model (or parameterization) the mean wPOMIC should be calculated and then the summarize.parameters function should be used. The dataset of results is ordered according to POMIC scores if order=TRUE (default).

See Also

summarize.parameters summarize.pomics pomic

Examples

dataset<-data.frame(P1=rep(1:5,each=5),P2=rep(seq(0,5,length=5),5),pomics=runif(25)*20)
analyse.pomics(dataset,c(1,2),3)