analyse.pomics {Pomic} | R Documentation |
This function calculate the POMIC weights values and do parameter averaging.
analyse.pomics(dataset,colnbparameters,colnbpomic,colnbvarparam=NULL,order=TRUE)
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 |
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).
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 |
Cyril Piou
cyril.piou@yahoo.fr
summarize.parameters
summarize.pomics
pomic
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)