Rank: 27 based on 1452 downloads (last 30 days) and 60 files submitted
photo

Giuseppe Cardillo

E-mail
Company/University
MeriGen Research
Lat/Long
40.86195755004883, 14.21667289733887

Personal Profile:

I'm a biochemist with a PhD in Biotechnologies

Professional Interests:
Biochemestry and Molecular Biology

 

Watch this Author's files

 

Files Posted by Giuseppe View all
Updated   File Tags Downloads
(last 30 days)
Comments Rating
28 Apr 2014 Screenshot PrimeSpiral Play into the world of prime numbers! Author: Giuseppe Cardillo sacks spiral, vogel spiral, prime numbers 9 0
12 Feb 2014 MyFisher23 A very compact routine for Fisher's exact test on 2x3 contingency table Author: Giuseppe Cardillo statistics, probability, fishers exact test 2x... 20 10
  • 4.0
4.0 | 4 ratings
12 Feb 2014 MyFisher33 A very compact routine to compute Fisher's exact test on a 3x3 matrix Author: Giuseppe Cardillo statistics, probability, fisher exact test 3x3..., test, freemanhalton test, fisherfreemanhalton t... 20 4
  • 3.0
3.0 | 2 ratings
21 Jan 2014 Screenshot ROC curve compute a ROC curve Author: Giuseppe Cardillo statistics, probability, roc, receiver operating ch..., curve, by definition 143 50
  • 4.0
4.0 | 26 ratings
19 Nov 2013 Screenshot CODIS Combined DNA Index System: a GUI tool for forensic paternity lawsuit Author: Giuseppe Cardillo gui, biotech 24 4
Comments and Ratings by Giuseppe View all
Updated File Comments Rating
15 May 2014 Variable Precision Integer Arithmetic Arithmetic with integers of fully arbitrary size. Arrays and vectors of vpi numbers are supported. Author: John D'Errico

Dear John,
I saw that you use the Miller-Rabin algorithm to test primality. considering that are known the strong pseudoprimes to the first 9 prime bases (http://arxiv.org/pdf/1207.0063.pdf), I suggest a little implementation to speed up your Miller-Rabin algorithm:

if N<2047
witnesses=2;
elseif N<1373653
witnesses=[2 3];
elseif N<25326001
witnesses=[2 3 5];
elseif N<3215031751
witnesses=[2 3 5 7];
elseif N<2152302898747
witnesses=[2 3 5 7 11];
elseif N<3474749660383
witnesses=[2 3 5 7 11 13];
elseif N<341550071728321
witnesses=[2 3 5 7 11 13 17 19];
elseif N<3825123056546413051
witnesses=[2 3 5 7 11 13 17 19 23];
else
witnesses=[2 3 5 7 11 13 17 19 23 29 31 37];
end

Regards

03 Apr 2014 MyFisher23 A very compact routine for Fisher's exact test on 2x3 contingency table Author: Giuseppe Cardillo

For negative values, gammaln would return INF so that table is not added to computation. Maybe in previous version of gammaln this is not allowed and it calls error. So rows in Tables matrix with negative values must be erased before to give it to gammaln.
[I J]=find(Tables<0); Tables(I,:)=[];
will do this

03 Apr 2014 MyFisher23 A very compact routine for Fisher's exact test on 2x3 contingency table Author: Giuseppe Cardillo

try to add this line before line 108
[I ~]=find(Tables<0); Tables(I,:)=[];

03 Apr 2014 MyFisher23 A very compact routine for Fisher's exact test on 2x3 contingency table Author: Giuseppe Cardillo

I checked but on my pc the version is the same and the result is

2x3 matrix Fisher's exact test: 540 tables were evaluated
-----------------------------------------------------------------
p-value (2-tails): 0.3537569193
-----------------------------------------------------------------
Mid-p correction: 0.3472754799
-----------------------------------------------------------------

12 Feb 2014 MyFisher33 A very compact routine to compute Fisher's exact test on a 3x3 matrix Author: Giuseppe Cardillo

in my PC version it is correct. Bah! I uploaded it again

Comments and Ratings on Giuseppe's Files View all
Updated File Comment by Comments Rating
26 Jun 2014 ROC curve compute a ROC curve Author: Giuseppe Cardillo abdul

09 Apr 2014 ROC curve compute a ROC curve Author: Giuseppe Cardillo Gemmeke, Jort

Nice. One issue with the verbose option: when set to 0, it doesnt fill the .co part of the output anymore. To fix, switch the lines

m=[table(CSe,1) table(CSp,1) CE table(CEff,1)];
end

near the end of the code (line 300 or so)

09 Apr 2014 ROC curve compute a ROC curve Author: Giuseppe Cardillo Gemmeke, Jort

03 Apr 2014 MyFisher23 A very compact routine for Fisher's exact test on 2x3 contingency table Author: Giuseppe Cardillo Cardillo, Giuseppe

For negative values, gammaln would return INF so that table is not added to computation. Maybe in previous version of gammaln this is not allowed and it calls error. So rows in Tables matrix with negative values must be erased before to give it to gammaln.
[I J]=find(Tables<0); Tables(I,:)=[];
will do this

03 Apr 2014 MyFisher23 A very compact routine for Fisher's exact test on 2x3 contingency table Author: Giuseppe Cardillo Cardillo, Giuseppe

try to add this line before line 108
[I ~]=find(Tables<0); Tables(I,:)=[];

Contact us