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. > &

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

ed from the original sample via jackknife as 1/6 the coefficient of skewness.) The script apparently works correctly if the number of resamples equals or exceeds the number of original data, but not otherwise. Robert A. LaBudd

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

argument 'eps' as a default 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, PAS,

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

uide.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 Dri

Re: [R] analysing a three level reponse

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

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

, 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 Virginia Bea

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

- within(dd, { + a <- as.numeric(as.character(a)) + b <- as.numeric(as.character(b)) +} ) > 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

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

his. See ?read.table for some of the options available. ============ 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 Vir

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 1

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

06547 22 -20H1 5.738683 37 -20H1 5.796394 4 14 -20H1 4.413691 50 45H1 6.406547 77 45H1 5.705433 but note that the values for 'Temp' in rows 5 and 7 are 45 and not 4, as expected, although the result is numeric. The

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

___ 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. =

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

___ 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. =

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

org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. __ 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 an

Re: [R] comparing proportions

value while model 1 (generalized mixed model) and model 3 (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: Rober

Re: [R] comparing proportions

ted in. So do I really need a mixed model here? Thanks again John From: Bert Gunter To: Robert A LaBudde Cc: array chip Sent: Thu, February 10, 2011 10:04:06 AM Subject: Re: [R] comparing proportions Robert: Yes, exactly. In an offlist email exchange, he clarified this for me, and I suggest

Re: [R] comparing proportions

n/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 Formulations, Ltd.

Re: [R] goodness-of-fit test

-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] How to do bootstrap for the complex sample design?

7;t be a problem. 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). ======== Robe

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

E 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.URL: http://lcfltd.com/ 824 Timb

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 L

Re: [R] Comparing distributions

__ 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. ===========

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

ship nearby you can see it happen. 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

Re: [R] Singularity in simple ANCOVA problem

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

[R] Singularity in simple ANCOVA problem

; 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 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? ===

Re: [R] logistic regression with 50 varaibales

ilman/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 Lea

Re: [R] glm poisson function

ng 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] When normality test violate and sample size is large(n=400) can I use t-test

e simplest 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://lcflt

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

quot; it. You need to give us the function and the appropriate limits of integration. Can anybody help please? Code and figures here: https://docs.google.com/fileview?id=0B1cQ7z9xYFl2ZTZhMmMyMjAtYTA3Zi00N2QyLTkxNzMtOGYyMjdiOGU2ZWE4&hl=en Thank you, Claudia ===========

Re: [R] Histogram Bin

thz.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...@lcfltd

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 wrote: Because if you use the sample standard deviation then it is a t test not a z test. I'm doubting that seriously... You calc

Re: [R] P values

a, 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

1.694695 ==== 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 pe

Re: [R] From THE R BOOK -> Warning: In eval(expr, envir, enclos) : non-integer #successes in a binomial glm!

read chapters 16 and 17 before 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

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 belie

Re: [R] Help with simple bootstrap test

mething to sum all the columns and rows of the table (or a subset should I desire). I've tried writing 'sum' for the function, but this gives me a huge figure which can't possibly be right. Robert A. LaBu

Re: [R] multiple-test correlation

) Then do the pairwise tests without familywise error control. I.e., this is similar to doing the F test in ANOVA before doing LSD testing. Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e-mail: r...@lcfltd.com Least Cost Formulations,

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

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

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

ain myself from asking 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. ============ Rober

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

e-tail p-value of 4/120 = 0.0333, or a 2-tail p-value of 2*4/120 = 0.067. This is not, however, what permTS() returns. The permTS() value of 0.05 appears 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

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

' I ended up using trim <- function(x) gsub("^[[:space:]]+|[[:space:]]+$", "", 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:::menu

[R] Problem accessing functions in package 'roxygen'

m" I have a similar problem with trim.right(), trim.left() and other functions I've tried. Any ideas? Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e-mail: r...@lcfltd.com Least Cost Formulations, Ltd.URL: http://lcfltd.com/ 824 Timber

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

nfo/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 Co

Re: [R] transformation/link glm in R

ling 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.

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

E 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] How to test frequency independence (in a 2 by 2 table) withmany missing values

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

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

osting-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 Drive

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

ning: "no point within 0.25 inches". Does identify() work with cloud()? If so, what is the correct syntax. Thanks. Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e-mail: r...@lcfltd.com Least Cost Formulations, Ltd.U

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

lly significant"? ============ 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-3239

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

used. If there are multiple responses all 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. =

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

e, as opposed to discrete times. "Trajectory" also implies an underlying causal model, as it is a term from kinematics. I hope this helps. Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e-mail: r...@lcfltd

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

