Re: [R] significance testing for the difference in the ratio of means

__ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. Robert

[R] Mininum number of resamples required to do BCa bootstrap?

. Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e-mail: r...@lcfltd.com 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

[R] Possible error in BCa method for confidence intervals in package 'boot'

0.001/n in the calculations makes me suspect the programming in the script. The script apparently works correctly if the number of resamples equals or exceeds the number of original data, but not otherwise. Robert A. LaBudde, PhD

Re: [R] coxph() - unexpected result using Crawley's seedlings data (The R Book)

. Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e-mail: r...@lcfltd.com Least Cost Formulations, Ltd.URL: http://lcfltd.com/ 824 Timberlake Drive Tel: 757-467-0954 Virginia Beach, VA 23464-3239Fax: 757

Re: [R] analysing a three level reponse

/mailman/listinfo/r-help PLEASE do read the posting 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: r...@lcfltd.com

Re: [R] Using MCMC sampling to estimate p values with a mixed model

, reproducible code. Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e-mail: r...@lcfltd.com Least Cost Formulations, Ltd.URL: http://lcfltd.com/ 824 Timberlake Drive Tel: 757-467-0954 Virginia Beach, VA 23464

[R] How to convert a factor column into a numeric one?

of the sort order in the recode() 2nd argument. Any hints? Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e-mail: r...@lcfltd.com Least Cost Formulations, Ltd.URL: http://lcfltd.com/ 824 Timberlake Drive Tel: 757

Re: [R] How to convert a factor column into a numeric one?

Exactly! Thanks. At 12:49 AM 6/5/2011, Jorge Ivan Velez wrote: Dr. LaBudde, Perhaps as.numeric(as.character(x)) is what you are looking for. HTH, Jorge On Sun, Jun 5, 2011 at 12:31 AM, Robert A. LaBudde wrote: I have a data frame: head(df) Time Temp Conc ReplLog10 10 -20H

Re: [R] How to convert a factor column into a numeric one?

Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e-mail: r...@lcfltd.com 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] How to convert a factor column into a numeric one?

)) +} ) str(de) 'data.frame': 20 obs. of 3 variables: $ a: num 1 1 1 1 2 2 2 2 3 3 ... $ b: num 1 1 2 2 1 1 2 2 1 1 ... $ y: num 0.6396 1.467 1.8403 -0.0915 0.2711 ... HTH, Dennis On Sat, Jun 4, 2011 at 9:31 PM, Robert A. LaBudde r...@lcfltd.com wrote: I have a data frame: head(df

Re: [R] Finding solution set of system of linear equations.

PLEASE do read the posting 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: r...@lcfltd.com Least Cost Formulations, Ltd

Re: [R] problem with glm(family=binomial) when some levels have only 0 proportion values

://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting 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: r

Re: [R] [R-sig-ME] Fwd: Re: ANOVA and Pseudoreplication in R

/mailman/listinfo/r-help PLEASE do read the posting 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: r...@lcfltd.com Least

Re: [R] comparing proportions

. Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e-mail: r...@lcfltd.com 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] comparing proportions

do I really need a mixed model here? Thanks again John From: Bert Gunter gunter.ber...@gene.com To: Robert A LaBudde r...@lcfltd.com Cc: array chip arrayprof...@yahoo.com Sent: Thu, February 10, 2011 10:04:06 AM Subject: Re: [R] comparing proportions Robert: Yes, exactly. In an offlist email

Re: [R] comparing proportions

(regular linear model) did not. Also as you can from the data, prop.success has some value outside the range 0.15 to 0.85, so maybe regular linear model may not be appropriate? Thank you, John From: Robert A LaBudde r...@lcfltd.com To: array chip arrayprof...@yahoo.com Cc: Bert Gunter

Re: [R] goodness-of-fit test

://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting 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: r

Re: [R] How to do bootstrap for the complex sample design?

. The fact that you choose a subsample size of 1 means you won't be able to estimate within-region variances unless you make some serious assumptions (e.g., UrbanCommunity effect independent of Region effect). Robert A. LaBudde

Re: [R] When to use bootstrap confidence intervals?

@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. Robert A. LaBudde

Re: [R] Function to compute the multinomial beta function?

)*Gamma(n3)/Gamma(n1+n2+n3) beta3- function (n1, n2, n3) exp(lgamma(n1)+lgamma(n2)+lgamma(n3)-lgamma(n1+n2+n3)) beta3(5,3,8) [1] 1.850002e-07 Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e-mail: r...@lcfltd.com Least Cost

Re: [R] Comparing distributions

-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. Robert A. LaBudde, PhD

[R] Singularity in simple ANCOVA problem

