I have been using LDA on several data sets. For some it seems to be successful and for other very similar data sets (all that has changed is my grouping variable of interest) the analysis does not work beyond the Effect Size calculator. I have pasted the error message below. Please let me know what the next best steps are.
Error Message:
/galaxy_venv/local/lib/python2.7/site-packages/rpy2/rinterface/init.py:185: RRuntimeWarning: Error in (function (file = “”, n = NULL, text = NULL, prompt = “?”, keep.source = getOption(“keep.source”), :
:1:47: unexpected numeric constant
1: effect.size <- abs(mean(LD[sub_d[,“class”]==""1
^
warnings.warn(x, RRuntimeWarning)
Traceback (most recent call last):
File “/shed_tools/testtoolshed.g2.bx.psu.edu/repos/george-weingart/lefse/a6284ef17bf3/lefse/run_lefse.py”, line 89, in
if params[‘rank_tec’] == ‘lda’: lda_res,lda_res_th = test_lda_r(cls,feats,class_sl,params[‘n_boots’],params[‘f_boots’],params[‘lda_abs_th’],0.0000000001,params[‘nlogs’])
File “/export/shed_tools/testtoolshed.g2.bx.psu.edu/repos/george-weingart/lefse/a6284ef17bf3/lefse/lefse.py”, line 199, in test_lda_r
robjects.r(‘effect.size <- abs(mean(LD[sub_d[,“class”]=="’+p[0]+’"]) - mean(LD[sub_d[,“class”]=="’+p[1]+’"]))’)
File “/galaxy_venv/local/lib/python2.7/site-packages/rpy2/robjects/init.py”, line 358, in call
p = _rparse(text=StrSexpVector((string,)))
rpy2.rinterface.RRuntimeError: Error in (function (file = “”, n = NULL, text = NULL, prompt = “?”, keep.source = getOption(“keep.source”), :
:1:47: unexpected numeric constant
1: effect.size <- abs(mean(LD[sub_d[,“class”]==""1
^
Thank you for any and all help!