not, stick with the fixed effects model. ============ 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-46

Re: [R] Physical Units in Calculations

already exists. ======== 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] p-values from bootstrap - what am I not understanding?

osting 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] Genstat into R - Randomisation test

ing-based approximate method. The name does create confusion in terminology for "randomization" tests for bootstrapping. Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e-mail: r...@lcfltd.com Least Cost Formulations, Ltd.URL: http://lcfltd.com/ 824 Tim

Re: [R] course in ecological data analysis

anyone can 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://lcflt

Re: [R] 'stretching' a binomial variable

- any 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://lcflt

[R] Unscaled factanal()

most instances, but occasionally I'd like to perform the analysis directly upon the covariance matrix itself. Thanks. ======== Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e-mail: r...@lcfltd.com Least Cost Formulations, Ltd.

Re: [R] Logistic regression problem

Tura, M.D,MPH,Ph.D National Institute of Cardiology Brazil __ 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 commente

Re: [R] Located Latent Class Analysis (Uebersax)

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

not take into account the stepwise >>>> procedure above that was used to torture the data until they >>>> confessed. >>>> >>>> Frank >>> Please book this as a fortune. >>> >>> Dieter >> Seconded! >> T

[R] Confidence interval for binomial variance

t. #09.26.08 02.50 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 (p<0.25 | p>0.75

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

ieve reasonable coverage. So I'm still 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 PROTE

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

#x27;? 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-3239Fax: 757-467-2947 "Vere

Re: [R] t-test between percentages

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

Re: [R] Averaging 'blocks' of data

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 PROTECTED] Least Cost Fo

Re: [R] Projecting Survival Curve into the Future

Richie. Mathematical Sciences Unit HSL ============ 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 Virg

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

point, the 2.7.2 has the same or new packages than 2.7.1, most or all of which 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. ====

[R] Specifying random effects distribution in glmer()

x27;t use nAGQ>1, not yet implemented summary(fit6) but this ends up 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 apprecia

Re: [R] Null and Alternate hypothesis for Significance 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. > [[alternative HTML version deleted]] _

Re: [R] shapiro wilk normality test

tribution. If you wish to test how the distribution is nonnormal, within some family of nonnormal distributions, you will have to specify such a null hypothesis and test for deviation from it. E.g., testing for coefficient of skewness = 0.

Re: [R] shapiro wilk normality test

ke you 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

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 suite

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

ut how to get a new data 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 Form

Re: [R] How to make R running faster

d minimize the number of lines of R. ============ 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-294

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 know how 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 Ti

Re: [R] quasi-random sequences

ex.ac.uk/atto > >__ >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. ===

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

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] 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: &

[R] Multidimensional contingency tables

7 black black yes 6 8 black blackno97 The xtabs() function doesn't seem appropriate, as it has no means of using 'count'. This must be a common problem. Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e-mail:

Re: [R] Problem installing and using package "tseries"

he problem and recompiled tseries which will appear on >CRAN master within 12 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

[R] Problem installing and using package "tseries"

;tseries' > There doesn't seem to be anything unusual. The package 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.

Re: [R] generate random numbers subject to constraints

the interior >of a general (k-r)-dimensional simplex embedded in >k dimensions, with (k+1) given vertices? > The method of "rejection": 1. Generate numbers randomly in the hypercube. 2. Test to see if the point falls within the prescribed area. 3. Accept the point if it does. 4.

Re: [R] generate random numbers subject to constraints

e random variates for X1, X2, based upon whatever unspecified distribution you wish. 2. Solve the two equations for X3 and X4. ============ Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e-mail: [EMAIL PROTECTED] Least Cost Formulations, Ltd.

Re: [R] peak finding

short function might be useful in picking 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 (&#x

Re: [R] stepAIC and polynomial terms

ge to see what the dependence is. Frequently curve is flat up to a certain age, and then linear thereafter. This gives rise 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 model

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

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

ction marginal proportions as well for some pairs. This must be a common problem in data mining, but I don't seem to be able to locate the proper library or function for doing this in R. Thanks for any help. ======== Robert A. LaBudde, PhD

[R] Getting 'tilting' confidence intervals in R

n the tilt.boot() call. This must be simple, but I'm not finding out the answer from the documentation. Also, use of any other packages to accomplish the tilting interval would also be useful. Thanks. Robert A. LaBudde, PhD,

[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.b

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

rest? Or, must these summary functions be programmed >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: [

Re: [R] Inferences for ratios of non-normal means

t;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.