Gives FDR correction using mafdr always different values?

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.

Tags

Asked:

on 7 Nov 2012

Closed:

on 20 Aug 2021

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!