similar to: Applying glm coefficients (Beginner Question)

Hello, I'm working with a dataset that has 2 columns and 1000 entries. Column 1 has either value 0 or 1, column 2 has values between 0 and 10. I would like to count how often Column 1 has the value 1, while Column 2 has a value greater 5. This is my attempt, which works but doesn't seem to be very efficient, especially when testing different values or columns. count=0 for (i in 1:1000) {

something like ... > x<-c(1,2,3,0,5,6,0) > is.na(x[x==0])<-T > x [1] 1 2 3 NA 5 6 NA ~~~~~~~~~~~~~~~~ Robert Kinley ~~~~~~~~~~~~~~~~ [[alternative HTML version deleted]]

It's often useful to view objects in time order. In Splus I can do this with > objects.summary(order = "dataset.date") which delivers this sort of thing ... data.class storage.mode extent object.size dataset.date reference data.frame list 25 x 4 1700 2004.09.13 15:43 x data.frame list 15 x 4 1175

I really like RStudio ... ... but I wish it wouldn't automatically reload the last .RData it had. Anyone know how to fix this ... ? Also - does anyone know is there an Rstudio-user email-list forum thingy out there ? ta. Robert Kinley [[alternative HTML version deleted]]

platform: R 2.14.2 - Windows XP - Rstudio When I use cat() or print() in a script, the output text is of course mixed up with the lines of the script. I can sink() the output to a textfile and then play the file back at the end of the script, but that's not really what I'm after. It would be nice if the lines of text output could appear in a separate window, in the same

Hi I'm having difficulty working out how to get what I think is the appropriate partitioning of variability in a repeated measures setup. I have G=5 treatment-groups, each containing n=6 subjects, and a response is measured on each subject on t=4 occasions. I think the anova degrees of freedom should partition as follows - Between-subjects: G*n-1=29 [ between-groups: g-1 = 4 ,

hi I've used R for many years on windows machines, but have now acquired an Asus eee 1000 linux machine. In order to get the best out of the machine, I used the 'pimpmyeee.sh' script, to get the full KDE desktop. The version of Linux is Xandros, which I believe is a close relative of Debian, but sadly I have only a nodding acquaintance with Linux at present. Naturally I want to

Hi does anyone know of any R functions or packages for the following :- . generating cube plots for displaying the results from a response surface experiment design . generating ternary plots from a mixture design . tetrahedral plot or quaternary plot for displaying results from a four factor mixture experiment thanks Bob Kinley [[alternative HTML version deleted]]

Hi folks I have inherited a package which was created before release 2.10 and I need to have it working in release 2.12 There is a folder containing the Man, Data and R subfolders, and the description file, and there is also the corresponding zip file. I can install the package from the zip file, but when I try to load it I get the message :- Error: package was built before R

{ I think this message got rejected at the 1st attempt - trying again} R 2.15.1 , windows XP I have a very non-stationary bivariate time-series - say {xt,yt} t=1 ... lots. I want to do a bivariate density contour-plot of the whole series and then step through the series 1 second at a time plotting that second's {x,y} subset on top of the contour plot and losing the previous

Hi For bureaucratic reasons beyond my control I need to rewrite an R function (for producing operating characteristic curves) as an Splus function ( version 6 , windows XP ). The R function makes extensive use of the fact that the student's t distribution function pt() has a non-centrality parameter built in ... sadly that parameter is not present in the Splus pt() function .

An embedded and charset-unspecified text was scrubbed... Name: not available URL: <https://stat.ethz.ch/pipermail/r-help/attachments/20091126/5252ee65/attachment-0001.pl>

I think the rule is that you can do anything as long as you don't complain. If you want to complain, you must follow the instructions. -- Jari Oksanen in Re: [Rd] Keeping up to date with R-devel -- Patrick Burns pburns at pburns.seanet.com twitter: @burnsstat @portfolioprobe http://www.portfolioprobe.com/blog http://www.burns-stat.com (home of: 'Impatient R' 'The R

My small brain is having trouble getting to grips with lme() I wonder if anyone can help me correctly set the random = argument to lme() for this kind of setup with (I think) 9 variance/covariance components ... Study.1 Study.2 ... Study.10 Treatment.A: subject: 1 2 3 4 5 6 etc. 28 29 30 Treatment.B: subject: 31

I'm interested in using a data frame as if it were a hash table. For instance if I had the following, > (d <- data.frame(key=seq(0.5, 3, 0.5), value=rnorm(6))) key value 1 0.5 -1.118665122 2 1.0 0.465122921 3 1.5 -0.529239211 4 2.0 -0.147324638 5 2.5 -1.531503795 6 3.0 -0.002720434 Then I'd like to be able to quickly retrieve the "value" of "key" 1.5

Hello, I just started learning R and have a very basic question: When I try ar model and extract residuals it returns Null. Could someone explain what's going on here? Does that mean the model is null? But it seems residual has a length=# observation of the original series according to the model summary. I am learning it by myself and have no one to ask here so please allow me to ask this

Hi, I have an n x m matrix of numerical observations. ie. stock prices I wish to convert the matrix of observations to a matrix of simple returns (by taking the differences between (column) observations.) Can any good soul suggest a function for this? -- View this message in context: http://r.789695.n4.nabble.com/R-beginner-matrix-algebra-tp4653335.html Sent from the R help mailing list

Is there a reason that 'which.min' and 'which.max' don't have an 'arr.ind' argument? The context in which I wanted that was a grid search optimization, which seems like it would be reasonably common to me. -- Patrick Burns pburns at pburns.seanet.com http://www.burns-stat.com (home of 'Some hints for the R beginner' and 'The R Inferno')

I've just realized that it could be handy to have a 'decreasing' argument in 'is.unsorted'. And I'm cheekily hoping someone else will implement it. It is easy enough to work around (with 'rev'), but would be less hassle with an argument. The case I have in mind uses 'is.unsorted' in 'stopifnot'. Pat -- Patrick Burns pburns at pburns.seanet.com

R-help is all about solving R problems. So here ya go: http://www.portfolioprobe.com/2011/09/12/solve-your-r-problems/ -- Patrick Burns pburns at pburns.seanet.com twitter: @portfolioprobe http://www.portfolioprobe.com/blog http://www.burns-stat.com (home of 'Some hints for the R beginner' and 'The R Inferno')