Gives FDR correction using mafdr always different values?
Info
This question is closed. Reopen it to edit or answer.
Show older comments
Hiya,
I ran mafdr on a column of p values, to check that I copy-pasted the result properly I re-ran. Then I noticed that pFDR values are different when I re-run after a short delay. It seems to iterate between the 3 results below, so is state depended, and if so, what state? And sorry for the display, I seem not to be able to get the values in single rows, although they show up as such in the message body (not in the preview though).
Thanks very much in advance for clarification.
Cheers,
Torsten
pvals=[ 0.280504000000000 0.0651080000000000 0.297689000000000 0.616244000000000 0.417778000000000 0.680408000000000 0.666234000000000 0.118566000000000 0.00645200000000000 0.299904000000000 0.247687000000000 0.606613000000000 0.640534000000000 0.0378090000000000 0.594998000000000 0.0452850000000000 0.509980000000000 0.00309100000000000 0.170900000000000 0.305498000000000 0.161217000000000 0.0581170000000000 0.902015000000000 0.654404000000000 0.520190000000000 0.494796000000000 0.582451000000000 0.203360000000000 0.00541900000000000 0.360566000000000 0.520378000000000 0.847790000000000 0.441881000000000 0.00530800000000000 0.645299000000000 0.00372100000000000 0.259240000000000 0.943877000000000 0.275483000000000 0.0672710000000000 0.234132000000000 0.0103720000000000 0.876059000000000 0.139472000000000 ]
1st run: PVALS=mafdr(pvals) PVALS= 0.168168204694041 0.0819706688302527 0.170358665835012 0.221670950432382 0.202300214097318 0.214157625934811 0.215073206086985 0.124395034923099 0.0162460759136447 0.164164234575248 0.173242655370554 0.224624394811544 0.224008145391406 0.0680019356366789 0.226999912558633 0.0712670051679132 0.221400843214628 0.0389155460702696 0.143441536145995 0.160258566447957 0.144979608790744 0.0812989248702634 0.264100624382673 0.216813576006476 0.218305595926815 0.222480664294786 0.229157217771095 0.160018308169889 0.0170562394172429 0.181580366022295 0.211339831904851 0.260332600570609 0.206046711636852 0.0222758242791967 0.219575253032771 0.0234236083674334 0.171779975195503 0.270076533617892 0.173415907118668 0.0769944325076646 0.173394763402751 0.0217638328394714 0.262608143263656 0.135072768123650
2nd run: PVALS= 0.217133877645733 0.105838135149918 0.219962137138880 0.286215061380766 0.261204131990530 0.276514076077769 0.277696246447601 0.160615238426732 0.0209764590522600 0.211964068290553 0.223685860253372 0.290028462562192 0.289232779297124 0.0878021145499622 0.293095661743988 0.0920178770324700 0.285866306824710 0.0502466172741287 0.185207524903453 0.206921183597700 0.187193439411268 0.104970798954691 0.340999018008195 0.279943361292197 0.281869813866376 0.287260540287121 0.295881111269952 0.206610969121377 0.0220225185222018 0.234451268897619 0.272875822668064 0.336133855637796 0.266041500266116 0.0287618941541114 0.283509158064634 0.0302438794689474 0.221797290301982 0.348714937596341 0.223909557853911 0.0994129640495253 0.223882257705032 0.0281008257504186 0.339071969947743 0.174402015888741
3rd run: PVALS= 0.297645330649711 0.145082043730845 0.301522285462696 0.392341248198031 0.358056472238936 0.379043217464892 0.380663726868713 0.220170045629203 0.0287543572572161 0.290558598576381 0.306626734430671 0.397568627121358 0.396477911056727 0.120358415279657 0.401773118491073 0.126137347759709 0.391863178325754 0.0688776490096520 0.253880949331382 0.283645853784518 0.256603224578016 0.143893106419855 0.467438644612951 0.383744053385144 0.386384818703018 0.393774382069831 0.405591389753729 0.283220614336032 0.0301882869607978 0.321383867910892 0.374055929664110 0.460769520116568 0.364687496819129 0.0394265675556166 0.388632018224483 0.0414580608160652 0.304037898300025 0.478015563617859 0.306933377258589 0.136274472119299 0.306895954439413 0.0385203804339540 0.464797063007400 0.239068846593688
Answers (0)
This question is closed.
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!