search for: labudd

...s problem in this range). The composite method gives reliable results over a wide range of rho. This should suit the purpose until someone comes up with a more theoretically sound method. Bootstrap is not reliable for n < 50, at least. #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 | x==0 | x==n) { #use proportion-based CI pCI<- binWilson(n,...

...h appears to be possibly suited to my needs, but I haven't been able to figure out 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: ral at lcfltd.com Least Cost Formulations, Ltd. URL: http://lcfltd.com/ 824 Timberlake Drive Tel: 757-467-0954 Virginia Beach, VA 23464-3239 Fax: 757-467-2947 "Vere scire est per causas scire"

...#39;Temp' in rows 5 and 7 are 45 and not 4, as expected, although the result is numeric. The same happens if I use the order given by levels(df$Temp) instead of the sort order in the recode() 2nd argument. Any hints? ================================================================ Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e-mail: ral at lcfltd.com Least Cost Formulations, Ltd. URL: http://lcfltd.com/ 824 Timberlake Drive Tel: 757-467-0954 Virginia Beach, VA 23464-3239 Fax: 757-467-2947 "Vere scire est per causas scire"

...yes 0 6 black white no 9 7 black black yes 6 8 black black no 97 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: ral at lcfltd.com Least Cost Formulations, Ltd. URL: http://lcfltd.com/ 824 Timberlake Drive Tel: 757-467-0954 Virginia Beach, VA 23464-3239 Fax: 757-467-2947 "Vere scire est per causas scire"

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: ral at lcfltd.com Least Cost Formulations, Ltd. URL: http://lcfltd.com/ 824 Timberlake Drive Tel: 757-467-0954 Virginia Beach, VA 23464-3239 Fax: 757-467-2947 "Vere scire est per causas scire"

...t work, because of the strong non-normal character of the sampling distribution for v (or v'). Does anyone know of an R package with R function that computes a reasonable confidence interval for v or v'? Thanks. ================================================================ Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e-mail: ral at lcfltd.com Least Cost Formulations, Ltd. URL: http://lcfltd.com/ 824 Timberlake Drive Tel: 757-467-0954 Virginia Beach, VA 23464-3239 Fax: 757-467-2947 "Vere scire est per causas scire"

...en' was built under R version 2.9.1 > trim(" 1234") Error: could not find function "trim" 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: ral at lcfltd.com Least Cost Formulations, Ltd. URL: http://lcfltd.com/ 824 Timberlake Drive Tel: 757-467-0954 Virginia Beach, VA 23464-3239 Fax: 757-467-2947 "Vere scire est per causas scire"

...kage 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, Dpl. ACAFS e-mail: ral at lcfltd.com Least Cost Formulations, Ltd. URL: http://lcfltd.com/ 824 Timberlake Drive Tel: 757-467-0954 Virginia Beach, VA 23464-3239 Fax: 757-467-2947 "Vere scire est per causas scire"

Hello all! I am once again analyzing patient survival data with chronic liver disease. The severity of the liver disease is given by a number which is continuously variable. I have referred to this number as "meld"--model for end stage liver disease--which is the result of a mathematical calculation on underlying laboratory values. So, for example, I can generate a Kaplan-Meier plot

...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, PAS, Dpl. ACAFS e-mail: ral at lcfltd.com Least Cost Formulations, Ltd. URL: http://lcfltd.com/ 824 Timberlake Drive Tel: 757-467-0954 Virginia Beach, VA 23464-3239 Fax: 757-467-2947 "Vere scire est per causas scire"

How do u calculated p values for a z test.. so far i ve done this A = read.table("cw3_data.txt") xbar = mean(A) s = 1 n = 20 mu = 0 z.test = (xbar-mu)/(s/sqrt(n)) p.value = pnorm(abs(z.test)) error = qnorm(0.99)*s/sqrt(n) left = xbar - error right = xbar + error and have got values off of it...but the values for p dont match up with other sites that i have used to check it

...or interaction 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, PAS, Dpl. ACAFS e-mail: ral at lcfltd.com Least Cost Formulations, Ltd. URL: http://lcfltd.com/ 824 Timberlake Drive Tel: 757-467-0954 Virginia Beach, VA 23464-3239 Fax: 757-467-2947 "Vere scire est per causas scire"

...the first arguments of identify(), but no matter what I use, when I click on a point I get warning: "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: ral at lcfltd.com Least Cost Formulations, Ltd. URL: http://lcfltd.com/ 824 Timberlake Drive Tel: 757-467-0954 Virginia Beach, VA 23464-3239 Fax: 757-467-2947 "Vere scire est per causas scire"

...d, family=Gamma(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 e-mail: ral at lcfltd.com Least Cost Formulations, Ltd. URL: http://lcfltd.com/ 824 Timberlake Drive Tel: 757-467-0954 Virginia Beach, VA 23464-3239 Fax: 757-467-2947 "Vere scire est per causas scire"

...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? ================================================================ Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e-mail: ral at lcfltd.com Least Cost Formulations, Ltd. URL: http://lcfltd.com/ 824 Timberlake Drive Tel: 757-467-0954 Virginia Beach, VA 23464-3239 Fax: 757-467-2947 "Vere scire est per causas scire"

...questions in the help systems. Obviously this can be done by writing new, brute force, functions scanning the subscripts, but there must be an (obvious?) direct way of doing this more elegantly. Thanks for any pointers. ================================================================ Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e-mail: ral at lcfltd.com Least Cost Formulations, Ltd. URL: http://lcfltd.com/ 824 Timberlake Drive Tel: 757-467-0954 Virginia Beach, VA 23464-3239 Fax: 757-467-2947 "Vere scire est per causas scire"

Hi, I have a dataset that has 2 groups of samples. For each sample, then response measured is the number of success (no.success) obatined with the number of trials (no.trials). So a porportion of success (prpop.success) can be computed as no.success/no.trials. Now the objective is to test if there is a statistical significant difference in the proportion of success between the 2 groups of

Hi I am using R version 2.4.1. How can I upgrade to version 2.5 without having to install all the packages again? Thanks Jason Dr. Iasonas Lamprianou Department of Education The University of Manchester Oxford Road, Manchester M13 9PL, UK Tel. 0044 161 275 3485 iasonas.lamprianou@manchester.ac.uk ___________________________________________________________ now. [[alternative HTML

Hi everybody, somehow i dont get the shapiro wilk test for normality. i just can?t find what the H0 is . i tried : shapiro.test(rnorm(5000)) Shapiro-Wilk normality test data: rnorm(5000) W = 0.9997, p-value = 0.6205 If normality is the H0, the test says it?s probably not normal, doesn ?t it ? 5000 is the biggest n allowed by the test... are there any other test ? ( i know qqnorm

...aking some type of trivial error here, but I don't know how to find it. I would guess that it has 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: ral at lcfltd.com Least Cost Formulations, Ltd. URL: http://lcfltd.com/ 824 Timberlake Drive Tel: 757-467-0954 Virginia Beach, VA 23464-3239 Fax: 757-467-2947 "Vere scire est per causas scire"