search for: trubisz

Hi... In Stata, there is the ability to display scatter plots with data points at the same (x,y) location, using the 'jitter' command of the twoway scatter stata command. Anyone know of a way that I can do the equivalent thing in R? For non-Stata readers, if jitter is enabled in stata, and n-points would be at the same (x,y) location, the points are offset a bit (according to

I have a table statement that returns the following: [10.839,10.841] (10.841,10.843] (10.843,10.846] (10.846,10.848] (10.848,10.85] 0 0 0 0 1 (10.85,10.852] (10.852,10.854] (10.854,10.857] (10.857,10.859] (10.859,10.861] 0 0 0 0 0 What I want to do is get the upper bound

OK...newbie question here. Either I'm reading the docs wrong, or I'm totally confused. Given the following: x<-c("aaa","bbb","ccc") y<-rep(0,3) z<-rep(0,3) is.character(x) [1] TRUE is.numeric(y) [1] TRUE Now...I want to create a data frame, but keep the data types. In reading the docs, I assume you do it this way:

Hi... Is this possible in R? I have 2-sets of data, that were collected simultaneously using 2- different data acquisition schemes. The x-values are the same for both. The y-values have different ranges (16.4-37.5 using one method, 557-634 using another). In theory, if you plot both plots on top of each other, the graphs should overlap. The problem I'm having is trying to have to

Hi... I'm working through the book, A Handbook of Statistical Analyses using R by Everitt, and I'm trying to do the following (p. 19 of his book): boxplot(log(marketvalue)~country, data = subset(Forbes2000, country %in% c("United Kingdom","Germany","India","Turkey")), ylab="log(marketvalue",

Hi... I'm trying to understand the following syntax: cor.test(~mortality + hardness,data=water,method="pearson") which is the same as: cor.test(water$mortality,water$hardness,data=water,method="pearson") Can anyone point me to the correct doc or explain to me how to interpret "~mortality + hardness"? Thanks, Joe

Hi... I have a rather large dataframe that I'm trying to remove rows from. I'm issuing the command: dtx[-which(dtx$rdate > "2008-06-16"),] and it tries to print out over 170,000 lines of output. So...I did: options(max.print=1e6) and ran it again. It worked, but when I did a: which(dtx$rdate > "2008-06-16") it tells me that all the data I thought I deleted

Hello all... I'm attempting to write my own GAM plot function, so I can overlay it on top of an already existing plot. Problem is that after I do the gam, e.g. m<-gam(...), I cannot match the graph that gam.plot outputs when I attempt to plot the values from m$residuals, m$linear.predictors or m$fitted.values. Kind of at a loss what variables to use or if I need to do something

Hi... I am issuing a: plot(x,y), where x is a factor. i expect a box-and-whisker plot and I do get it. The data is 'bucketed' into 10 buckets on the x-axis. When I process some data, I get data in buckets 1,7,8,9 and 10. Problem is when I plot it, it shows only those buckets on the x-axis. if I subsequently do a points(4,3,pch-16) for example, it places the point not at x=4, but at x=9.