similar to: R-beta: Re: S Compatibility

At 03:28 30/04/97, ihaka at stat.auckland.ac.nz wrote: >Bill Venables writes: > (As a complete side-issue, Brian Ripley and I have a kind of > convention: we refer to the language as "S" and the commercial > product as "S-PLUS". There is a useful distinction to be made.) > >This is generally what I try to do too. > >However, I suspect though that most

Bill Venables writes: > Are the scoping differences between R and S set out precisely and > definitively somewhere? This would be useful. In the source code perhaps? :-) You can find a pretty precise description in the article Robert and I did in JCGS. Actually its pretty simple. Functions have access to the variables which were in effect when the function was defined. f <-

I have a small problem. I am running glm() in R-0.61.0 on Redhat 4.2. I want to get the standard errors from the output. If I do out <- glm(....) summary(out) I get the coefficients printed as well as their correlation matrix. If I do out$coefficients I get the coefficients out$fitted gives me the fitted values I can then assign the fitted values or the value of the estimated

Hello everybody, I am doing a plot and I would like to know which parameter I have to use in order to get an equal scale for the axes x and y: Means, if the intervals in the x-axis is 0,5,10,15.... I would like to have the same scale for the y-axis. Thank you so much for your help. Halima from Leiden where the sun is shining. Holland

Hello, I want to try the Design and Hmisc library from FE Harrell but, i have trouble with the copyright for the library Hmisc : in the home of the author : http://fharrell.biostat.virginia.edu/s/unix/ in the file Hmisc.README ###COPYRIGHT NOTICE ###You may not port code in the Hmisc library to R. Aie, Why ???? ###You may distribute these functions freely as long as you do so without

I want to do something like this in R. If I have three vectors > a1 [1] 1 2 3 > a2 [1] 4 5 6 > a3 [1] 9 10 7 I want to compute 1. A vector that is the mean at each serial position of a1, a2, and a3. so in this example it would have the contents 4.667, 5.667, 5.333333 2. A vector that is the SD at each serial position of a1, a2, and a3. so in this example it would have the contents

A non-text attachment was scrubbed... Name: not available Type: text Size: 1538 bytes Desc: not available Url : https://stat.ethz.ch/pipermail/r-help/attachments/19980220/040a76f5/attachment.pl

A non-text attachment was scrubbed... Name: not available Type: text Size: 1480 bytes Desc: not available Url : https://stat.ethz.ch/pipermail/r-help/attachments/19980326/974231c9/attachment.pl

Hi, I have just seen Stefan Steinhaus' web page : http://www.uni-franfurt.de/~stst/ncrunch.html I think it would be nice to include "R" as well. I have taken Forrest Young's email on stat-lisp list and changed the stuff for "R" :) Here it is: (someone please check this so we can also send it to Stefan Steinhaus.

It seems that the following should work: > x<-read.table('1A.tab',header=T,skip=1); > mode(x[1]) [1] "list" > x[1] p01 R1-00 0.0295 R1-01 0.0283 R1-02 0.0145 R1-03 0.0235 R1-04 0.0339 R1-05 0.0239 R1-06 0.0414 R1-07 0.0259 > hist(x[1]) Error: hist: x must be numeric > as.numeric(x[1]) Error: (list) object cannot be coerced to vector type 14 >

Dear all, Could you please help me how to get the output as I described in the following example? x<-c(543, 543, 543, 543, 551 , 551 ,1128 ,1197, 1197) diff<-x-lag(x) diff [1] NA 0 0 0 8 0 577 69 0 How to index the occasions in x repeatedly if the diff<15? if diff>=15, it will give a new index. I want the output be like y. y<-c(1,1,1,1,1,1,2,3,3) Thank you so

Dear all, Could you please help me how to get the output as I described in the following example? x<-c(543, 543, 543, 543, 551 , 551 ,1128 ,1197, 1197) diff<-x-lag(x) diff [1] NA 0 0 0 8 0 577 69 0 How to index the occasions in x repeatedly if the diff<15? if diff>=15, it will give a new index. I want the output be like y. y<-c(1,1,1,1,1,1,2,3,3) Thank you so

Hi All, I have a vector of data on which I am operating. The equation with which I want to operate on the vector has a value k. I want to run the equation and output a new vector, each time replacing k with each value from the sequence I defined. I have thought about using for loops and such, but this seems like overkill. I am wondering if there is a simple solution that would allow me to

Hi syslinux list, Some months ago, hpa was kind enough to help me troubleshoot booting Extlinux in a system with the following particulars: * Windows Vista x64 bootloader in the MBR and Vista installed in /dev/sda1 * Extlinux 3.72 installed in the first sector of /dev/sda2 and Foresight Linux installed in /dev/sda2 (at the time I tried, 3.72 was the default bootloader in Foresight) The problem

Hi, there, Suppose I have two vector say x=c(1 2 3 4 5) and y=(2 3 6 7). Then I want to combine these two vector together and get z=c(1 2 3 4 5 6 7) with 2 and 3 only appear once. I want to extend this one to a general case(say more than 100 elements in x and y and each time I don't know which elements are the same). Do you happen to know how to do this and which command should use? Thank

I want to use question numbers from my survey data (fake data below) as markers in a scatterplot. I'm using "as.character" to convert question numbers to characters. However, plot truncates the 0 (zero) off of question 10. How can I get the ending zero so I can add questions 11 - 20? # Sample code below Question <- c(1,2,3,4,5,6,7,8,9,10) Performance <-

Hi. I can quickly create a by object that selects rows from a data frame. After that, though, I don't know how to merge the rows back into a data frame that I can use. Here is an example where there is a data frame with three columns, a, b, and c. I update it so that there are two rows for each combination of a and b. I use by to select the subgroups of rows that share the same a and b

I thnk what you want is ?outer. e.g.: outer(Data -min(Data),value,FUN = "+") Whether this works for your real task, however, may depend on details and complexities that you have omitted. Cheers, Bert Bert Gunter "The trouble with having an open mind is that people keep coming along and sticking things into it." -- Opus (aka Berkeley Breathed in his "Bloom

Hello, I run R on Windows 2000. When I load lattice library and plot some lme class objects, the graphic window displays a dark-green color. Graphics look ok on the screen, but lines and points become invisible whem I send the plot to a printer because of the background color. Are there ways to reset the background color? Both par(bg=) and windows() do not work. Any ideas? Richard

Hello Everyone, I have three variables: Waist circumference (WC), serum triglyceride (TG) level and gender. Waist circumference and serum triglyceride is numeric and gender (male and female) is categorical. From these three variables, I want to calculate the "Lipid Accumulation Product (LAP) Index". The equation to calculate LAP is different for male and females. I am giving both