Residual standard error: 0.1787 on 6 degrees of freedom Multiple R-squared: 0.7464, Adjusted R-squared: 0.535 F-statistic: 3.532 on 5 and 6 DF, p-value: 0.07811 Why the NA's here? Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e-mail: r

Re: [R] Singularity in simple ANCOVA problem

of the groups. Rich Doh! Thanks. Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e-mail: r...@lcfltd.com Least Cost Formulations, Ltd.URL: http://lcfltd.com/ 824 Timberlake Drive Tel: 757-467-0954 Virginia

Re: [R] Popularity of R, SAS, SPSS, Stata...

are in a ship nearby you can see it happen. snip Once thing both whales and R users have in common is that, when you sight one, you say R, Matey! Thar she blows! Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e-mail: r...@lcfltd.com Least Cost

Re: [R] logistic regression with 50 varaibales

. Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e-mail: r...@lcfltd.com 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] glm poisson function

/mailman/listinfo/r-help PLEASE do read the posting 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: r...@lcfltd.com Least

Re: [R] When normality test violate and sample size is large(n=400) can I use t-test

distribution-free tests should have near perfect power, so which test you use is not important. Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e-mail: r...@lcfltd.com Least Cost Formulations, Ltd.URL: http://lcfltd.com/ 824

Re: [R] Histogram Bin

/posting-guide.html and provide commented, minimal, self-contained, reproducible code. Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e-mail: r...@lcfltd.com Least Cost Formulations, Ltd.URL: http://lcfltd.com/ 824 Timberlake

Re: [R] Cubic B-spline, how to numerically integrate?

limits of integration. Can anybody help please? Code and figures here: https://docs.google.com/fileview?id=0B1cQ7z9xYFl2ZTZhMmMyMjAtYTA3Zi00N2QyLTkxNzMtOGYyMjdiOGU2ZWE4hl=en Thank you, Claudia Robert A. LaBudde, PhD, PAS, Dpl

Re: [R] P values

At 07:10 AM 5/7/2010, Duncan Murdoch wrote: Robert A LaBudde wrote: At 01:40 PM 5/6/2010, Joris Meys wrote: On Thu, May 6, 2010 at 6:09 PM, Greg Snow greg.s...@imail.org wrote: Because if you use the sample standard deviation then it is a t test not a z test. I'm doubting that seriously

Re: [R] P values

. In Asymptopia, the Z and T distribution are identical. And it is only in Utopia that any P-value less than 0.01 actually corresponds to reality. Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e-mail: r...@lcfltd.com Least Cost Formulations, Ltd

Re: [R] Random numbers with PDF of user-defined function

Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e-mail: r...@lcfltd.com 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] From THE R BOOK - Warning: In eval(expr, envir, enclos) : non-integer #successes in a binomial glm!

deciding which example code to run. Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e-mail: r...@lcfltd.com Least Cost Formulations, Ltd.URL: http://lcfltd.com/ 824 Timberlake Drive Tel: 757-467-0954 Virginia

Re: [R] Nonparametric generalization of ANOVA

A search on bluesky...@gmail.com shows the user is in Norfolk, VA, USA. At 01:26 PM 3/5/2010, John Sorkin wrote: The sad part of this interchanges is that Blue Sky does not seem to be amiable to suggestion. He, or she, has not taken note, or responded to the fact that a number of people

Re: [R] Help with simple bootstrap test

, but this gives me a huge figure which can't possibly be right. Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e-mail: r...@lcfltd.com Least Cost Formulations, Ltd.URL: http://lcfltd.com/ 824 Timberlake Drive

Re: [R] multiple-test correlation

. Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e-mail: r...@lcfltd.com Least Cost Formulations, Ltd.URL: http://lcfltd.com/ 824 Timberlake Drive Tel: 757-467-0954 Virginia Beach, VA 23464-3239Fax: 757

Re: [R] random effects in mixed model not that 'random'

read the posting 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: r...@lcfltd.com Least Cost Formulations, Ltd

Re: [R] Implementation of gamma function for large number

the obvious question of what possible use gamma(17000) could be to you, and why a simple infinity would not work just as well. I'm sure you must have good reason for your request. Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e-mail: r

[R] Confusion on use of permTS() in 'perm'

to correspond to 3 patterns, not 4. I am misunderstanding how to solve this simple problem, or is something going on with permTS() that I'm missing. Thanks. Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e-mail: r...@lcfltd.com Least

Re: [R] Problem accessing functions in package 'roxygen'

