similar to: Loading Problem

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

Does an R package exist that is similar to SAS/Intrnet (http://support.sas.com/documentation/onlinedoc/intrnet/index.html) that will allow users to take parameters passed to it by a web request, produce results and return them to the web-browser in HTML format?  I tried looking through the package descriptions, but didnt' see anything like this.  Any direction or suggestions you could provide

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

Does anyone know a command to obtain the expected values of residuals under normality? Better yet, is there a command to perform a correlation test for normality (not Shapiro-Wilk Test)? I've looked everywhere, and can't seem to find how to do this. Thanks. Mark J. Lamias Department of Statistics Department of Political Science Survey Methodology Program/Survey Research Center

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 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

[These things don't seem to be of too great importance and may well not meet the criteria of a "bug", but, for what it's worth,I thought I'd submit them to r-bugs anyway:] This is the R version: > version _ platform i386-pc-mingw32 arch x86 os Win32 system x86, Win32 status major 1

Hi all! What is the equivalent R function of SPLUS crosstabs? Thank you. Danar. Stat. Inst. Umea Univ. Umea, Sweden -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- 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:

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

g'day R friends, I was hoping someone could provide me with advice as to my dysfunctional 'arrow' keys. Instead of browsing thru my command history they just move my cursor. What have I done wrong? My TERM seems ok for other unix apps. I'm running under hpux 10.20 with R1.0.0 thanks, John Strumila -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-

Dear friends. I'm in windows version R1.0.0. The code below produced a very beautiful and efficient graph in R Version 0.99.0 Patched (February 9, 2000) but in the newest, the x-axis is unreadable with dates written unformatted. I tried to set the numbers using par(xaxp) or xaxs but to no avail. Since it was OK just in the last version I guess there might perhaps be something wrong

Hi, I am not on this mailing list, but here is a terrible bug that has stopped me in my tracks. I am unable to remove observations from a data matrix. temp is the original matrix. Notice that there are 288 entries with a 104 in the first column. I attempt to remove these entries, but R does not do it. brad ACTUAL COMMANDS: > dim(temp) [1] 30528 11 > table(temp[,1]) 1 3 4

Hi, R allows me to run a one sample Wilcoxon test like this: wilcox.test(c(1,3.5,2.1,4,1.5,5), mu=2, exact=TRUE) The function 'wilcoxsign_test' from the package 'coin' should (I suppose) be able to calculate exact p values even if there are ties in the ranks. However, I couldn't find information on how to run a one sample test using 'wilcoxsign_test' like in the

Dear List, This problem is more a statistic one than a R one. Any one can recommend me some references website or online paper on maximum likelihood estimation?I'm now working on that,while still doubt how to prove that the estimated parameters are normal distributed. Thanks for your time and help! Best, Ning

Hi list, Could anyone recommend some good website for formatting R code? For example, when you copy paste R code to gmail and back to R, it loses its format, the dash symbol causes errors. I've had someone used it to format my code here on the list, but can't find it anymore. Mike [[alternative HTML version deleted]]

Hi all, I want to plot a series "x" and I want to put on title "Time evolution with alpha=0.2", e.g., plot(x, main="Time evotion with alpha=0.2",col="orange") Is there any posibility to put the greek letter and avoid written the word? Many thanks in advance. [[alternative HTML version deleted]]