Search Comments and Ratings

go

Comments and Ratings

   
Date File Comment by Comment Rating
27 Oct 2009 ROC curve compute a ROC curve Author: Giuseppe Cardillo Kaftan, Jens

Hi Giuseppe.
I have had a look at the new release today and I think it is still not perfectly correct. I have validated the scripts using the example data of Hanley and McNeil's 1982 paper: "The Meaning and Use of the Area under a Receiver Operating Characteristic (ROC) Curve", which seems to be the basis for the calculations (such as the approximation of Q_1 and Q_2) anyways. To my opinion the problem is that when integrating over the ROC curve to compute the AUC, the data point (sensitivity=1, specificity=0) is not considered when using the trapezoidal rule. Consequently the AUC value (and all AUC dependent measures) differ slightly from the example in the mentioned article (which becomes more severe for non-continous tests with only a few cut-off points).

Best,
Jens.

23 Oct 2009 MyFisher22 A very compact routine for Fisher's exact test on 2x2 matrix; power and sample sizes calculation Author: Giuseppe Cardillo Giaccari, Luigi

20 Oct 2009 McNemar est Perform the McNemar test on a 2x2 matrix Author: Giuseppe Cardillo Skynet

Thanks for this file. I used it to confirm the results of my own calculations. Note that errors that you throw probably shouldn't use the term "Warning" in the error text, because that is misleading.

14 Oct 2009 ISLEAP Function Calculate leap year Author: Giuseppe Cardillo Simon, Jan

Descriptive help, reliable error checking, new version is really efficient, does exactly what the name says.

10 Oct 2009 Clinical Test Performance The test calculate the performance of a clinical test based on the Bayes theorem Author: Giuseppe Cardillo Simon, Jan

Exhaustively commented, so you get not only the resulting numbers, but also the exact description of what they mean, if you read the source.

05 Oct 2009 ISLEAP Function Calculate leap year Author: Giuseppe Cardillo Simon, Jan

CLOCK -> DATE -> DATESTR -> STR2DOUBLE waste a lot of time: str2double(datestr(date,10))
Faster: c = clock; default = {c(1)};

Checking the number of inputs can be done easier and faster without VARARGIN, e.g.:
  function x = isleap(Year)
  switch nargin
  case 1, % Your error checks...
  case 0, c = clock; Year = c(1);
  otherwise error(...);
  end

The permanent disabling of the warning for logical conversion is not secure. You can omit this by simplifying the engine - instead of:
  warning('OFF','MATLAB:conversionToLogical')
  A=logical(mod(Year,4));
  B=logical(mod(Year,100));
  C=logical(mod(Year,400));
  x=ismember(Year,Year(~C | (~A & B)));
call just:
  x = ~mod(Year, 4) & (mod(Year, 100) | ~mod(Year, 400));
And suddenly the function is 150 times faster (400 times if called without argument).
Good help, H1-line, date, author and reference included. I'll rate this, if its efficiency is encreased.

25 Sep 2009 Fleiss'es kappa Compute the Fleiss'es kappa for multiple raters Author: Giuseppe Cardillo ADEWUNMI, Adrian

With much pleasure, I am pleased to announce I have solved the problem .

Thanks to Giuseppe Cardillo for this matlab function.....good job.

25 Sep 2009 Fleiss'es kappa Compute the Fleiss'es kappa for multiple raters Author: Giuseppe Cardillo ADEWUNMI, Adrian

Whenever I imput any other matrix than a 5 x 10 matrix into matlab, using your function "fleiss(X)"it gives an error message as follows:

EDU>> fleiss(X)
??? Error using ==> fleiss at 107
The raters are not the same for each rows

Can you tell me how to fix this?
 Thx

09 Sep 2009 Logrank Comparing survival curves of two groups using the log rank test Author: Giuseppe Cardillo Cardillo, Giuseppe

At the moment, there is not on FEX the new uploaded file (this takes 1-2 days). If you contact me in private, I'll send you the file as an attach

09 Sep 2009 Logrank Comparing survival curves of two groups using the log rank test Author: Giuseppe Cardillo cyclist, the

Thanks for the really quick reply. I appreciate it.

