similar to: plot - central limit theorem

Dear All I am totally new to R and I would like to know whether R is able and appropriate to illustrate to my students the Central Limit Theorem, using for instance 100 independent variables with uniform distribution and showing that their sum is a variable with an approximated normal distribution. Thanks in advance, Paul

Hello guys, I am stuck here: How do I make 1000 samples of n = 10 observations from an Exponential distribution and then compute the mean for all those 1000 samples? Basically I need to prove the Central Limit theorem, which states: http://www.nabble.com/file/p22664113/d175f06cbf200bd52a2c27a2e56dc594.png Where the Sn is sum of random variables, n we have from the question, mu is mean and

my data looks like this: PM10 Ref UZ JZ WT RH FT WR 1 10.973195 4.338874 nein Winter Dienstag ja nein West 2 6.381684 2.250446 nein Sommer Sonntag nein ja Süd 3 62.586512 66.304869 ja Sommer Sonntag nein nein Ost 4 5.590101 8.526152 ja Sommer Donnerstag nein nein Nord 5 30.925054 16.073091 nein Winter Sonntag nein nein Ost 6

In complex analysis, Cauchy's integral theorem states (loosely speaking) that the path integral of any entire differentiable function, around any closed curve, is zero. I would like to see this numerically, using R (and indeed I would like to use the residue theorem as well). Has anyone coded up path integration? -- Robin Hankin Uncertainty Analyst Southampton Oceanography Centre

I want to understand ANOVA better. But a few textbook that I have do not describe Cochran's Theorem in details. Could somebody recommend a book for me?

Dear Siegfried, I believe your boss is wrong saying that: >He also tried to explain me that the monthly means >(based on the daily measurements) must follow a >log-normal distribution too then over the course of a year. every statistician know that increasing the sample size the sample distribution of the mean is proxy to a gaussian distribution (Central Limit Theorem) independently

The code below is a small reproducible example of a much larger problem. While the script below works, it is really slow on the true dataset with many more rows and columns. I'm hoping to get the same result to examp, but with significant time savings. The example below is setting up a data.frame for an ensuing regression analysis. The purpose of the script below is to appends columns to

It is interesting to graph the distribution of the standardized average as n increases. Do this when the Xi are uniform on [0; 1]. Look at the histogram when n is 1, 5, 10 and 25. Do you see the normal curve taking shape? (A rule of thumb is that if the Xi are not too skewed, then n > 25 should make the average approximately normal. You might want > f=function(n,a=0,b=1) { mu=(b+a)/2

Hello to the R people, within ward clustering the distance calculated to decide the clustering of 2 subsets (h1 and h2) is the variation of inertia : d(h1,h2)=I(h1Uh2)-I(h1)-I(h2); i've been said that a way to calculate faster this d(h1,h2) is using the huygens theorem decomposing the inertia into "the inertia to the centroid + the distance to an axe" (that's my version ...). My

Hi, If I have following vector; x <- c(1,1,1,2,2,3,4,4,5) and I want to change values in the range of 1 to 3 into the value 1, how can I do that? I tried x[x == c(1:3)] <- c(1) but than I get; x [1] 1 1 1 2 1 1 4 4 5 R doesn't change the 2 into a 1. But why?

Hi, How can I make the distance between an axis-label and the axis bigger? I haven't found anything in par()...

Ein eingebundener Text mit undefiniertem Zeichensatz wurde abgetrennt. Name: nicht verf?gbar URL: <https://stat.ethz.ch/pipermail/r-help/attachments/20080929/e8ff7cb6/attachment.pl>

Hi, I tried to calculate the formula for the birthday problem (the probability that at least two people out of a group of n people share the same birthday) But the factorial-function allows me only to calculate factorials up to 170. So is there a way to push that limit? to solve this formula: (factorial(365) / factorial((365-23))) / (365^23) (n=23)

Hello, Today I updated one of our computers running CentOS 5.3 32bit. After a reboot I notice the following with lsof -i sshd 3638 root 3u IPv6 17317 TCP pc46.astro.up.pt:ssh->119.1.193.205:zephyr-clt (ESTABLISHED) sshd 3639 sshd 3u IPv6 17317 TCP pc46.astro.up.pt:ssh->119.1.193.205:zephyr-clt (ESTABLISHED) What is this connection?

R-lings, I know the algo for Laplace generation, but arent there any packages available. Also, i'm trying to check for the Bernoulli sums obey the CLT, but though i calculate 2000 values of S(n)-n*p ---------- sqrt(n*p*(1-p)) for p=.5 and n=45, this fails both KS test and shapiro test. Why?? Thanks sapsi ISI,Delhi please reply to sapsi at bigfoot.com

Hi, I have a data frame and would like to have summary statistics for grouped data. With summary() I get the central tendencies for the overall data. How can I get descriptive statistics with variances and standard deviations? for example my data.frame: group x y exp 2 4 exp 3 5 exp 2 4 control 1 2 control 2 3 control 1 2 now I want tables with summary statistics (variances

On Wed, Jun 12, 2013 at 03:56:36PM +0100, Ben Hutchings wrote: > On Mon, 2013-06-10 at 10:07 +0300, Michael S. Tsirkin wrote: > > On Thu, Feb 07, 2013 at 01:14:20PM -0500, David Miller wrote: > > > From: Ben Hutchings <bhutchings at solarflare.com> > > > Date: Thu, 7 Feb 2013 16:20:46 +0000 > > > > > > > If the consensus is still that we must

On Wed, Jun 12, 2013 at 03:56:36PM +0100, Ben Hutchings wrote: > On Mon, 2013-06-10 at 10:07 +0300, Michael S. Tsirkin wrote: > > On Thu, Feb 07, 2013 at 01:14:20PM -0500, David Miller wrote: > > > From: Ben Hutchings <bhutchings at solarflare.com> > > > Date: Thu, 7 Feb 2013 16:20:46 +0000 > > > > > > > If the consensus is still that we must

Hi, Is there a function which counts the frequencies of the occurence of a number within an interval? for example I have this vector: x <- c(1, 3, 1.2, 5, 5.9) and I want a vector that gives me the frequencies within an interval of 2, beginning at 0 (so the intervals are 0-2, 2-4, 4-6 and so on) so I get these frequencies: 2, 1, 2 Which functions do I have to use for this purpose?

From: Eric Dumazet <edumazet at google.com> Straightforward patch to add GRO processing to virtio_net. napi_complete_done() usage allows more aggressive aggregation, opted-in by setting /sys/class/net/xxx/gro_flush_timeout Tested: Setting /sys/class/net/xxx/gro_flush_timeout to 1000 nsec, Rick Jones reported following results. One VM of each on a pair of OpenStack compute nodes with