similar to: cumulative distribution functions

Hi, I have a question for the group, perhaps someone can help me figure this out. I've already looked in the help files and they were no help to me. I have a vector of values and I am plotting an ecdf graph. 1. How can i draw a continuous line through the ecdf points? (lines and type for the plot with an ecdf object does not work) 2. Supposing I have this line drawn. I can add a

Hello again, There is any package which does Z aware (real 3D) interpolations? It can be any method (IDW, kriging or spline) but it should take into consideration not only x and y coordinates for interpolation, but also z coordinate. I looked into different packages but it seems i didn't find the right one. The ultimate goal is to import the output into a GIS (Geographical Information

Hi, I am a beginner user of R. I have a trivial question ? I am almost ashamed I cannot figure it out does not matter how many times I am reading the help. I have a table in .txt format, tab delimited. I can read it with ?read.delim()? with no problems. Afterwards I would like to use boxplot function to see if there are any outliers in the column 5 of my data called TPAH16.ppm In the

Dear R experts, beginners and everyone else, I'm calculating "cumulative periodogram" using the command "cpgram" [1] from the MASS library. Here is a short example with the "lh" (hormone level) dataset: library(MASS) plot(lh,type="l",ylab="value",xlab="time", main="Hormone Levels (lh)") spectrum(lh,

Hi, I am trying to increase the memory R can use. I am running R under Windows on a machine with 2 GB of physical RAM and 4GB of paged memory. I wrote in the R property windows --sdi --max-mem-size=4094M, but the R itself when it is doing a bayesian modelling (geoR) it stops at 1,096K and i get memory errors because it cannot allocate a new segment of about 500K of memory. I don't have

Hi, I would like to write a little automatic routine in R, but i am a too much of a beginner for that. I will appreciate any help regarding this particular problem. Let?s suppose I have a data.frame with j columns (from 1 to n) and i rows (from 1 to p). I would like to write a procedure which reads every column j (j from 1 to n) and compare each value with the interval [0,1]. If z(i,j) is

Hi, Thank you so much for all your rapid answers. I am impressed. What i didn't know was that i have to assign my data to an object to work further on. It was not clear from the help (at least for me) that 'data()' itself is calling data already in R packages. All of you make that clear. Now, if you can suggest any good package to use for identifying outliers it will be great

Hi everybody, I would like to plot only a part of a data frame, and identify the IDs for all the points with values higher than a certain value. I will try to be more explicit: I did a qqnorm plot of my data. It is obvious from the plot that all the sample quantiles which are in Theoretical Quantiles = or > than 1 belong to another group. To be sure about that I plot the cooks.distance

Hi, Very so often when i am plotting something, doing a histogram, or whatever i am struggling to find out which are the numbers for different colors, palette names, types of lines, symbols, etc. Is there any documentation on line with all these numbers / names and the associated symbol / color??? For example if i am using the command image it uses a palette from red to yellow, with red the

Hi Again, First of all thank you for all the responses to my previous query. Your answers were very helpful and I did the job ;-). Now I hope you can answer as quick the following (sorry I am invading you with trivial questions): Let?s use again the following data.frame example: DF <- data.frame(x=rnorm(5), y=rnorm(5)) I want to obtain a new data.frame (or matrix) that contains only n

Hi, I have two variables ranging both from 0 to 1 (n=500 each). Now I am interested in plotting them both in one plot (using ggplot2). So far I used ecdf() (from an example I found with google) to get values for the cumulatice distribution function which gives a relative curve. I also want to do the same plot but using absolute cumulative values instead of relative. Can that be done with ecdf or

Hi, First of all - thank you for the answers regarding shared library. If i am attempting to translate the library from Linux to Windows - i think i will need some heavy baby-sitting ;-))) It is way out of my league - for now. Meanwhile i have this new question: suppose i have a data.frame with x and y columns and 10 rows, 1 to 10. I also have a variable m (or an array if you like) with 5

Hi y'all, First of all many thanks to Christian, Petr and Spencer for their replies. Your answers helped me to learn few more triks since then ;-) Now i have a question relating "aov". When i use aov i end up with 4 plots. How do i "see" the data behind those plots? I know about summary - but this gives me only some statistical info. Also, if i want to identify which

Hello, I have a data-set containing 22 variables, after appropriate transformations etc I ran a kmeans cluster analysis for 4 clusters , I ran it 20 times to find a result with the lowest within sum of squares. My question is how best do I go about finding out what the characteristics are of each cluster? Is one cluster dominated by a particular set of variables or by a particular

Hi everybody, I am using quite frequently the "robust" package and until now i never had any problems. Actually last time i used it was last Friday very successfully. Anyway, today anytime i want to use the function "fit.models" i get the following error even if i use the example form the help file: data(woodmod.dat) woodmod.fm <- fit.models(list(Robust =

Hi, I am doing a barplot, and the fist bar is very big (high values) and the rest of the bars are quite small (small values). So - is there any way to make the Y scale logarithmic so that i have a wider distance from 0 to 50 for example than for 50 to 100, and so on? Thanks in advance for any help, Monica

Hello, My specific system is > version _ platform i386-pc-mingw32 arch i386 os mingw32 system i386, mingw32 status major 1 minor 5.1 year 2002 month 06 day 17 language R Each of the hmisc and stepfun packages has a function

Hi, I am trying to do a bayesian modelling and my Windows PC seems not to have enough memory to do it. So now i would like to run GeoR on an IRIX 6.5 platform (UNIX related as i understood). Do you know if R is compatible with this type of platform? I am not allowed to try it until i have an answer to that. Thank you so much, Monica

Hi, I am trying to calculate mahalanobis distances for a matrix x with n*p variables. I am getting the following error: md2 <- mahalanobis(x, center, cov) Error in solve.default(cov, tol = tol.inv) : system is computationally singular: reciprocal condition number = 2.11165e-009 What does it means? Thank you so much for any help, Monica

Hi, This may be a very basic question but it seems i cannot figure it out does not matter what. how do you sort a table (ascending or descending) after the values in one particular column? I want to do something like the sort function in Xcel. Also, is there any other plot function that accepts log for y axes like the parplot2() from gregmisc? Thanks for the tip Marc about barplot2. very