The hardware and bandwidth for this mirror is donated by METANET, the Webhosting and Full Service-Cloud Provider.
If you wish to report a bug, or if you are interested in having us mirror your free-software or open-source project, please feel free to contact us at mirror[@]metanet.ch.

How to find the best fitting distribution

library(peppwR)

The package has one function for estimating the fit of a range of distributions find_fits(). It takes a data frame with at least three columns:

• id_col a unique id for the phosphopeptide
• group_col a category that groups measurements from the same bio replicate together
• value_col the measurment abundance/counts etc

df_path <- system.file("raw_data.csv", package="peppwr")
df <- readr::read_csv(df_path)

fits_df <- find_fits(df, id_col="pep_id", group_col="group", value_col="value")

It returns a nested dataframe with elements for AIC and LogLikelihood of fit of each peptide’s counts to a range of distributions.

The result can be plot with evaldist()

evaldist(fits_df)

These binaries (installable software) and packages are in development.
They may not be fully stable and should be used with caution. We make no claims about them.