:]]+$, , x) instead. At 01:42 PM 9/3/2009, Uwe Ligges wrote: Robert A. LaBudde wrote: I have Vista Home with R-2.9.0, and installed and tried to test the package 'roxygen': utils:::menuInstallPkgs() trying URL 'http://lib.stat.cmu.edu/R/CRAN/bin/windows/contrib/2.9/roxygen_0.1.zip

[R] Problem accessing functions in package 'roxygen'

? Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e-mail: r...@lcfltd.com 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

Re: [R] online classes or online eduction in statistics? esp. time series analysis and cointegration?

read the posting 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: r...@lcfltd.com Least Cost Formulations, Ltd

Re: [R] transformation/link glm in R

/listinfo/r-help PLEASE do read the posting 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: r...@lcfltd.com Least Cost

Re: [R] bootstrapped correlation confint lower than -1 ?

. Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e-mail: r...@lcfltd.com 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] How to test frequency independence (in a 2 by 2 table) with many missing values

. Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e-mail: r...@lcfltd.com Least Cost Formulations, Ltd.URL: http://lcfltd.com/ 824 Timberlake Drive Tel: 757-467-0954 Virginia Beach, VA 23464-3239Fax: 757-467-2947 Vere

Re: [R] How to test frequency independence (in a 2 by 2 table) withmany missing values

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: r...@lcfltd.com Least Cost Formulations, Ltd.URL: http

[R] Getting identify() to work with lattice::cloud()

is the correct syntax. Thanks. Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e-mail: r...@lcfltd.com Least Cost Formulations, Ltd.URL: http://lcfltd.com/ 824 Timberlake Drive Tel: 757-467-0954 Virginia Beach

Re: [R] Comparing model fits for NLME when models are not nested

? Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e-mail: r...@lcfltd.com 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] time series, longitudinal data or trajectories

measured at the same times for the subjects, the term panel data is used. For controlled experiments, the terms repeated measures and time series are common. Longitudinal could be used, but generally is not. Robert A. LaBudde

Re: [R] time series, longitudinal data or trajectories

, as it is a term from kinematics. I hope this helps. Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e-mail: r...@lcfltd.com Least Cost Formulations, Ltd.URL: http://lcfltd.com/ 824 Timberlake Drive Tel: 757-467-0954

Re: [R] Block factor as random or fixed effect?

. Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e-mail: r...@lcfltd.com 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] p-values from bootstrap - what am I not understanding?

, minimal, self-contained, reproducible code. Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e-mail: r...@lcfltd.com Least Cost Formulations, Ltd.URL: http://lcfltd.com/ 824 Timberlake Drive Tel: 757-467-0954

Re: [R] Physical Units in Calculations

. Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e-mail: r...@lcfltd.com 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

Re: [R] Genstat into R - Randomisation test

. Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e-mail: r...@lcfltd.com 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] course in ecological data analysis

help me I will be most grateful. Best regards - J. Capelle Try http://www.statistics.com/ Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e-mail: r...@lcfltd.com Least Cost Formulations, Ltd.URL: http://lcfltd.com/ 824

Re: [R] 'stretching' a binomial variable

advice much appreciated! nicola Work with events instead of proportions, and use a Poisson model. Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e-mail: r...@lcfltd.com Least Cost Formulations, Ltd.URL: http://lcfltd.com

Re: [R] Logistic regression problem

, reproducible code. 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] Located Latent Class Analysis (Uebersax)

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] FW: logistic regression

://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting 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

Re: [R] R function which finds confidence interval for binomial variance

looking for a reliable method for all p and for reasonable n. The proportion-based method is the best I've found, so far. Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e-mail: [EMAIL PROTECTED] Least Cost Formulations, Ltd.URL

[R] Confidence interval for binomial variance

binomVarCI.r #copyright 2008 by Robert A LaBudde, all rights reserved #CI for binomial sample variance #created: 09.26.08 by r.a. labudde #changes: require('binGroup') binomVarCI- function (n, x, conf=0.95) { p- x/n #proportion if (p0.25 | p0.75 | x==0 | x==n) { #use proportion-based CI pCI

[R] R function which finds confidence interval for binomial variance

. 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] Averaging 'blocks' of data

the posting 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

Re: [R] Projecting Survival Curve into the Future

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

Re: [R] Upgrading R means I lose my packages

will work. 5. I use the Packages|Update Package ... to update packages to 2.7.2. 6. Then I delete the 2.7.1 subfolder. You need Administrator rights to do this. Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e-mail: [EMAIL PROTECTED

[R] Specifying random effects distribution in glmer()

with Error in pgamma(id, shap, scal) : object shap not found. My questions are: 1. Can this be done? 2. Am I using the right package and function? 3. What am I doing wrong? Any help would be appreciated. Thanks. Robert A. LaBudde, PhD

Re: [R] Null and Alternate hypothesis for Significance test

. 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] shapiro wilk normality test

