[R] S-Plus resample package and associated functions

2007-09-10 Thread Robert A. LaBudde
Are there any packages in R that reproduce the package resample of S-Plus? The sample() function in R doesn't provide equivalent flexibility of bootstrap() and bootstrap2(). Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e-mail: [EMAIL

Re: [R] Excel

2007-08-28 Thread Robert A LaBudde
At 01:21 AM 8/28/2007, David wrote: On Tue, 28 Aug 2007, Robert A LaBudde wrote: If you format the column as Text, you won't have this problem. By leaving the cells as General, you leave it up to Excel to guess at the correct interpretation. Not true actually. I had converted the column to Text

Re: [R] Excel

2007-08-27 Thread Robert A LaBudde
. Clicking on that entry it showed 1/08/2699. In a column of character data, Excel had interpreted AUG2699 as a date. The .csv did not actually have a date in that cell, but if I had saved the .csv file it would have. David Scott Robert

Re: [R] chi-square statistics

2007-08-23 Thread Robert A LaBudde
the corresponding pvalues for each values in x. ? pchisq E.g. pchisq(x, df, lower.tail=FALSE) Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e-mail: [EMAIL PROTECTED] Least Cost Formulations, Ltd.URL: http://lcfltd.com/ 824

[R] Conditional logistic regression on n:m matched cohort data

2007-07-22 Thread Robert A. LaBudde
. Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e-mail: [EMAIL PROTECTED] Least Cost Formulations, Ltd.URL: http://lcfltd.com/ 824 Timberlake Drive Tel: 757-467-0954 Virginia Beach, VA 23464-3239Fax: 757-467-2947 Vere scire est per causas scire

[R] Conditional logistic regression on n:m matched cohort data [Corrected]

2007-07-22 Thread Robert A. LaBudde
? Thanks. Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e-mail: [EMAIL PROTECTED] Least Cost Formulations, Ltd.URL: http://lcfltd.com/ 824 Timberlake Drive Tel: 757-467-0954 Virginia Beach, VA 23464-3239

Re: [R] Effect size, interactions, and main effects (Stats Question) [ffmanova]

2007-07-21 Thread Robert A LaBudde
the main question is fair, are there general guidlines concerning the relationship between level of significance and effect size for interaction terms. Thank you in advance for your help, Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e

Re: [R] No convergence using ADAPT

2007-07-07 Thread Robert A LaBudde
guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e-mail: [EMAIL PROTECTED] Least Cost Formulations, Ltd.URL: http

Re: [R] outlying

2007-06-19 Thread Robert A LaBudde
Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e-mail: [EMAIL PROTECTED] Least Cost Formulations, Ltd.URL: http://lcfltd.com/ 824 Timberlake Drive Tel: 757-467-0954 Virginia Beach, VA 23464-3239Fax: 757-467-2947 Vere scire est per causas scire

Re: [R] Optimization

2007-06-18 Thread Robert A LaBudde
. Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e-mail: [EMAIL PROTECTED] Least Cost Formulations, Ltd.URL: http://lcfltd.com/ 824 Timberlake Drive Tel: 757-467-0954 Virginia Beach, VA 23464-3239Fax: 757-467-2947 Vere scire est per causas scire

Re: [R] Iterative Solver [Converting Matlab's solve()]

2007-06-18 Thread Robert A LaBudde
place a ) after y). 3. To find numerical roots, see uniroot(). 4. If you want more help, you'll have to specify domains or values for x and y. Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e-mail: [EMAIL PROTECTED] Least Cost Formulations

Re: [R] [Not R question]: Better fit for order probit model

2007-06-17 Thread Robert A LaBudde
model, you would need to include a full set of interaction terms with class. Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e-mail: [EMAIL PROTECTED] Least Cost Formulations, Ltd.URL: http://lcfltd.com/ 824 Timberlake Drive

Re: [R] [Not R question]: Better fit for order probit model

