[Numpy-discussion] QQ Plot in Python and Noncentral Chisquare

[josef.pktd at gmail.com]

On Wed, Aug 15, 2012 at 10:59 AM, Yi, Chuang
<chuang.yi at bankofamerica.com> wrote:
> Hello,
>
>
>
> I am a new user of Python. I have a couple of questions that would
> appreciate your guidance!
>
>
>
> QQ Plot in Python: I could not find any functions in either Numpy or Scipy
> to do QQ Plot of two vectors of data. For example, in Matlab, one could just
> call qqplot(X,Y), which will generate the QQ plot of vector X against vector
> Y. Would you please let me know which python package has similar
> functionality?
> Noncentral Chisquare: it looks to me that in Numpy and Scipy, the noncentral
> chisquare distribution only allows Integer degree of freedom. For example,
> the following code will produce an error message. In R, however, it also
> allows non-integer degree of freedom. Would you please let me know which
> python package has similar functionality?
>
> data = numpy.random.noncentral_chisquare(0.5, 2, 100)
> ValueError: df <= 0

restriction is >1 not integer

>>> np.random.noncentral_chisquare(1, 2, 100)
Traceback (most recent call last):
  File "<stdin>", line 1, in <module>
  File "mtrand.pyx", line 1957, in
mtrand.RandomState.noncentral_chisquare
(numpy\random\mtrand\mtrand.c:9847)
ValueError: df <= 0

>>> np.random.noncentral_chisquare(1.0001, 2, 100)
array([  5.75083221,   1.08050491,   0.69267684,   1.37806056,
         2.39899   ,   0.31415666,   2.9202386


I don't know if the >1 restriction is justified or could be dropped.

Josef

>
>
>
> Thank you very much,
>
>
>
> Regards,
>
> Chuang
>
> ________________________________
> This message w/attachments (message) is intended solely for the use of the
> intended recipient(s) and may contain information that is privileged,
> confidential or proprietary. If you are not an intended recipient, please
> notify the sender, and then please delete and destroy all copies and
> attachments, and be advised that any review or dissemination of, or the
> taking of any action in reliance on, the information contained in or
> attached to this message is prohibited.
> Unless specifically indicated, this message is not an offer to sell or a
> solicitation of any investment products or other financial product or
> service, an official confirmation of any transaction, or an official
> statement of Sender. Subject to applicable law, Sender may intercept,
> monitor, review and retain e-communications (EC) traveling through its
> networks/systems and may produce any such EC to regulators, law enforcement,
> in litigation and as required by law.
> The laws of the country of each sender/recipient may impact the handling of
> EC, and EC may be archived, supervised and produced in countries other than
> the country in which you are located. This message cannot be guaranteed to
> be secure or free of errors or viruses.
>
> References to "Sender" are references to any subsidiary of Bank of America
> Corporation. Securities and Insurance Products: * Are Not FDIC Insured * Are
> Not Bank Guaranteed * May Lose Value * Are Not a Bank Deposit * Are Not a
> Condition to Any Banking Service or Activity * Are Not Insured by Any
> Federal Government Agency. Attachments that are part of this EC may have
> additional important disclosures and disclaimers, which you should read.
> This message is subject to terms available at the following link:
> http://www.bankofamerica.com/emaildisclaimer. By messaging with Sender you
> consent to the foregoing.
>
> _______________________________________________
> NumPy-Discussion mailing list
> NumPy-Discussion at scipy.org
> http://mail.scipy.org/mailman/listinfo/numpy-discussion
>