to the asymptotic region. Use sample() to select the data at random from within your data set to avoid bias in using the test. E.g., shapiro.test(sample(mydata, 1000, replace=TRUE)) Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e-mail: [EMAIL

Re: [R] shapiro wilk normality test

distributions, you will have to specify such a null hypothesis and test for deviation from it. E.g., testing for coefficient of skewness = 0. Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e-mail: [EMAIL PROTECTED] Least Cost Formulations

[R] Is there a good package for multiple imputation of missing values in R?

frame with the imputed values replaced (I don't have Herrell's book). Any pointers would be appreciated. Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e-mail: [EMAIL PROTECTED] Least Cost Formulations, Ltd.URL: http

Re: [R] Is there a good package for multiple imputation of missing values in R?

At 03:02 AM 6/30/2008, Robert A. LaBudde wrote: I'm looking for a package that has a start-of-the-art method of imputation of missing values in a data frame with both continuous and factor columns. I've found transcan() in 'Hmisc', which appears to be possibly suited to my needs, but I

Re: [R] How to make R running faster

. 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] Solving 100th order equation

. Is this a homework problem? 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] Simulation study in R

to sample X in the first place. Is its distribution uniform, or something else? Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e-mail: [EMAIL PROTECTED] Least Cost Formulations, Ltd.URL: http://lcfltd.com/ 824 Timberlake Drive

Re: [R] quasi-random sequences

list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting 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

Re: [R] bootstrap for confidence intervals of the mean

-contained, reproducible code. 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

[R] Multidimensional contingency tables

, as it has no means of using 'count'. This must be a common problem. Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e-mail: [EMAIL PROTECTED] Least Cost Formulations, Ltd.URL: http://lcfltd.com/ 824 Timberlake Drive

Re: [R] Multidimensional contingency tables

Now that is simple and elegant. Thanks! PS. Is there a course available for learning how to read R help information? :) At 10:52 PM 4/21/2008, Gabor Grothendieck wrote: xtabs(count ~., prob1) On Mon, Apr 21, 2008 at 10:46 PM, Robert A. LaBudde [EMAIL PROTECTED] wrote: How does one ideally

[R] Problem installing and using package tseries

is installed like any other, and the required file is at C:\Program Files\R\R-2.6.2\library\tseries\libs\tseries.dll. It is 36,864 bytes and looks superficially okay. Any ideas as to what's wrong here? Thanks. Robert A. LaBudde, PhD, PAS

Re: [R] Problem installing and using package tseries

hours. Best wishes, Uwe Ligges Robert A. LaBudde wrote: I have R 2.6.2, and have tried downloading and installing the package tseries. I get the same error when I use two different mirror sites: utils:::menuInstallPkgs() trying URL 'http://cran.mirrors.hoobly.com/bin/windows/contrib/2.6

Re: [R] generate random numbers subject to constraints

falls within the prescribed area. 3. Accept the point if it does. 4. Repeat if it doesn't. Efficiency depends upon the ratio of volumes involved. Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e-mail: [EMAIL PROTECTED] Least Cost Formulations

Re: [R] peak finding

off peaks (and symmetrically modified, valleys): #findpeaks() #Copyright 2007 by Robert A LaBudde, all rights reserved #find peaks in (x,y) curve findpeaks- function(x, y) { nx- length(x) ny- length(y) if (nx != ny) { print (' findpeaks01: x and y must be same length!') stop

Re: [R] stepAIC and polynomial terms

to a pseudo-quadratic relationship. You should be able to fit it better with the split plus a linear term. 5. Think about how age should affect your response before trying models. Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e-mail: [EMAIL

Re: [R] How to do multi-factor stratified sampling in R

At 03:54 PM 3/8/2008, David Winsemius wrote: Robert A. LaBudde [EMAIL PROTECTED] wrote in news:[EMAIL PROTECTED]: Given a set of data with a number of variables plus a response, I'd like to obtain a randomized subset of the rows such that the marginal proportions of each variable

[R] [R} Getting 'tilting' confidence intervals in R

[Resend with proper subject line] I am trying to compute bootstrap confidence intervals for a sample using R 2.5.1 for Windows. I can get Normal, Basic, Percentile, BCa and ABC from boot.ci() and boot() in the Davison Hinkley boot package. But I can't figure out how to use tilt.boot() to get

Re: [R] applying math/stat functions to rows in data frame

separately to work on a row? Try using t() to transpose the matrix, and then apply the column function of interest. Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e-mail: [EMAIL PROTECTED] Least Cost Formulations, Ltd.URL: http