2007-06-16 Thread Robert A LaBudde
model are you fitting? Etc. Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e-mail: [EMAIL PROTECTED] Least Cost Formulations, Ltd.URL: http://lcfltd.com/ 824 Timberlake Drive Tel: 757-467-0954 Virginia Beach

Re: [R] [Not R question]: Better fit for order probit model

2007-06-15 Thread Robert A LaBudde
utilized. And don't use so many categories. People have trouble even selecting responses on a 5-level scale. Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e-mail: [EMAIL PROTECTED] Least Cost Formulations, Ltd.URL: http://lcfltd.com

Re: [R] Appropriate regression model for categorical variables

2007-06-12 Thread Robert A LaBudde
. Consider a change to white, nonwhite. Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e-mail: [EMAIL PROTECTED] Least Cost Formulations, Ltd.URL: http://lcfltd.com/ 824 Timberlake Drive Tel: 757-467-0954 Virginia

Re: [R] pretty report

2007-06-12 Thread Robert A LaBudde
snip Note that files with a .csv extension are also associated with Excel and can opened with a double-click. Comma-separated-value files also can be unambiguously loaded by Excel without parsing. Robert A. LaBudde, PhD, PAS, Dpl

Re: [R] p-value from GEE why factor 2*pnorm?

2007-06-11 Thread Robert A LaBudde
Wald confidence interval, which obviously assumes a normal distribution applies. Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e-mail: [EMAIL PROTECTED] Least Cost Formulations, Ltd.URL: http://lcfltd.com/ 824 Timberlake Drive

Re: [R] Windows vista's early terminate Rgui execution

2007-06-10 Thread Robert A LaBudde
with. Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e-mail: [EMAIL PROTECTED] Least Cost Formulations, Ltd.URL: http://lcfltd.com/ 824 Timberlake Drive Tel: 757-467-0954 Virginia Beach, VA 23464-3239Fax: 757-467-2947

Re: [R] What ECDF function?

2007-06-09 Thread Robert A LaBudde
are not very accurate, so the differences between these formulae are difficult to justify in practice. Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e-mail: [EMAIL PROTECTED] Least Cost Formulations, Ltd.URL: http://lcfltd.com/ 824

Re: [R] What ECDF function?

2007-06-09 Thread Robert A LaBudde
At 06:36 PM 6/9/2007, Marco wrote: On 6/9/07, Robert A LaBudde [EMAIL PROTECTED] wrote: At 12:57 PM 6/9/2007, Marco wrote: snip snip Hmmm I'm a bit confused, but very interested! So you don't use the R ecdf, do you? Only when an i/n edf is needed (some tests, such as ks.test() are based

Re: [R] pnorm how to decide lower-tail true or false

2007-06-08 Thread Robert A LaBudde
should use pnorm(z, lower.tail=FALSE) instead of 1-pnorm(z) because the former returns a more accurate answer for large z. This is really simple issue, and has no inherent complexity associated with it. Robert A. LaBudde, PhD, PAS, Dpl

[R] glm() for log link and Weibull family

2007-06-08 Thread Robert A. LaBudde
(link='log')) but this doesn't quite do the trick. The standardized deviance residuals are still curved away from normal at the tails. Thanks for any info you can give on this nonstandard model. Robert A. LaBudde, PhD, PAS, Dpl. ACAFS

Re: [R] Finding density curve peak

2007-06-02 Thread Robert A LaBudde
. Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e-mail: [EMAIL PROTECTED] Least Cost Formulations, Ltd.URL: http://lcfltd.com/ 824 Timberlake Drive Tel: 757-467-0954 Virginia Beach, VA 23464-3239Fax: 757-467-2947 Vere scire est per causas scire

[R] How to reference or sort rownames in a data frame

2007-05-27 Thread Robert A. LaBudde
. Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e-mail: [EMAIL PROTECTED] Least Cost Formulations, Ltd.URL: http://lcfltd.com/ 824 Timberlake Drive Tel: 757-467-0954 Virginia Beach, VA 23464-3239Fax: 757-467-2947 Vere scire est per causas scire

Re: [R] How to reference or sort rownames in a data frame

