similar to: Shade area under curve?

Hi! With the following script, I'm trying to make a demonstration of a Confidence Interval, but I'm observing some differences on tails. # Teste de média entre uma amostra e uma população normal # Autor: Raphael de Freitas Saldanha # Agosto de 2008 n <- 200 # Sample size xbar <- 100 # Sample mean s <- 2 # Sample SD nc <- 0.95 # Confidence level (95%

Hi.. I'm not sure why polygon returns an area above the standard normal curve. z <- pretty(c(-3,3), 100) ht <- dnorm(z) data <- data.frame(z=z, ht=ht) zc <- 1.645 plot(data, type="l") lines(data) t <- subset(data, z>zc) polygon(t, col="red") Thanks, Lance [[alternative HTML version deleted]]

Dear all, is there an easy way to know at which position in the source code browser() (or stop()) was called? This seems to be important for me when I try to debug a program (e.g. using options(error=quote(dump.frames()))) traceback() returns some information, but I miss something like a line number, e.g. "stopped at line 999 in foo.R". Marcus --

Hi List, is sth. like the following possible: legend(x,y,c("A","B","C"), lty=list(1,"42","11")) ? Or: is there a possibility to define a solid line using the "string" mode for lty? Thanks Marcus -- +-------- ><> ------------------------------------------- | E-Mail: eger.m at gmx.de (NEW) | marcus.eger at

What is the best way (if possible) to pass the name of a workspace file to R when it starts so it is loaded instead of .RData ? I would like to know this for Linux and windows version. For windows I would like to associate the .RData extension so doubleclicking the file starts R and loads it. Thanks -- Victor Moreno V.Moreno at ico.scs.es Servei

I have tried to find a neat way to include graphs from R in LaTeX documents, but have not succeeded (I work with a WinEdt/MikTeX combination). The two roads I have stumbled along so far are the following: *Generate postscript files and convert them into EPS files by means of GhostScript or other not so straightforward tools. *Generate pictex files and include these. None of these solutions have

Dear R'ers, when I e.g. install the integrate2.2-2 package under R1.1.1 (debian2.2) atlas is linked into the library: Installing package `integrate' ... libs gcc -shared -lc -shared -o /usr/lib/R/library/integrate/libs/integrate.so adapt.o adapt_callback.o bsrl.o flfm.o nxprt.o symrl.o wht.o -L/usr/local/lib -lf77blas -lcblas -latlas -lg2c -lm -L/usr/lib/gcc-lib/i386-linux/2.95.2

Hello all, I'm using R1.2.2 and wonder whether there is already a new (unofficial) version that supports fill patterns for polygon(). thanks M.E. -- +------------------------------------------------------- | Marcus Eger +-------- ><> ------------------------------------------- | E-Mail: eger.m at gmx.de (NEW) | marcus.eger at physik.uni-marburg.de (OLD) | WWW:

Dear all, after building the R-1.1.1 from source on our alpha everything seems to work, but locfit still hangs... I installed locfit with install.packages(...), the installation process went fine, but calling locfit() blocks the R environment. Has anybody got an idea howto get to the roots of this problem? Thanks M. Eger -- +------------------------------------------------------- | Marcus

Dear Rs, quite often in I had segmentation faults of the following kind: Error in brf.infobase(nn, Pc, knams = knams, plot = plot, dev = dev, quiet = quiet, : unused argument(s) (df ...) Process R segmentation fault at Fri Sep 29 15:46:27 2000 This happened when I (wrongly) called a routine with an argument which the routine did not support. platform i386-pc-linux-gnu arch i386

Hi, many thanks for all suggestions as to shading of areas. Now to another small detail: how do I add items to the x-axis? Given a plot of "dnorm", I'd like to add some text to the x-axis, preferably with a corresponding tick mark. Imagine the typical \lambda demarcating a tail area. So far, I've only been able to add text in the plot region, not in the one of the axes.

Dear all, usually I work under linux, but I tried to install R-1.1.1 on our alpha's which worked well. However, installation of the locfit package produced the following warning: ---------%<---------------- Installing package `locfit' ... libs /vol/bin/gcc -I/vol/math/R/include -I/usr/local/include -mieee -g -O2 -c S_e nter.c -o S_enter.o S_enter.c: In function `basis':

Hello, how can I change (enlarge) the history buffer? I use R.64.1 and ESS5.1, emacs. Thanks - Marcus Eger -- +------------------------------------------------------- | Marcus Eger | Philipps Universitaet Marburg | Fachbereich Physik AG Neurophysik | Renthof 7 | 35032 Marburg | Tel. +49 6421 284169 | Fax +49 6421 287034 +------------------------------------------------------- | E-Mail:

Hello, is the following somehow possible with Mod() being a primitive: "Mod" <- function(x) UseMethod("Mod") "Mod.default" <- function(x) { .Primitive("Mod") ### I know it?s wrong. How can I access the original Mod ? } Thanks for help M. Eger -- +------------------------------------------------------- | Marcus Eger | Philipps Universitaet

Dear R-ers, somehow it is not possible to increase nsize to more than 20000k. When I specify e.g. > R --vsize=10M --nsize=21000K the result is: free total (Mb) Ncells 99658 350000 6.7 Vcells 1219173 1310720 10.0 Maybe I have overlooked s.th.... Marcus -- +------------------------------------------------------- | Marcus Eger | E-Mail: eger.m at gmx.de (NEW) |

Hello, I faced some difficulties using time series objects. ---------------------------------------------------------------------------- 1. Why is the time data changed, when I, e.g., apply sqrt(arrts), but not when arrts^0.5 ? This is my examples time series: ---------------------------------------------------------------------------- > arrts Time Series: Start = 3 End = 23 Frequency

Marcus Eger <marcus.eger@physik.uni-marburg.de> writes: > > time(sqrt(arrts)) > Time Series: > Start = c(1, 1) > End = c(5, 1) > Frequency = 1 > [1] 1 2 3 4 5 Looks like a bug... > ---------------------------------------------------------------------------- > 2. (At least) boolean indexing with matrices does not seem to work > properly: >

Hi. If I call the quantiles function as follows: qvec = quantiles(dvals,probs=seq(0,1,0.1)) the results will return a vector something like the following example: 0% 10% 20% 30% 40% 50% 60% 70% 80% 90% 100% 56.0 137.3 238.4 317.9 495.8 568.5 807.4 1207.7 1713.0 2951.1 8703.0 Now I want to assign the deciles, 1 - 10, to each

I've recently installed Debian/GNU Linux 2.2. The R package on the debian website and Official CD is 0.90.x, so I checked cran for a more recent deb package, with the following result. homebrew:~# apt-get install r-base Reading Package Lists... Done Building Dependency Tree... Done Some packages could not be installed. This may mean that you have requested an impossible situation or if you

Every now and then I realize that my attempts to document what all dataframes consist of are unsufficient. So far, I have been writing notes in an external file. Are there any better ways to do this within R? One possibility could be to set up the data as packages, but I would like to have a solution on a lower level, closer to data. I can't find any pointers in the standard manuals.