However, the program still crashes for me, with the same error, when I try the syntax you suggest. Maybe the version on the File Exchange is not quite the same as yours?

09 Sep 2009 Logrank Comparing survival curves of two groups using the log rank test Author: Giuseppe Cardillo Cardillo, Giuseppe

You are partially right. If you use only the first column, the routine crashes because the informations in the second column are mandatory. If you want no censored data:
>> x1(:,2)=0; x2(:,2)=0; logrank(x1,x2);

Anyway, a bug was present but I fixed it and upload the new version.

08 Sep 2009 Logrank Comparing survival curves of two groups using the log rank test Author: Giuseppe Cardillo cyclist, the

It seems to me that this function does not work for the case when none of the data are censored. (That shouldn't be a problem, should it?) For example, using only the first columns of your sample x1 and x2:

>> logrank(x1(:,1),x2(:,1))

causes a crash for me.

04 Sep 2009 MyRegression A simple function on LS linear regression with many informative outputs Author: Giuseppe Cardillo Ben

21 Aug 2009 Testt Calculate the Student t Test for unequal or equal samples size, unpaired or paired samples. Author: Giuseppe Cardillo Ben

Thanks

04 Aug 2009 Testt Calculate the Student t Test for unequal or equal samples size, unpaired or paired samples. Author: Giuseppe Cardillo Chandler, john

can the p-value, etc be passed onto variables?

29 Jul 2009 ROC curve compute a ROC curve Author: Giuseppe Cardillo cabrego

I tested the new release and it is agreeing with other codes now. Michael, you may also wish to verify that the new version is working correctly.

I also think adding the cut off points as an additional x or y axis would be very useful to understand the trade off between sensitivity and sensibility.

24 Jul 2009 ROC curve compute a ROC curve Author: Giuseppe Cardillo Michael

Agree with cabrego, this algorithm does not work correctly. Depending on the input data, it generates ROC curves with specificity and sensitivity backward. I believe this is because elements that fall below a cutoff value (I in the code) are called "true positives" when they should be "false positives". The convention is that higher values of a test are abnormal (positive).

I confirmed that other software (online ROC calculator, ROCR in R, STATA) does not behave this way with the same input data and all others produce correct results.

Use at your own risk.

16 Jul 2009 ROC curve compute a ROC curve Author: Giuseppe Cardillo cabrego

Nice function, but I think it may have a bug, I am getting different results for significantly overlaping distributions when I compare to medcalc, and online calculators.

email me for more details: cpabrego@gmail.com

25 Jun 2009 ROC curve compute a ROC curve Author: Giuseppe Cardillo Vo, Phong

Thank you very much!

12 Jun 2009 Logrank Comparing survival curves of two groups using the log rank test Author: Giuseppe Cardillo Bl, En

Thanks.

08 Jun 2009 Logrank Comparing survival curves of two groups using the log rank test Author: Giuseppe Cardillo Cardillo, Giuseppe

I have corrected the bug and uploaded the new version. Anyway, with all these censored data, it is difficult to assess differences....

01 Jun 2009 Equinoxes and Solstices compute the dates of equinoxes and solstices Author: Giuseppe Cardillo Burkey, Jeff

Great function! While the level of accuracy does not affect my needs, you may want to let others know if they are in need to to the minute or second, may want to validate the results. Comparing this functions results to the US Naval Observatory, the results differ by approximately 1-minute. I'm guessing this is an artifact of matlab and complier precision. Anything dealing with astronomy and floating point precision becomes significant.

http://aa.usno.navy.mil/data/docs/EarthSeasons.php

30 May 2009 Logrank Comparing survival curves of two groups using the log rank test Author: Giuseppe Cardillo Bl, En

Here is the data:
x1 =
3010 1
1733 1
2797 1
334 0
3678 1
3314 1
4378 1
3253 1
1672 0
3709 1
1642 1
3040 1
2827 1
 and x2 =
1955 1
1808 1
2037 1
3136 1
938 0
3497 1
2858 1
3832 1
3256 0
2933 0
5054 1
927 0
3795 0
4155 1
3304 0

sorry for the late reply.

22 May 2009 Logrank Comparing survival curves of two groups using the log rank test Author: Giuseppe Cardillo Cardillo, Giuseppe

Perhaps, if you email me the data I can understand the error. I think that your input data matrix are in a wrong format.

19 May 2009 Logrank Comparing survival curves of two groups using the log rank test Author: Giuseppe Cardillo Bl, En

I get this error:
Matrix dimensions must agree.

Error in ==> logrank at 157
    K=find((table12(:,1)<table(c(I),1) & table12(:,1)>=table(J,1)),1,'last');

I debugged it and J=[] .

15 Apr 2009 Testt Calculate the Student t Test for unequal or equal samples size, unpaired or paired samples. Author: Giuseppe Cardillo Cardillo, Giuseppe

You have partially right. I have fixed the bug and uploaded the new file. Thank you very much!

14 Apr 2009 Testt Calculate the Student t Test for unequal or equal samples size, unpaired or paired samples. Author: Giuseppe Cardillo Berg, Rune W

There's a small bug in the procedure that prevents you chose 2-tailed test. Its in line 109.

As far as i can see you need to change it to:

109: if tail ~= 0 && tail ~= 1 && tail ~= 2 %check if tst is 0 or 1 or 2

02 Feb 2009 Testt Calculate the Student t Test for unequal or equal samples size, unpaired or paired samples. Author: Giuseppe Cardillo keihaninejad, shiva

this m file is good,easy and useful.many thanks Giuseppe

20 Jan 2009 ISLEAP Function Calculate leap year Author: Giuseppe Cardillo Petaccio, Vince

Cool, just what I needed. Thanks! :o]