2007-05-27 Thread Robert A LaBudde
in R. At 10:29 PM 5/27/2007, Gabor wrote: On 5/27/07, Robert A. LaBudde [EMAIL PROTECTED] wrote: As I was working through elementary examples, I was using dataset plasma of package HSAUR. In performing a logistic regression of the data, and making the diagnostic plots (R-2.5.0) data(plasma

[R] Curve crosses back to origin in plot

2007-05-27 Thread Robert A. LaBudde
to be determined by inspection. Thanks again for answering these simple questions that seem to be the hardest to find answers for. Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e-mail: [EMAIL PROTECTED] Least Cost Formulations, Ltd.URL

Re: [R] Log-likelihood function

2007-05-02 Thread Robert A LaBudde
sure you are performing the same fit in both cases. Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e-mail: [EMAIL PROTECTED] Least Cost Formulations, Ltd.URL: http://lcfltd.com/ 824 Timberlake Drive Tel: 757

Re: [R] upgrade to 2.5

2007-05-02 Thread Robert A LaBudde
. Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e-mail: [EMAIL PROTECTED] Least Cost Formulations, Ltd.URL: http://lcfltd.com/ 824 Timberlake Drive Tel: 757-467-0954 Virginia Beach, VA 23464-3239Fax: 757

Re: [R] Is R's fast fourier transform function different from fft2 in Matlab?

2007-05-02 Thread Robert A LaBudde
= 6.7500-4.2500 -5.2500 3.7500 I also tried mvfft with inverse but can't get same result with ifft2. Does any function work? Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e-mail: [EMAIL PROTECTED] Least Cost Formulations

Re: [R] Survival statistics--displaying multiple plots

2007-05-02 Thread Robert A LaBudde
. I am unsure about how to go about this. Any suggestions would be appreciated. Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e-mail: [EMAIL PROTECTED] Least Cost Formulations, Ltd.URL: http://lcfltd.com/ 824 Timberlake Drive

Re: [R] Survival statistics--displaying multiple plots

2007-05-02 Thread Robert A LaBudde
),subset(tips,meld10 meld 20)) And for patients with meld20, plot(survfit(Surv(days,status==1),subset(tips,meld20)) But how do I display the curves in each cadre (meld10, 10meld20, and meld20) on the same chart? -Original Message- From: Robert A LaBudde [mailto:[EMAIL PROTECTED] Sent

Re: [R] bivariate normal distribution in likelihood

2007-04-30 Thread Robert A LaBudde
. Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e-mail: [EMAIL PROTECTED] Least Cost Formulations, Ltd.URL: http://lcfltd.com/ 824 Timberlake Drive Tel: 757-467-0954 Virginia Beach, VA 23464-3239Fax: 757-467-2947 Vere scire est per causas scire

Re: [R] How to solve difficult equations?

2007-04-25 Thread Robert A LaBudde
this numerically? Is it a single instance of a larger, more general problem? Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e-mail: [EMAIL PROTECTED] Least Cost Formulations, Ltd.URL: http://lcfltd.com/ 824 Timberlake Drive

[R] Use of Lexis function to convert survival data to counting format

2007-04-24 Thread Robert A. LaBudde
something to do with 'start' and 'stop', of which arguments I am apparently clueless as to the meaning. Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e-mail: [EMAIL PROTECTED] Least Cost Formulations, Ltd.URL: http://lcfltd.com/ 824

Re: [R] Using R to create pdf's from each file in a directory

2007-04-20 Thread Robert A LaBudde
will amplify gecko951's question to include How you you concatenate two strings in R?. I.e., x-abc y-def What operator or function returns abcdef as a result? Thanks. Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e-mail: [EMAIL PROTECTED] Least

Re: [R] Using R to create pdf's from each file in a directory

2007-04-20 Thread Robert A LaBudde
. Replace pdf(/tmp/graph/x.pdf) with pdf(paste(paste('/tmp/graph/',x,sep=''),'.pdf',sep='')) and the program should now work for the different file names. Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e-mail: [EMAIL PROTECTED] Least