similar to: y-axis with "break"

Hello I'm trying to calculate 2d confindence bounds into a scatterplot using the function "kde2d" (package MASS) and a contour plot. I found a similar post providing a solution - unfortunatly I do not realy understand which data I have to use to calculated the named "quantile": Post URL: http://tolstoy.newcastle.edu.au/R/help/03b/5384.html > (...) > >> Is

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 would like to announce the availability of a new package on CRAN: RNetCDF: R Interface to NetCDF Datasets This package provides an interface to Unidata's NetCDF library functions (version 3) and furthermore access to Unidata's udunits calendar conversions. The routines and the documentation follow the NetCDF and udunits C interface, so the corresponding manuals can be

Dear all, I would like to announce the availability of a new package on CRAN: RNetCDF: R Interface to NetCDF Datasets This package provides an interface to Unidata's NetCDF library functions (version 3) and furthermore access to Unidata's udunits calendar conversions. The routines and the documentation follow the NetCDF and udunits C interface, so the corresponding manuals can be

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

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 all I have a big matrix of standardized values (dimensions 285x5829) and R fails to calculate the principal components using prcomp() with the following error message: pc <- prcomp(my.matrix) Error in La.svd(x, nu, nv, method) : error code 1 from Lapack routine dgesdd Is the matrix too big? I'm using R-1.8.1 under Unix (Solaris8) and Linux(Suse 8.2). I tried to perform a principal

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

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 all, 'Apply' is a great thing for running functions on rows or columns of a matrix: X <- rnorm(20, mean = 0, sd = 1) dim(X) <- c(5,4) apply(X,2,sum) Is there a way to use apply for excluding rows or columns from a matrix to run functions on the remaining rows or columns? I know, I could do this with a 'for' loop, but 'apply' would be much easier and

Dear List, I''m thinking about using xen productive in our datacenter, i''m still testing around with it. Now I got some questions, just for basic understanding, we got for example this environment: 2 Nodes 1 SCSI Pool server (Connected via scsi to both nodes) Now I want to build a "cluster" so i would like to make this: Node 1 -> Primary -| | --> domU

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:

am I right in the assumption, that for calculation of the condition numbers I have to use the correlation matrix of X, and not t(x) %*% x? > e <- eigen(t(x) %*% x) better (x must not have a first column of ones): > e <- eigen(cor(x)) > e$val [1] 6.6653e+07 2.0907e+05 1.0536e+05 1.8040e+04 2.4557e+01 2.0151e+00 > sqrt(e$val[1]/e$val) [1] 1.000 17.855 25.153 60.785 1647.478

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

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

Hello everybody, I have a question. I would like to get a correlation between constitutive and induced plant defence which I messured on 30 plant species. So I have table with Species, Induced defence (ID), and constitutive defence (CD). Since Induced and constitutive defence are not independant (so called spurious correlation) I should do a randomisation test. I have a syntax of my

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 = ""),