06 Jan 2009 Sequential analysis Perform a sequential analysis to test two conditions. Author: Giuseppe Cardillo Zhang, Kay

I cannot use it

12 Dec 2008 MyFisher24 A very compact routine to compute Fisher's exact test on a 2x4 matrix Author: Giuseppe Cardillo Jon

Giuseppe thank you for your quick response.
Your routine takes all my data now.
Jon

12 Dec 2008 MyFisher24 A very compact routine to compute Fisher's exact test on a 2x4 matrix Author: Giuseppe Cardillo Cardillo, Giuseppe

Jon was right. I have fixed the bug and uploaded the new version

10 Dec 2008 MyFisher24 A very compact routine to compute Fisher's exact test on a 2x4 matrix Author: Giuseppe Cardillo Jon

I get error messages for some data.
For instance:
y =
    27 13 9 15
    12 10 3 16
>> myfisher24(y)
??? Attempted to access idxBstart(15); index out of bounds because
numel(idxBstart)=14.
Error in ==> myfisher24 at 116
            Tables(idxBstart(J+1):idxBstop(J+1),2)=J;

I hope this can be fixed?
Jon

28 Nov 2008 Cohen's kappa Compute the Cohen's kappa Author: Giuseppe Cardillo Karelitz, tzur

I think there is an error in the quadratic loop. small caps j needs to be J.

20 Nov 2008 ROC curve compute a ROC curve Author: Giuseppe Cardillo Pawel

17 Aug 2008 The chaotic rhythm of life Explore the May-Oster-Yorke law Author: Giuseppe Cardillo sanchez, nathaly

good

11 Aug 2008 ROC curve compute a ROC curve Author: Giuseppe Cardillo Eibl, Günther

easy to use!

13 Jul 2008 ROC curve compute a ROC curve Author: Giuseppe Cardillo García Bascuñán, Diego

Good function

13 Jul 2008 Testt Calculate the Student t Test for unequal or equal samples size, unpaired or paired samples. Author: Giuseppe Cardillo García Bascuñán, Diego

If you get the powerStudent, testt works. It's fantastic do it so easily!!!

24 Jun 2008 ROC curve compute a ROC curve Author: Giuseppe Cardillo Phan, Truc

10 Jun 2008 Dunn's test Dunn's procedure for multiple, non parametric, comparisons Author: Giuseppe Cardillo Makhlof, Abdallah

thank you

04 Jun 2008 ROC curve compute a ROC curve Author: Giuseppe Cardillo HAMADICHAREF, Brahim

