similar to: Excluding data with apply

Dear all, I've got many responses to my initial question, which is stated below. However, from those responses it has become clear that I need to rephrase my problem. All responses dealt with subscripting the data matrix before 'apply' is run on it. But this is not want I wanted to do. 'apply' cycles through rows or columns of a matrix, and runs a function on each row or

Dear R people I came across a strange thing: sd(rep(0.01, 15)) #0 sd(rep(0.001, 15)) #4.489023e-19 sd(rep(0.00001, 15)) #0 sd(rep(0.00000001,15)) #1.712427e-24 sd(rep(0.01, 13)) #1.805557e-18 sd(rep(0.001, 13)) #4.513894e-19 sd(rep(0.00001, 13)) #0 sd(rep(0.00000001,13)) #0 sd(rep(5.01, 15)) #0 sd(rep(5.001, 15)) #4.489023e-19 sd(rep(5.00001, 15))

Dear R-users I would like to draw a barplot with a special y-axis. Most of my data points are in the range from -50 to 50, apart from one value that is 550. I would now like to have a y-axis from -50 to 50, then a "break" and the "rest" of the y-axis, let's say from 500 to 600. In order to make clear what I'm trying to do, I attached a figure. Is it possible to do

Dear all, I am having exactly the same problem as Andy on an Intel Mac (see below / https://stat.ethz.ch/pipermail/r-sig-mac/2009-May/006208.html) with the degree symbol. I am using R version 2.11.1 (2010-05-31) on Xubuntu 10.04 with Windows fonts copied from /WINDOWS/Fonts and included via 'sudo dpkg-reconfigure fontconfig' and 'sudo fc-cache'. Any help would be highly

hello! i don't want to test my sample data for normality, but exponential- lognormal- or gammadistribution. as i've learnt the anderson-darling-test in R is only for normality and i am not supposed to use the kolmogorov-smirnov test of R for parameter estimates from sample data, is that true? can you help me, how to do this anyway! thank you very much! nadja

Hi I tried to customize help(), so it starts the compiled html help, and searches in all packages. I wrote the following function: > my.help function(topic){ help(topic, package=.packages(all=T), lib.loc=.lib.loc, chmhelp=T) } When I use it with >my.help(help) give the following error: Error in help(topic, package = .packages(all = T), lib.loc = .lib.loc, : No documentation

Dear all, I am interested in correctly testing effects of continuous environmental variables and ordered factors on bacterial abundance. Bacterial abundance is derived from counts and expressed as percentage. My problem is that the abundance data contain many zero values: Bacteria <-

Hello there I ` m a newbie in R and have a question about the output of write.table My command looks like this: write.table(datatr$V1[data.chull],file="outr.txt", quote = FALSE, sep = "\t") In the outputfile (outr.txt) I have an "x" in the left corner and a new row with numbers for each column(1,2,3...). I just want the values in my file. What should I change in

hello everybody i want to use the maximum likelihood method to estimate FRECHET parameters of my sample data. Should it work with fitdistr in the package MASS? I only find how to do it for GEV, Gumbel, and almost all other distributions, but FRECHET? I would be very happy if somebody can tell me how to do fit the FRECHET distribution! Thanks Nadja Riedwyl

I started the graphics device with X11(). When I tried to plot a tree object with > plot (tree) It gives me the error message: Error in plot.window(xlim, ylim, log, asp, ...) : invalid xlim I was looking in the manual under "plot", "plot.window", etc. but I didn't really get the resolution of my problem..... What's the problem? I'm runnig R on a SuSE 6.0

hello my data is data2:2743 4678 21427 6194 10286 1505 12811 2161 6853 2625 14542 694 11491 14924 28640 17097 2136 5308 3477 91301 11488 3860 64114 14334 by calculating shape<-(mean(data2))^2/var(data2) scale<-var(data2)/mean(data2) i get the idea what the parameters of the gammadistribution would be. but if i try using the method mle() i get stock and i don't know, how to

Dear R-users, I'm trying to extract the exact location of the maximum value in a vector or a matrix, meaning I want not only the maximum itself but also at which position (e.g row 6, column 12) it is located. I appreciate any suggestions. Thanks, Patrick -------------- next part -------------- A non-text attachment was scrubbed... Name: patrick.vcf Type: text/x-vcard Size: 339 bytes Desc:

Dear R Community, When making a summary with a vector, the result is "numeric" and I can recall its components (like median, mean, etc.) for further use. However, if I use summary with a matrix, the result is of mode "character" (like "Median : 1.2127") and I cannot extract the results for direct further use. Thanks for any hints, Patrick -------------- next part

Dear all, I am using Cairo version 1.4-5 for plotting R graphs on a Gentoo system that does not have X11 installed. I have full control on the fonts when running R in interactive mode. But when running R in batch mode I cannot change e.g. font size. Letters (including symbols) are, however, drawn. I am exporting the following environemntal variables when running R in batch mode to point to

Dear all, I've sent this question to the r-help mailing list, but maybe it is better suited for the r-sig-debian mailing list, although this is not Debian (but still Linux): I am using Cairo version 1.4-5 for plotting R graphs on a Gentoo system that does not have X11 installed. I have full control on the fonts when running R in interactive mode. But when running R in batch mode I cannot

Hi there Is there an easy function available in R for calculating the running mean? Reto -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the "body", not the subject !) To: r-help-request at

Dear all, I am stuck with a problem that might be trivial for most of you (and therefore is a bit embarrassing for me...): I want to calculate a generalized least squares regression using two time series (Y depending on X) with an autoregressive correlation structure of order two (the data along time are given below). I use 'gls' from package 'nlme': Calib.gls <- gls(Y~X,

Dear R Community, We have recently installed R 1.5.0 on several workstations. Ever since, if I try to extract certain numbers or lines from a file (which contains many NA's), I don't only get my desired line but all other lines which contain NA's. In a file with 15000 rows and 7500 NA's this is not very practical. In the earlier R-Version this worked when I typed in the

I've had a question a few moments ago about how to create multiple objects from multiple files within a loop. Thanks for the quick answers, it worked with "assign", like this: for(i in seq(1,nfn,1)){ fin<-paste("/home/klimet/patrick/LAEGEREN/NEBEL/FOGEVENT2000/",fn[i],sep="") assign(paste("f", i, sep = ""),

Dear all, I use R 2.0.1 on Windows XP professional. When I want to load the 'Gregmisc' library I get the following error message: Error in library(pkg, character.only = TRUE) : 'gregmisc' is not a valid package -- installed < 2.0.0? Can anybody tell me what's wrong with this package? Cheers, Christian