def run_msms_data(fragment, neutral_loss, mzdiff, ms1, ms2):
if len(sys.argv) > 1:
K = int(sys.argv[1])
else:
K = 300
print "Cross-validation for K=" + str(K)
n_folds = 4
n_samples = 500
n_burn = 250
n_thin = 5
alpha = 50.0 / K
beta = 0.1
is_num_samples = 10000
is_iters = 1000
ms2lda = Ms2Lda.lcms_data_from_R(fragment, neutral_loss, mzdiff, ms1, ms2)
df = ms2lda.df
vocab = ms2lda.vocab
cv = CrossValidatorLda(df, vocab, K, alpha, beta)
cv.cross_validate(n_folds,
n_burn,
n_samples,
n_thin,
is_num_samples,
is_iters,
method="with_mixture")
def run_msms_data(fragment, neutral_loss, mzdiff,
ms1, ms2):
if len(sys.argv)>1:
K = int(sys.argv[1])
else:
K = 300
print "Cross-validation for K=" + str(K)
n_folds = 4
n_samples = 500
n_burn = 250
n_thin = 5
alpha = 50.0/K
beta = 0.1
is_num_samples = 10000
is_iters = 1000
ms2lda = Ms2Lda.lcms_data_from_R(fragment, neutral_loss, mzdiff, ms1, ms2)
df = ms2lda.df
vocab = ms2lda.vocab
cv = CrossValidatorLda(df, vocab, K, alpha, beta)
cv.cross_validate(n_folds, n_burn, n_samples, n_thin,
is_num_samples, is_iters, method="with_mixture")