Good stuff.

07 May 2008 MyFisher23 A very compact routine for Fisher's exact test on 2x3 contingency table Author: Giuseppe Cardillo Gandhi, Shilpa

Fast, provides good output display. Exactly what I needed. Thanks!

27 Mar 2008 Coins Is there an algorithm to use the max numbers of coins in your pocket? Author: Giuseppe Cardillo Shvorob, Dimitri

Is there an algorithm to use the max numbers of coins? Of course there is, as anybody who has taken half a course in dynamic programming will tell you. Your own code solves the problem in less than 10 lines :)

13 Mar 2008 MyRegression A simple function on LS linear regression with many informative outputs Author: Giuseppe Cardillo C, C

It might already be in stat toolbox, but for those of use that can't afford the toolbox, this is great.

07 Mar 2008 Wilcoxon test non parametric Wilcoxon test for two paired samples Author: Giuseppe Cardillo Trujillo-Ortiz, Antonio

Only to correct and clarify my previous comment, and to recall to all that, with respect to the non parametric tests between two samples, two classes exist: (1) Independent, that can be both of the same size or not, being been able to utilize the Mann-Whitney’s test or the Wilcoxon ranks test; and (2) Dependent, that necessarily them should be of the same size, since is the same sample observed in two different circumstances or times, where applies the Wilcoxon signs and ranks test.

Also, I should clarify that in the personal thing I do not use the term unpaired to assign independent samples, since unpaired is synonym of unbalanced; neither paired for assign dependent, since paired is synonym of balanced. This by the fundamental reason that the concepts of independent or dependent arise for the itself nature of the data, and not because they are unpaired or paired.

The m-file here developed by Giuseppe it is a Wilcoxon’s test for two dependent samples for it is correct his comment that the data vectors must have the same length, and, if it is not, delete those uncoupled data. There isn't any other solutions.

So, my previous review comment is obsolete and should be ignored.

Thanks,

Prof. Antonio Trujillo-Ortiz

05 Mar 2008 Wilcoxon test non parametric Wilcoxon test for two paired samples Author: Giuseppe Cardillo Trujillo-Ortiz, Antonio

Hi Giuseppe,

You are wrong. According to the non-parametric statistics theoretic fundamentals, there are a two Wilcoxon tests:(1)Wilcoxon Rank Test for two unpaired samples, and(2) Wilcoxon Sign-Rank Test for two paired samples. You must to review them.

Yours,

Prof. Antonio Trujillo-Ortiz

19 Jan 2008 Dunn's test Dunn's procedure for multiple, non parametric, comparisons Author: Giuseppe Cardillo Akkus, Ozan

The code is writen very clearly, instructions are detailed and to the point. It did not take me more than 5 minutes to get run the test.

27 Dec 2007 Holm-Sidak t-test Holm-Sidak t-test for multiple comparisons Author: Giuseppe Cardillo Norrell, Gilbert

This file is similar to multcompare matlab function, but provides much more clear outputs.

22 Dec 2007 Tolkiencal Convert a date into the Shire calendar Author: Giuseppe Cardillo Observer, A FEX

I agree with Urs. This submission does not provide anything valuable or useful (in terms of contents or coding) to the Matlab community. Yet it is rated amazingly high by a few of those people, who are usually very quick in annihilating such surrealistic contributions. The fact that a guru(?) has an emotional bond with the idea behind a program certainly cannot be sufficient to give out a 'reward'.
Walter.

 

MATLAB Central Terms of Use

NOTICE: Any content you submit to MATLAB Central, including personal information, is not subject to the protections which may be afforded information collected under other sections of The MathWorks, Inc. Web site. You are entirely responsible for all content that you upload, post, e-mail, transmit or otherwise make available via MATLAB Central. The MathWorks does not control the content posted by visitors to MATLAB Central and, does not guarantee the accuracy, integrity, or quality of such content. Under no circumstances will The MathWorks be liable in any way for any content not authored by The MathWorks, or any loss or damage of any kind incurred as a result of the use of any content posted, e-mailed, transmitted or otherwise made available via MATLAB Central. Read the complete Terms prior to use.

Contact us at files@mathworks.com