search for: statisic

...ention I couldn't cut and paste yesterday because I was working R off a non-network computer while asking for help on a network computer. Essentially, I have this user-defined function: test <- function(X){ chisq <- svychisq(~X + SEX, design=audit, statisic="adjWald", round=4) } test(con) "con" is a data variable in my design object audit. When I just run: chisq <- svychisq(~X + SEX, design=audit, statisic="adjWald", round=4) It works just fine. It's only when it's nested in the function test() t...

I'm failing to understand how getrusage() works, which is a bit perplexing, because it doesn't seem like it would be terribly complicated. I've attached the code. My aim is to verify that I can use getrusage() to do (admittedly crude) instrumentation of which functions in my program are allocating lots of memory[1]. So I figure I can call getrusage() at various points and look at

> From: Robin Hankin > > [stuff about the CLT deleted] > > > > > So you can use R usefully to eveluate general statisical > > issues of this kind! > > > > absolutely! R is excellent for this sort of thing. I use it for > teaching stats all the time. > I'd say that without a tool like R you cannot learn statistics. > > > Consider an exponential distribution, which is very sk...

Hi wizards, I need to know how can I generate this numbers I have n numbers N1,N2 , Nn , but Nn = 1- sum(N1+N2+...+Nn-1) and sum(N1+N2+..Nn)=1 and N1,N2,..Nn with 0<N<1 . Does somebody know how to generate it ? Some distribution or algorithm. Thanks in advance. -- Web Page http://www.geocities.com/ricardo_rios_sv/index.html

I'd like to thank John Fox and Chuck Cleland for their help in resovling this issue. It turned out to be something simple, but perhaps others have had similar problems In my original data frame, I had 4 categories of race/ethnicity. One of the categories (other) was very small, and not similar to any of the other three categories, so I created a new data frame deleting those people.

Hi Everyone, I have a very long list of data-points (+2300) and i know from my histogram that there are outliers which are affecting my mean. I was wondering if anyone on here knows a way i can quickly get R to calculate and remove data which is 3 standard deviations from the mean? I am hoping this will tidy my data and give me a repeatable method of tidying for future data collection.

...2508^(1/3) > z [1] -2.205296 as you can see if Wq = -10.72508, Wq^(1/3) is NaN but z<- -10.72508^(1/3) returns a number. If someone can explain what I am doing incorrectly I would be most grateful. Below is the complete code for the function for reference: HallBoot <- function (x) { #statisics from the original data psd <- sqrt((sum((x-mean(x))^2)/length(x))) pmean <- mean(x) n <- length(x) k <- ((sum((x-mean(x))^3)/(length(x)*psd^3))) #the calulation of Q Qstat <- function (x,mean){ nb = length(x) bmean <- mean(x) Sb <- sqrt((sum((x-mean(x))^2)/nb)) Kb <- (...

Dear community, I've check it while working, but just to reassure myself. Let's say we have 2 models: model1 <- lm(vdep ~ log(v1) + v2 + v3 + I(v4^2) , data = mydata) model2 <- lm(vdep ~ log(v1) + v2 + v3 + v4^2, data = mydata) So in model1 you really square v4; and in model2, v4*^2 *doesn't do anything, does it? Model2 could be rewritten: model2b <- lm(vdep ~

Dear All, I have a problem understanding the difference between the outcome of a fisher exact test and a chi-square test (with simulated p.value). For some sample data (see below), fisher reports p=.02337. The normal chi-square test complains about "approximation may be incorrect", because there is a column with cells with very small values. I therefore tried the chi-square with

Dear R users, Can anyone please tell me how to generate a large number of samples in R, given certain distribution and size. For example, if I want to generate 1000 samples of size n=100, with a N(0,1) distribution, how should I proceed? (Since I dont want to do "rnorm(100,0,1)" in R for 1000 times) Thanks for help Debbie

...isting environment. as.environment() is now generic and uses list2env() as list method. ? There are several small changes to output which ?zap? small numbers, e.g. in printing quantiles of residuals in summaries from "lm" and "glm" fits, and in test statisics in print.anova(). ? Special names such as "dim", "names", etc, are now allowed as slot names of S4 classes, with "class" the only remaining exception. ? File .Renviron can have architecture-specific versions such as .Renviron.i386 on system...

...isting environment. as.environment() is now generic and uses list2env() as list method. ? There are several small changes to output which ?zap? small numbers, e.g. in printing quantiles of residuals in summaries from "lm" and "glm" fits, and in test statisics in print.anova(). ? Special names such as "dim", "names", etc, are now allowed as slot names of S4 classes, with "class" the only remaining exception. ? File .Renviron can have architecture-specific versions such as .Renviron.i386 on system...

...with x0 <- rep(1, length(y)) Message-ID: <58726081-253A-4326-9989-0D00037C3CCC@gmail.com> Content-Type: text/plain; charset=us-ascii Dear list, the following came up in an introductory class. Please help me understand the -1 (or 0+) syntax in formulae: Why do the enumerator dfs, F-statisics etc. differ between the models lm(y ~ x1) and lm(y ~ x0 + x1 - 1), if x0 is a vector containing simply ones? Example: N <- 40 x0 <- rep(1,N) x1 <- 1:N vare <- N/8 set.seed(4) e <- rnorm(N, 0, vare^2) X <- cbind(x0, x1) beta <- c(.4, 1) y <- X %*% beta + e summary(lm(y...