similar to: R version of SAS/IntrNet

R-users, I have the following piece of code which I am trying to run on a dataframe (aga2) with about a half million records.  While the code works, it is extremely slow.  I've read some of the help archives indicating that I should allocate space to the p1 and ags1 vectors, which I have done, but this doesn't seem to improve speed much.  Would anyone be able to provide me with advice on

I'm using R0.90.1 on Win98. I'm trying to load an image I saved yesterday but I keep getting the following error: Error in load("C:/temp/245amData.R.RDa") : restore file corrupted -- no data loaded The image is 4,890KB. The image Can anyone tell me how to solve this problem? I'm sure I saved the image correctly. Thanks. Mark Mark J. Lamias Department of Statistics

I am trying to graph two planes on the same graph using persp(). I can only get one plane to plot at a time. Can someone explain how I can graph two planes on the same graph using persp? I've looked throught the documentation, but cannot find any references to appending a persp plot. Thanks. Sincerely yours, Mark J. Lamias Statistical Consultant

Does anyone have a piece of code or know how I can use R to generate a table of critical values for the Mann-Whitney (aka Wilcoxon Rank Sum) test. Ideally, I'd like a table that contains the critical values for any two samples of size 3 through 30. I could use Monte Carlo simulation or the normal approximation when n1 and n2 are greater than, 10, but I figured someone may know how to

Hi, I am trying to convert a colleague from using SPSS to R, but am having trouble generating a result that is similar enough to a bootstrapped cox regression analysis that was run in SPSS. I tried unsuccessfully with bootcens, but have had some success with the bootcov function in the rms package, which at least generates confidence intervals similar to what is observed in SPSS. However, the

Full_Name: Mark J. Lamias Version: 1.7.0 OS: Windows 2000 Pro Submission from: (NULL) (65.222.84.72) I am running the qwilcox procedure and it is producing incorrect results. For example, dwilcox(.025, 3, 5) should equal 6, but it is equal to 1. Similarly, dwilcox(.025, 3, 6) should equal 7, but it equals 2. The critical values are not set being returned with the correct values. I've

I would like to be able to submit a list of URLs of various webpages and extract the "content" i.e. not the mark-up of those pages. I can find plenty of examples in the XML library of extracting links from pages but I cannot seem to find a way to extract the text. Any help would be greatly appreciated - I will not know the structure of the URLs I would submit in advance. Any

I wish to create a matrix of all possible percentages with two decimal place percision. I then want each row to sum to 100%. I started with the code below with the intent to then subset the data based on the row sum. This works great for 2 or 3 columns, but if I try 4 or more columns the number of rows become to large. I would like to find a way to break it down into some kind of for loop, so

I am looking for advice on an editor to use with R (windows) that has an autosave feature. I typically write scripts using the RGui (and tried TinnR yesterday) but I am having continuing problems with BSODs (non R related) and have in the past have had issues with R crashes and would really like a system that does not require me to remember to hit the save button on my script every 10 minutes so

Is there an equivalent in R to the crosstabs command S-Plus? I've been searching for a while and can't seem to find anything. Thanks. Mark J. Lamias Department of Statistics Department of Political Science Survey Methodology Program/Survey Research Center Institute for Social Research - University of Michigan 426 Thompson Street, Room 315 Ann Arbor, Michigan 48104-2321 (734) 647-5381

On Tuesday, March 07, 2000 5:40 PM, Richard Bilonick wrote: >I need to merge several data.frames into one data.frame. In S-Plus I would >use >"merge" but I don't see a merge command in R. What is the best way to >accomplish >this? The easiest way to to this, I think, is as follows: if you have several data frames, dataframe1, dataframe2, . . . , dataframen, you

I have been trying to get Rgraphviz to work (I know it is from Bioconductor) unsuccessfully. Since I have no experience with Bioconductor I thought I would ask here if anyone has advice. I have installed Graphviz 2.20.3 as is recommended on the Bioconductor site but basically R cannot seem to find the needed dll files. So, even though I have added the appropriate directories to the system path R

I am trying to run a partialPlot with Random Forest (as I have done many times before). First I run my forest... Cell is a 6 level factor that is the dependent variable - all other variables are predictors, most of these are factors as well. predCell<-randomForest(x=tempdata[-match("Cell",names(tempdata))],y=tempdata$Cell,importance=T) Then I try my partial plot to look at the

We occasionally utilize the coxph function in the survival library to fit multinomial logit models. (The breslow method produces the same likelihood function as the multinomial logit). We then utilize the predict function to create summary results for various combinations of covariates. For example:

Has there been a recent change in the behavior of binary operators? In SparseM it was, until quite recently ok to do scalar multiplication but now, > A*4 Error in A * 4 : non-numeric argument to binary operator > 4*A Error in 4 * A : non-numeric argument to binary operator > A%*%A An object of class "matrix.csr" Slot "ra": [1] 1 1 1 Slot "ja": [1] 1 2 3

There are loads of resources for users of any other statistics package who are learning R. For example http://www.google.com/search?q="r+for+sas-users" The reverse isn't the case http://www.google.com/search?q="sas+for+r-users" Having heard rumors of how unelegant other statistics packages are, I don't see why any R-user would ever want to switch to anything else.

I am trying to create a dendrogram object to be able to display a multiway table. I understand how the dendrogram object is structured, a nested list of lists, but I am having trouble trying to create that structure from a table. Basically, what I would like to do is pass a multiway table, say table(a,b,c) to a function and have it display the table as a dendrogram. If a,b,c, and d are all

I am just starting to experiment with the party package and I am getting strange results. In the examples, the "statistic" and "criterion" seem related, i.e. criterion is a 1-p.value and statistic is the test statistic. Higher statistics are associated with higher criteria values. When I run these models on my own dataset, the highest statistic ends up getting a 0.00

Hello, I am starting to use R for various analyses, for example I use the ca package to do Correspondence Analysis. I am also looking to use packages such as: pls Partial Least Squares plspm Partial Least Squares Path Modelling However, although I can use packages such as these on un-weighted data there does not appear to be a facility to take account of weighted data. I am a

I am trying to run a very large Bradley-Terry model using the BradleyTerry2 package. (There are 288 "players" in the BT model). My problem is that I ran the model below successfully. WLMat is a win-loss matrix that is 288 by 288 WLdf<-countsToBinomial(WLMat) mod1<-BTm(cbind(win1,win2),player1,player2,~player,id="player",data=WLdf) Then I needed to run the same model