Metropolis step for the rate parameters

Metropolis step for the rate parameters

metropolisRate(
  res,
  pilong,
  qlong,
  dlong,
  currentb,
  currentSigma,
  bp,
  mean = 0,
  sd = 1,
  shape = 1,
  rate = 5,
  multivariate = FALSE,
  cormat = NULL,
  timepoints = timepoints,
  tCoef,
  index,
  numExch,
  density,
  phi
)

Arguments

res

data frame of uptake data

pilong

numeric vector of mixing parameters

qlong

numeric vector of stretch parameters

dlong

numeric vector of rate parameters

currentb

numeric vector of current rate parameters

currentSigma

numeric current variance parameter

bp

numeric vector of breakpoints

mean

numeric mean of the proposal distribution

sd

numeric standard deviation of the proposal distribution

shape

numeric shape parameter of the prior gamma distribution

rate

numeric rate parameter of the prior gamma distribution

multivariate

logical flag for multivariate model

cormat

numeric correlation matrix

timepoints

numeric vector of time points

tCoef

numeric vector of scaling variance of time points

index

list of indices for each peptide residues

numExch

list of number of exchangeable amides for each peptide

density

character string of density function

phi

numeric maximum uptake possible

Value

numeric updated rate parameters