Metropolis step for second rate parameter

Metropolis step for second rate parameter

metropolisdRate(
  res,
  pilong,
  qlong,
  blong,
  currentd,
  currentSigma,
  bp,
  mean = 0,
  sd = 1,
  dshape = 1,
  drate = 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

blong

numeric vector of rate parameters

currentd

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

dshape

numeric shape parameter of the prior gamma distribution

drate

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