[R] comparing columns and printing overlapping rows

2017-06-02 Thread Aanchal Sharma
Hi All, I have two files. 1. with only one column 2. data matrix I need to compare first columns of both files and print the rows from second file for the overlapping entries. I have solutions for awk and sed, but I need how to do it in R. Thanks Regards Anchal -- Anchal Sharma, PhD Postdoctora

[R] modEvA D-squared for gamma glm

2017-06-02 Thread Aanchal Sharma
Hi All, I am running a generalized linear model with gamma distribution in R (glm, family=gamma ) for my data (gene expression as response variable and few predictors). I want to calculate r-squared for this model. I have been reading online about it and found there are multiple formulas for calc

[R] ERROR : unable to start device PNG on server

2016-10-13 Thread Aanchal Sharma
Hi I realize that this question has been asked earlier also, but my problem is slightly different. I am running a rscript which runs for multiple datasets in loop and generates png for every file in the loop. When I run this only local machine (mac) it runs fine, just that that it takes very long

Re: [R] NaN Log-lik value in EM algorithm (fitting Gamma mixture model)

2016-09-16 Thread Aanchal Sharma
where TRUE/FALSE needed > > > Bill Dunlap > TIBCO Software > wdunlap tibco.com > > On Thu, Sep 15, 2016 at 3:04 PM, Aanchal Sharma < > aanchalsharma...@gmail.com> wrote: > >> I am using a function gammamixEM where it does it by default. I do not >> have >>

Re: [R] NaN Log-lik value in EM algorithm (fitting Gamma mixture model)

2016-09-15 Thread Aanchal Sharma
14, 2016 at 8:04 PM, Duncan Murdoch wrote: > On 14/09/2016 4:46 PM, Aanchal Sharma wrote: > >> Hi, >> >> I am trying to fit Gamma mixture model to my data (residual values >> obtained >> after fitting Generalized linear Model) using gammamixEM. It is part o

[R] NaN Log-lik value in EM algorithm (fitting Gamma mixture model)

2016-09-14 Thread Aanchal Sharma
Hi, I am trying to fit Gamma mixture model to my data (residual values obtained after fitting Generalized linear Model) using gammamixEM. It is part of the script which does it for multiple datasets in loop. The code is running fine for some datasets but it terminates for some giving following err

Re: [R] Fitting Mixture distributions

2016-09-13 Thread Aanchal Sharma
nal). > >> > >> For a while now that small package is providing an alternative > >> to the EM, namely direct MLE, simply using optim() where the > >> likelihood uses a somewhat smart parametrization. > >> > >> Of course, *as the EM*,

Re: [R] Fitting Mixture distributions

2016-09-12 Thread Aanchal Sharma
) > works considerably faster and more reliable than the EM (which I > also provide asnor1mix::norMixEM() . > > Apropos 'starting value': The help page shows how to use > kmeans() for "somewhat" reliable starts; alternatively, I'd > recommend using cluste

Re: [R] Fitting Mixture distributions

2016-09-07 Thread Aanchal Sharma
Hi Simon I am facing same problem as described above. i am trying to fit gaussian mixture model to my data using normalmixEM. I am running a Rscript which has this function running as part of it for about 17000 datasets (in loop). The script runs fine for some datasets, but it terminates when i