Re: [R] glm: offset

2008-03-02 Thread Ted Harding
On 03-Mar-08 03:19:01, Wensui Liu wrote: > HI, John, > my understanding is that you should use log(...) instead of its > original scale. Below is the logic in the case of poisson reg. > log(y / offset) = x'b > => log(y) - log(offset) = x'b > => log(y) = x'b + log(offset) Well, this is where it get

[R] Studdy Missing Data, differentiate between a percent with in the valid answers and with in the different missing answers

2008-03-02 Thread Ericka Lundström
Hi R experts I'm trying to emigrate from SPSS to R, thou I have some problems whit getting R to distinguish between the different kind of missing. I want to distinguish between data that are missing because a respondent refused to answer and data that are missing because the question didn't

[R] Formating a zoo dataset .

2008-03-02 Thread Megh Dal
Suppose I have following dataset : > head(data1) Date Return 1 03/31/00 0.14230650 2 04/28/00 -0.03276228 3 05/31/00 -0.06527890 4 06/30/00 -0.04999873 5 07/31/00 -0.01447902 6 08/31/00 0.22265729 Now I convert it to zoo object : > data11 = zoo(data1[,2], as.Date(data1[,1],

[R] Attempting to connect to an Empress RDBMS via RODBC 1.2-3 causes R 2.6.2.pat to segfault

2008-03-02 Thread Michael Manning
Hi I've experienced some unpleasant behaviour while attempting to connect to an Empress 8.6.2 RDBMS via RODBC 1.2-3 and a freshly minted R 2.6.2.pat on a SLED linux system where a call to odbcConnect() to initiate the odbc connection causes R 2.6.2.pat to segfault. The odbcConnect call is CHA

Re: [R] R data Export to Excel

2008-03-02 Thread Hans-Peter
> > it doesn't seem indicate any installation instruction. > http://cran.r-project.org/web/packages/xlsReadWrite/index.html As Wensui Liu already mentioned, installation is the same as with other packages: In the RGui menu select Packages->Install Packages...->(Select a CRAN mirror)->(Select the

Re: [R] unable to start device PNG and unable to open connection to X11 display

2008-03-02 Thread Prof Brian Ripley
Please consult the help page for png. On Mon, 3 Mar 2008, Ng Stanley wrote: > Hi, > > I have installed R on a computational cluster, and am using putty to access > R. Please help on how to solve the problem of saving png files. > >> png(file="myplot.png", bg="transparent") > Error in X11(paste("p

Re: [R] Problem with PCA

2008-03-02 Thread Richard Rowe
This is a homework problem. I know how to do a PCA, you need to learn. I suggest you visit your textbook, then check the documentation for R's various PCA implementations to work out how to effect the analysis. phthao05 wrote: > I have an exercise. With 3 kinds of yohourt a,b,c. There are 25 >

[R] handling big data set in R

2008-03-02 Thread shu zhang
Hello R users, I'm wondering whether it is possible to manage big data set in R? I have a data set with 3 million rows and 3 columns (X,Y,Z), where X is the group id. For each X, I need to run 2 regression on the submatrix. I used the function "split": datamatrix<-read.csv("datas.csv", header=F,

[R] handling big data set in R

2008-03-02 Thread shu zhang
Hello R users, I'm wondering whether it is possible to manage big data set in R? I have a data set with 3 million rows and 3 columns (X,Y,Z), where X is the group id. For each X, I need to run 2 regression on the submatrix. I used the function "split": datamatrix<-read.csv("datas.csv", header=F,

[R] Constrained regression

2008-03-02 Thread Carlos Alzola
Dear list members, I am trying to get information on how to fit a linear regression with constrained parameters. Specifically, I have 8 predictors , their coeffiecients should all be non-negative and add up to 1. I understand it is a quadratic programming problem but I have no experience in the su

[R] Problem with PCA

2008-03-02 Thread phthao05
I have an exercise. With 3 kinds of yohourt a,b,c. There are 25 participatients estimate 3 norms: taste (va,vb,vc), structure (ca,cb,cc) and price (ga,gb,gc) and give the mark from 1 to 5. I don't know how to PCA this data. Please help me! I attached the data file follow: Va Vb Vc C

Re: [R] insert vector image into a plot

2008-03-02 Thread Paul Murrell
Hi William Morris wrote: > How can I insert a vector image in svg or pdf format into a plot. > Basically i need the equivalent of what the 'pixmap' package does for > bitmap images. 'grImport' See ... http://cran.r-project.org/web/packages/grImport/index.html http://www.stat.auckland.ac.nz

Re: [R] plotting to stdout (while reading from stdin?)

2008-03-02 Thread Gene Selkov
On Sun, 2 Mar 2008, Gabor Grothendieck wrote: > You can do this: > > Lines <- "A,B > 1,2 > 3,4" > DF <- read.csv(textConnection(Lines)) > > which is slightly simpler than the examples there. Thank you Gabor, I made it even simpler by replacing the Lines object with a string contstant. It wor

[R] R function to convert a number to text

2008-03-02 Thread lin tang
hi, Dear R users - I wonder is there a written R function to convert a number to a text, say convert 1 to "one" , 100 to "one hundred". I know in xls. has such a function BAHTTEXT, does anybody know is there a similar function in R ? Thanks. Lin -

Re: [R] R data Export to Excel

2008-03-02 Thread Wensui Liu
i think you simply install it in the way you install other R packages. On Sun, Mar 2, 2008 at 10:40 PM, Christophe Lo <[EMAIL PROTECTED]> wrote: > thanks for your response. How do i install it? I try looking at the manual > it doesn't seem indicate any installation instruction. I also download a >

Re: [R] R data Export to Excel

2008-03-02 Thread Christophe Lo
thanks for your response. How do i install it? I try looking at the manual it doesn't seem indicate any installation instruction. I also download a windows version but it doesn't have an exe file. http://cran.r-project.org/web/packages/xlsReadWrite/index.html Newbie, Kei On 3/3/08, Wensui Liu <

Re: [R] R data Export to Excel

2008-03-02 Thread Wensui Liu
hi, did you try write.xls in xlsReadWrite package? On Sun, Mar 2, 2008 at 9:59 PM, Keizer_71 <[EMAIL PROTECTED]> wrote: > > Here is my R Code > > x<-1:2 > y<-2:141 > data.matrix<-data.matrix(data[,y])#create data.matrix > variableprobe<-apply(data.matrix[x,],1,var) > variableprobe #outp

Re: [R] glm: offset

2008-03-02 Thread Simon Blomberg
Yes, use the log. I've had the same problem in the past, too. Try it on a toy example to confirm it for yourself. Cheers, Simon. On Sun, 2008-03-02 at 22:01 -0500, John Sorkin wrote: > R 2.6.0 > Windows XP > > A question about running a generalized linear model. > > I am running a glm with > (

Re: [R] glm: offset

2008-03-02 Thread Wensui Liu
HI, John, my understanding is that you should use log(...) instead of its original scale. Below is the logic in the case of poisson reg. log(y / offset) = x'b => log(y) - log(offset) = x'b => log(y) = x'b + log(offset) On Sun, Mar 2, 2008 at 10:01 PM, John Sorkin <[EMAIL PROTECTED]> wrote: > R 2.

[R] unable to start device PNG and unable to open connection to X11 display

2008-03-02 Thread Ng Stanley
Hi, I have installed R on a computational cluster, and am using putty to access R. Please help on how to solve the problem of saving png files. > png(file="myplot.png", bg="transparent") Error in X11(paste("png::", filename, sep = ""), width, height, pointsize, : unable to start device PN

Re: [R] R data Export to Excel

2008-03-02 Thread jim holtman
If you are asking how to convert to multiple columns in Excel, look at the "text to column" option in I think the data tab. On Sun, Mar 2, 2008 at 9:59 PM, Keizer_71 <[EMAIL PROTECTED]> wrote: > > Here is my R Code > > x<-1:2 > y<-2:141 > data.matrix<-data.matrix(data[,y])#create data.matrix >

[R] glm: offset

2008-03-02 Thread John Sorkin
R 2.6.0 Windows XP A question about running a generalized linear model. I am running a glm with (1) a poisson distribution and a log link: family=poisson(link = "log") and an offset. I would like to know if I should express the offset as the log of the offset value, i.e. offset=log(NumUniqPt)

[R] R data Export to Excel

2008-03-02 Thread Keizer_71
Here is my R Code x<-1:2 y<-2:141 data.matrix<-data.matrix(data[,y])#create data.matrix variableprobe<-apply(data.matrix[x,],1,var) variableprobe #output variance across probesets hist(variableprobe) #displaying histogram of variableprobe write.table(cbind(data[1], Variance=apply(data[,y],1,

Re: [R] insert vector image into a plot

2008-03-02 Thread William Morris
Thanks Paul. As always you spend an hour searching and find nothing, if i'd used the search term "import" instead of "insert" I would have found it in 2 secs. Lesson learnt, again. Cheers, Will Morris On 03/03/2008, at 12:50 PM, Paul Murrell wrote: > Hi > > > William Morris wrote: >> How ca

Re: [R] Need help to locate my mistake

2008-03-02 Thread Rolf Turner
On 3/03/2008, at 1:58 PM, Louise Hoffman wrote: > What if I need to calculate the variance for the fuel data? > > Are there a 'R' way to do that? > > I have derived > variance = (Y-x*theta)^T * (Y - x*theta) / (n-p) This should be (something like) tr((Y - X%*%theta))%*%(Y-X%*%theta)/(n-

Re: [R] emacs and R

2008-03-02 Thread Wensui Liu
Hi, John, you don't have to switch to linux in order to use ess + emacs with R. just follow the installation instruction of ess and it will take you 5 minutes at most. i also feel that xemacs seems more friendly than gnuemacs for windows user. On Sun, Mar 2, 2008 at 3:40 PM, John Sorkin <[EMAIL PR

[R] insert vector image into a plot

2008-03-02 Thread William Morris
How can I insert a vector image in svg or pdf format into a plot. Basically i need the equivalent of what the 'pixmap' package does for bitmap images. Cheers, Will Morris __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-

Re: [R] emacs and R

2008-03-02 Thread Michael Bibo
John Sorkin grecc.umaryland.edu> writes: > Can anyone suggest a document that I might read? I would suggest the ess manual itself: http://stat.ethz.ch/ESS/ess.pdf and John Fox's document: http://socserv.mcmaster.ca/jfox/Books/Companion/ESS/ess-xemacs.pdf. Even though the latter deals wit

Re: [R] About R-Project

2008-03-02 Thread Spencer Graves
I'm not familiar with "the Schuster system", but RSiteSearch("Schuster") produced 34 hits and RSiteSearch("river", "fun") produced 133 for me just now. If this does not lead you to an answer to your question, PLEASE do read the posting guide http://www.R-project.org/posting-guide.h

Re: [R] Recommended Packages

2008-03-02 Thread Wensui Liu
I am using the same get-em-all approach as yours. FYI. it is about 2G with all packages installed. On Sun, Mar 2, 2008 at 4:52 PM, Murray Jorgensen <[EMAIL PROTECTED]> wrote: > Having just update to R 2.6.2 on my old Windows laptop I notice that the > number of packages is growing exponentially

Re: [R] Need help to locate my mistake

2008-03-02 Thread Louise Hoffman
What if I need to calculate the variance for the fuel data? Are there a 'R' way to do that? I have derived variance = (Y-x*theta)^T * (Y - x*theta) / (n-p) __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read

[R] reducing RODBC odbcQuery memory use?

2008-03-02 Thread Thomas Pujol
1. Can I avoid having RODBC use so much memory (35 times the data size or more) making a data.frame & then .rda file via. sqlQuery/save? 2. If not, is there some more appropriate way from w/in R to pull large data sets (2-5GB) into .rda files from sql? [R] reducing RODBC odbcQuery memory use?

Re: [R] plotting to stdout (while reading from stdin?)

2008-03-02 Thread Gabor Grothendieck
You can do this: Lines <- "A,B 1,2 3,4" DF <- read.csv(textConnection(Lines)) which is slightly simpler than the examples there. On Sun, Mar 2, 2008 at 5:59 PM, Gene Selkov <[EMAIL PROTECTED]> wrote: > Hi Dirk, > > I didn't at first pay attention to your comment about littler, as my > original

Re: [R] Need help to locate my mistake

2008-03-02 Thread Louise Hoffman
> Unfortunately, x^(-1) is not the inverse of x: > > x<-matrix(c(2,4,4,5),nrow=2) > x > # [,1] [,2] > # [1,]24 > # [2,]45 > > x^(-1) > # [,1] [,2] > # [1,] 0.50 0.25 > # [2,] 0.25 0.20 > > i.e. it is the matrix which you get by applying the operation > (...)^(-

Re: [R] Need help to locate my mistake

2008-03-02 Thread Louise Hoffman
> This is certainly ***NOT*** correct. (If you really got those numbers > from Matlab, then Matlab is up to Puttee.) It was my mistake =) I had calculated the straight line so the edge of the plot was the y-axis =) > Have you plotted your data? > > (1) Fitting a straight line is ridicu

Re: [R] Recommended Packages

2008-03-02 Thread Gavin Simpson
On Mon, 2008-03-03 at 10:52 +1300, Murray Jorgensen wrote: > Having just update to R 2.6.2 on my old Windows laptop I notice that the > number of packages is growing exponentially and my usual approach of > get-em-all may not be viable much longer. Has any thought been given to > dividing "contr

Re: [R] plotting to stdout (while reading from stdin?)

2008-03-02 Thread Gene Selkov
Hi Dirk, I didn't at first pay attention to your comment about littler, as my original problem of plotting to stdout was solved. But it was just part of the larger problem: I actually need to be piping the data with the code for making the picture in, and getting the picture out without opening

Re: [R] Problem plotting curve on survival curve (something silly?)

2008-03-02 Thread Peter Dalgaard
Calum wrote: > OK this is bound to be something silly as I'm completely new to R - > having started using it yesterday. However I am already warming to its > lack of 'proper' GUI... I like being able to rerun a command by editing > one parameter easily... try and do that in a Excel Chart Wizzar

[R] Runtime error

2008-03-02 Thread Gwénaël Leday
Hi everybody! I try to implements Self-Organizing Maps with R and the Kohonen package. I have made an algorithm that try differents width and height. For each step, Silhouette index, Davies-Bouldin index and Dunn index are computed. However, there is always after severals minutes a *Runtime error*

[R] help.start() on linux (fedora 8) when firefox alreading running - a way to open a new tab?

2008-03-02 Thread Dale Steele
Using linux fedora 8 (x86_64) I get the following when firefox is already running. Is there a way to adjust settings in either R or firefox to open a new tab when help.start() is invoked? Thanks. --Dale > help.start() Making links in per-session dir ... If '/usr/bin/firefox' is already running,

Re: [R] discrete variable

2008-03-02 Thread David Winsemius
"Pete Dorothy" <[EMAIL PROTECTED]> wrote in news:[EMAIL PROTECTED]: > Hello, > > I am sorry for asking such a basic question. I could not find an > answer to it using google. > > I have a discrete variable (a vector x) taking for example the > following values : 0, 3, 4, 3, 15, 5, 6, 5 > >

[R] Problem plotting curve on survival curve (something silly?)

2008-03-02 Thread Calum
OK this is bound to be something silly as I'm completely new to R - having started using it yesterday. However I am already warming to its lack of 'proper' GUI... I like being able to rerun a command by editing one parameter easily... try and do that in a Excel Chart Wizzard! I eventually want

Re: [R] [OT] "normal" (as in "Guassian")

2008-03-02 Thread Ted Harding
Thanks, Katherine! Now I wonder what, in particular, Peirce might have had in mind (he was a particularly sharp philosophical thinker, and might be expected to pay attention to the "semantic baggage" of what he said). I'm also enjoying the other delightful "OT" (= "On Tangent") responses that my q

[R] Recommended Packages

2008-03-02 Thread Murray Jorgensen
Having just update to R 2.6.2 on my old Windows laptop I notice that the number of packages is growing exponentially and my usual approach of get-em-all may not be viable much longer. Has any thought been given to dividing "contributed" binaries into a recommended set, perhaps a couple of hundr

Re: [R] listing components of an object

2008-03-02 Thread Charilaos Skiadas
On Mar 2, 2008, at 4:12 PM, Nair, Murlidharan T wrote: > Is there a method to list the components of an object, instead of > looking at the help for that method? Let me be more clear with an > example > > data(iris) > ## tune `svm' for classification with RBF-kernel (default in svm), > ##

Re: [R] emacs and R

2008-03-02 Thread Peter Dalgaard
Louise Hoffman wrote: >> At the suggestion of many people, I have installed emacs on my linux (Fedora >> 8.0) computer with the intention of using emacs as window interface to R >> (2.6.0). I have gone though the emacs tutorial and don't see any information >> about how I should use emacs to run

Re: [R] emacs and R

2008-03-02 Thread Gavin Simpson
On Sun, 2008-03-02 at 21:45 +0100, Louise Hoffman wrote: > > At the suggestion of many people, I have installed emacs on my linux > (Fedora 8.0) computer with the intention of using emacs as window > interface to R (2.6.0). I have gone though the emacs tutorial and > don't see any information about

Re: [R] [OT] "normal" (as in "Guassian")

2008-03-02 Thread Katharine Mullen
There is some information and references regarding the name 'normal' in the internet article 'Earliest Known Uses of Some of the Words of Mathematics (N)', http://members.aol.com/jeff570/n.html, by John Aldrich. It contains the comment, "Galton does not explain why he uses the term "normal" but th

[R] listing components of an object

2008-03-02 Thread Nair, Murlidharan T
Is there a method to list the components of an object, instead of looking at the help for that method? Let me be more clear with an example data(iris) ## tune `svm' for classification with RBF-kernel (default in svm), ## using one split for training/validation set obj <- tune(svm, Species

Re: [R] Need help to locate my mistake

2008-03-02 Thread Ted Harding
On 02-Mar-08 20:18:29, Louise Hoffman wrote: > Dear readers > > I would like to make General Linear Model (GLM) for the > following data set http://louise.hoffman.googlepages.com/fuel.csv > > The code I have written is > > fuelData<-read.table('fuel.csv',header=TRUE, sep=',') > n<-dim(fuelData)[

Re: [R] Need help to locate my mistake

2008-03-02 Thread Rolf Turner
On 3/03/2008, at 9:18 AM, Louise Hoffman wrote: > Dear readers > > I would like to make General Linear Model (GLM) for the following > data set > http://louise.hoffman.googlepages.com/fuel.csv > > The code I have written is > > fuelData<-read.table('fuel.csv',header=TRUE, sep=',') > n<-dim(fuel

Re: [R] discrete variable

2008-03-02 Thread Andrew Robinson
Pete, try x <- c(0, 3, 4, 3, 15, 5, 6, 5) table(x) or length(table(x)) Cheers Andrew On Sun, Mar 02, 2008 at 09:27:20PM +0100, Pete Dorothy wrote: > Hello, > > I am sorry for asking such a basic question. I could not find an answer to > it using google. > > I have a discrete variable (a

Re: [R] [OT] "normal" (as in "Guassian")

2008-03-02 Thread Ken Knoblauch
Johannes Hüsing huesing.name> writes: > > Am 02.03.2008 um 17:44 schrieb Gabor Csardi: > > > I'm not a statistician, but do i remember well that among all > > distributions with a given mean and variance, the normal distribution > > has the highest entropy? This is good enough for me to call it

Re: [R] [PS] discrete variable

2008-03-02 Thread Ben Fairbank
Try table(), with the name of your vector inside the parentheses. Ben -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Pete Dorothy Sent: Sunday, March 02, 2008 2:27 PM To: r-help@r-project.org Subject: [PS] [R] discrete variable Hello, I am sorry for a

Re: [R] discrete variable

2008-03-02 Thread Ted Harding
On 02-Mar-08 20:27:20, Pete Dorothy wrote: > Hello, > > I am sorry for asking such a basic question. I could not find > an answer to it using google. > > I have a discrete variable (a vector x) taking for example the > following values : 0, 3, 4, 3, 15, 5, 6, 5 > > Is it possible to know how man

Re: [R] discrete variable

2008-03-02 Thread Henrique Dallazuanna
See > length(unique(x)) [1] 6 > table(x) On 02/03/2008, Pete Dorothy <[EMAIL PROTECTED]> wrote: > Hello, > > I am sorry for asking such a basic question. I could not find an answer to > it using google. > > I have a discrete variable (a vector x) taking for example the following > values : 0, 3,

Re: [R] emacs and R

2008-03-02 Thread Louise Hoffman
> At the suggestion of many people, I have installed emacs on my linux (Fedora > 8.0) computer with the intention of using emacs as window interface to R > (2.6.0). I have gone though the emacs tutorial and don't see any information > about how I should use emacs to run R. Can anyone suggest a d

[R] emacs and R

2008-03-02 Thread John Sorkin
At the suggestion of many people, I have installed emacs on my linux (Fedora 8.0) computer with the intention of using emacs as window interface to R (2.6.0). I have gone though the emacs tutorial and don't see any information about how I should use emacs to run R. Can anyone suggest a document

Re: [R] summarizing matrix data

2008-03-02 Thread Henrique Dallazuanna
On 02/03/2008, Emilio Gagliardi <[EMAIL PROTECTED]> wrote: > Henrique, > > Again, thank you very much for your code snippets, I am learning from them, > but they are also creating as many questions as answers... > > On Sun, Mar 2, 2008 at 11:48 AM, Henrique Dallazuanna <[EMAIL PROTECTED]> > wrote:

[R] discrete variable

2008-03-02 Thread Pete Dorothy
Hello, I am sorry for asking such a basic question. I could not find an answer to it using google. I have a discrete variable (a vector x) taking for example the following values : 0, 3, 4, 3, 15, 5, 6, 5 Is it possible to know how many different values (modalities) it takes ? Here it takes 6 di

[R] Need help to locate my mistake

2008-03-02 Thread Louise Hoffman
Dear readers I would like to make General Linear Model (GLM) for the following data set http://louise.hoffman.googlepages.com/fuel.csv The code I have written is fuelData<-read.table('fuel.csv',header=TRUE, sep=',') n<-dim(fuelData)[1] xOnes<- matrix(1,nrow=n,ncol=1) x<-cbind(xOnes,fuelData[,3])

Re: [R] Variance Calculation in R

2008-03-02 Thread Keizer_71
unfortunately, it is not showing probeID Henrique Dallazuanna wrote: > > Try this: > > write.table(cbind(data.matrix[1], Variance = apply(data.matrix[,-1], > 1, var)),file='file.xls') > > > On 02/03/2008, Keizer_71 <[EMAIL PROTECTED]> wrote: >> >> sorry...in step 4-i need the R code to

Re: [R] summarizing matrix data

2008-03-02 Thread Emilio Gagliardi
Henrique, Again, thank you very much for your code snippets, I am learning from them, but they are also creating as many questions as answers... On Sun, Mar 2, 2008 at 11:48 AM, Henrique Dallazuanna <[EMAIL PROTECTED]> wrote: > Then you can try this: > > apply(array(as.vector(sapply(patt, get)),

Re: [R] [OT] "normal" (as in "Guassian")

2008-03-02 Thread Johannes Hüsing
Am 02.03.2008 um 17:44 schrieb Gabor Csardi: > I'm not a statistician, but do i remember well that among all > distributions with a given mean and variance, the normal distribution > has the highest entropy? This is good enough for me to call it > "normal" There's more. Among all rotation-sym

[R] About R-Project

2008-03-02 Thread Ion Andronache
Hello, My name is Ion Andronache and I am a candidate for a doctor's degree on geomorphology in Braila, Romania. I would like to ask for some advice. I would like to make an analysis on the Danube's flow capacity and to determine the circles and periodicity using the Schuster system. I understood

Re: [R] Make plots with GNUplot. Have anyone tried that?

2008-03-02 Thread Paul Murrell
Hi Louise Hoffman wrote: > [snip] >> Seriously. Be specific if you have a problem. (read the posting guide). R >> can >> also plot. If you don't like R's plots (which I could not understand) you >> can >> export data and import them to gnuplot. So what? > > Okay, my post was not very good.

Re: [R] Variance Calculation in R

2008-03-02 Thread Henrique Dallazuanna
Then you can try: rownames(data.matrix) <- as.character(data.matrix$ProbeID) data.matrix <- data.matrix[-1] as.matrix(apply(data.matrix1, 1, var)) or out <- apply(data.matrix1, 1, var) data.frame(ProbeID = names(out), Variance = unname(out)) Works for me On 02/03/2008, Keizer_71 <[EMAIL PROT

[R] new to latex to pdf

2008-03-02 Thread Mike Babyak
Dear All, I'm trying to teach myself latex along with the latex function in Hmisc and have hit a roadblock that I can't seem to get around. I'd greatly appreciate any pointers. I'm running R 2.6.0 on Windows XP and have Miktex 2.7 installed. I've reproduced the code below, taken from Frank Ha

Re: [R] Imputation Packages

2008-03-02 Thread Dieter Menne
ArunPrasad uark.edu> writes: > I am looking for a package in R which can help me in using the > imputation technique to find the missing values for my regression analysis. The recommended search for imputation in r-project gave me 400 citations. I assume you have tried this. Diete

Re: [R] summarizing matrix data

2008-03-02 Thread Henrique Dallazuanna
Then you can try this: apply(array(as.vector(sapply(patt, get)), dim = c(dim(get(patt[1])), ngroups, nsubs)), 1:3, mean) For create the matrices in a loop one option is: lapply(seq(1, 51, by=10), function(x)matrix(seq(x, x+8), nrow=3, ncol=3)) On 02/03/2008, Emilio Gagliardi <[EMAIL PROTECTE

Re: [R] summarizing matrix data

2008-03-02 Thread Emilio Gagliardi
Hi Henrique, Thank you for your help, but it doesn't produce what I'm looking for... On Sun, Mar 2, 2008 at 10:50 AM, Henrique Dallazuanna <[EMAIL PROTECTED]> wrote: > One options is: > > data.1 <- matrix(seq(from=1,to=9,by=1),nrow=3,ncol=3) > data.2 <- matrix(seq(from=11,to=19,by=1),nrow=3,ncol

Re: [R] Variance Calculation in R

2008-03-02 Thread Keizer_71
Hi Henrique, It is definitely better, but it doesn't show me the ProbeID which identify the probes name Here was the result when i export to excel with your code. "Variance" 1 2.425509867 21.6216446425273 any suggestions? thanks, Kei Keizer_71 wrote: > > Hello, > > Thanks everyon

Re: [R] elegant way to minus on each row of a matrix

2008-03-02 Thread Bo Zhou
Hi Dimitris and everyone I tried this but now I know why it didn't work out for me initially. > t1=matrix(1:10,5,2) > t2=matrix(1,5,1) > t2 [,1] [1,]1 [2,]1 [3,]1 [4,]1 [5,]1 > t1 [,1] [,2] [1,]16 [2,]27 [3,]38 [4,]49 [5,]5 10 >

Re: [R] Variance Calculation in R

2008-03-02 Thread Henrique Dallazuanna
Try this: write.table(cbind(data.matrix[1], Variance = apply(data.matrix[,-1], 1, var)),file='file.xls') On 02/03/2008, Keizer_71 <[EMAIL PROTECTED]> wrote: > > sorry...in step 4-i need the R code to output in this format when i export > to > excel. > > ProbeID Variance > 1

Re: [R] elegant way to minus on each row of a matrix

2008-03-02 Thread Bo Zhou
Brilliant! > t1=matrix(1:15,5,3) > t1 [,1] [,2] [,3] [1,]16 11 [2,]27 12 [3,]38 13 [4,]49 14 [5,]5 10 15 > t2=1:5 > t2 [1] 1 2 3 4 5 > sweep(t1, 1, t2) [,1] [,2] [,3] [1,]05 10 [2,]05 10 [3,]05 10 [4,]0

Re: [R] elegant way to minus on each row of a matrix

2008-03-02 Thread Bo Zhou
Hi Mark, (CC'ing r-help) Only need to change the parameter from 1 to 2 after that it worked great. Thanks > t1=matrix(1:15,5,3) > t1 [,1] [,2] [,3] [1,]16 11 [2,]27 12 [3,]38 13 [4,]49 14 [5,]5 10 15 > t2=1:5 > t2 [1] 1 2 3 4 5 >

Re: [R] elegant way to minus on each row of a matrix

2008-03-02 Thread Dimitris Rizopoulos
try this: mat <- matrix(1:9, 3, 3, TRUE) dat <- as.data.frame(mat) vec <- 1:3 result.mat <- mat - vec result.dat <- dat - vec result.mat result.dat I hope it helps. Best, Dimitris Dimitris Rizopoulos Biostatistical Centre School of Public Health Catholic University of Leuven Address: Ka

Re: [R] elegant way to minus on each row of a matrix

2008-03-02 Thread Henrique Dallazuanna
Try this: sweep(mat, 1, vec) On 02/03/2008, Bo Zhou <[EMAIL PROTECTED]> wrote: > > How to do this in an elegant way formatrix/data frame/zoo? > > mat= > 1 2 3 > 4 5 6 > 7 8 9 > > vector= > 1 > 2 > 3 > > > result= > 0 1 2 > 2 3 4 > 4 5 6 > > ie > 1-1 2-1 3-1 > 4-2 5-2 6-2 > 7-3 8-3 9-3 > >

Re: [R] summarizing matrix data

2008-03-02 Thread Henrique Dallazuanna
One options is: data.1 <- matrix(seq(from=1,to=9,by=1),nrow=3,ncol=3) data.2 <- matrix(seq(from=11,to=19,by=1),nrow=3,ncol=3) data.3 <- matrix(seq(from=21,to=29,by=1),nrow=3,ncol=3) data.4 <- matrix(seq(from=31,to=39,by=1),nrow=3,ncol=3) data.5 <- matrix(seq(from=41,to=49,by=1),nrow=3,ncol=3) data

[R] Imputation Packages

2008-03-02 Thread ArunPrasad
Hi everyone, I am looking for a package in R which can help me in using the imputation technique to find the missing values for my regression analysis. Any help would be appreciated. Cheers Arun -- View this message in context: http://www.nabble.com/Imputation-Packages-tp15790749p

[R] regression output to latex

2008-03-02 Thread Luca Braglia
hello everybody i was seeking a converter beetween R regression output (eg with summary) and the conventional way to present regression output in paper: every model as a vertical vector with \beta, t beetween parenthesis below the first, and other statistics (R^2 etc) . I've seen hmisc and xtable

Re: [R] Variance Calculation in R

2008-03-02 Thread Keizer_71
sorry...in step 4-i need the R code to output in this format when i export to excel. ProbeID Variance 1 224588_at 21.58257457 thanks Keizer_71 wrote: > > Hello, > > Thanks everyone for helping me with the previous queries. > > step 1: Here is the orginal data: shor

[R] elegant way to minus on each row of a matrix

2008-03-02 Thread Bo Zhou
How to do this in an elegant way formatrix/data frame/zoo? mat= 1 2 3 4 5 6 7 8 9 vector= 1 2 3 result= 0 1 2 2 3 4 4 5 6 ie 1-1 2-1 3-1 4-2 5-2 6-2 7-3 8-3 9-3 Thanks in advance. _ 08 [[alternative HTML version

[R] Variance Calculation in R

2008-03-02 Thread Keizer_71
Hello, Thanks everyone for helping me with the previous queries. step 1: Here is the orginal data: short sample ProbeID Sample_1_D Sample_1_CSample_2_D Sample_2_C 1 224588_at 2.425509867 11.34031409 11.46868531 11.75741478 step 2: i calculate the v

[R] summarizing matrix data

2008-03-02 Thread Emilio Gagliardi
Hi everyone, I'm sure this is simple, but I can't seem to figure this out. Situation. 3 different groups of subjects each submit n X n matrices of scores. What I want to do is aggregate each group of scores into a summary n X n matrix. I need the result to be a matrix so that I can calculate a

Re: [R] COMPAR.GEE Output

2008-03-02 Thread Ben Bolker
Charles Willis gmail.com> writes: > > Hello, > > I am running the program COMPAR.GEE within the package APE. My dependent > variable is binomial, while my independent variable is a multi-state > categorical variable. The output reports an estimate for each state of the > independent variable ex

Re: [R] [OT] "normal" (as in "Guassian")

2008-03-02 Thread Ben Bolker
roger koenker ysidro.econ.uiuc.edu> writes: > > A nice survey of this territory is: > > http://books.google.com/books?id=TN3_d7ibo30C&pg=PA85&lpg=PA85&dq=stigler+normal+oxymoron&source=web&ots=OwGhmnDk3O&sig=J7ou_L8-_Mu4L14c3KJAhefrD4I&hl=en > > I particularly like the phrase: "[normal] is in

Re: [R] [OT] "normal" (as in "Guassian")

2008-03-02 Thread Gabor Csardi
I'm not a statistician, but do i remember well that among all distributions with a given mean and variance, the normal distribution has the highest entropy? This is good enough for me to call it "normal" Gabor On Sun, Mar 02, 2008 at 10:10:21AM -0600, roger koenker wrote: > A nice survey of

Re: [R] Newbie:Export Data into Excel from R

2008-03-02 Thread David Winsemius
Keizer_71 <[EMAIL PROTECTED]> wrote in news:[EMAIL PROTECTED]: > All i want is to export my list into c: drive and save it as csv > file and manually import into Excel. > > I have the read the article but i am having issues > http://pbil.univ-lyon1.fr/library/base/html/write.table.html > > >>

Re: [R] [OT] "normal" (as in "Guassian")

2008-03-02 Thread roger koenker
A nice survey of this territory is: http://books.google.com/books?id=TN3_d7ibo30C&pg=PA85&lpg=PA85&dq=stigler+normal+oxymoron&source=web&ots=OwGhmnDk3O&sig=J7ou_L8-_Mu4L14c3KJAhefrD4I&hl=en I particularly like the phrase: "[normal] is in this respect a rare one-word oxymoron." url:www.econ.

Re: [R] Newbie:Export Data into Excel from R

2008-03-02 Thread John Kane
R uses unix type slashes that is "/" rather than "\" even if you're using Windows You probably need either to write file ="c:/foo.csv" or escape the "\" thusly file = > "c:\\foo.csv" Either way should work. Personally find it easier to use '/". --- Keizer_71 <[EMAIL PROTECTED]> wrote: > > H

Re: [R] Newbie:Export Data into Excel from R

2008-03-02 Thread Charilaos Skiadas
You will likely need to escape that backslash, i.e. "c:\\foo.csv". Haris Skiadas Department of Mathematics and Computer Science Hanover College On Mar 2, 2008, at 10:12 AM, Keizer_71 wrote: > > Hi, > > All i want is to export my list into c: drive and save it as csv > file and > manually impor

Re: [R] Make plots with GNUplot. Have anyone tried that?

2008-03-02 Thread Louise Hoffman
> There is a very basic interface between R and gnuplot in the TeachingDemos > package. Look at the help for gp.plot. [snip] This looks mighty interesting =) Is it possible to plot "with lines" (gnuplot syntex) so all the data points are connected? Also is it possible to make it write the outpu

Re: [R] Idioms for a timeseries operation - moving window

2008-03-02 Thread Bo Zhou
Thanks. Problem solved. > rollapply(x,3,force, align="left") 1 1 2 3 2 2 3 4 3 3 4 5 4 4 5 6 5 5 6 7 6 6 7 8 7 7 8 9 8 8 9 10 > x<-1:10 > embed(x, 3)[, 3:1] [,1] [,2] [,3] [1,]123 [2,]234 [3,]345 [4,]456 [5,]567 [6,

[R] Newbie:Export Data into Excel from R

2008-03-02 Thread Keizer_71
Hi, All i want is to export my list into c: drive and save it as csv file and manually import into Excel. I have the read the article but i am having issues http://pbil.univ-lyon1.fr/library/base/html/write.table.html > excel<-write.table(probe_gene, file = "c:\foo.csv", sep = ",", col.names =

Re: [R] Idioms for a timeseries operation - moving window

2008-03-02 Thread Gabor Grothendieck
Try: rollapply(x, 3, force) or embed(x, 3)[, 3:1] On Sun, Mar 2, 2008 at 9:53 AM, Bo Zhou <[EMAIL PROTECTED]> wrote: > > I needed that matrix for following calculations. So I don't think roll* > would work for me. > > Re: embed > > > x <- 1:10 > > embed (x, 3) > [,1] [,2] [,3] > [1,]3

Re: [R] Idioms for a timeseries operation - moving window

2008-03-02 Thread Bo Zhou
I needed that matrix for following calculations. So I don't think roll* would work for me. Re: embed > x <- 1:10 > embed (x, 3) [,1] [,2] [,3] [1,]321 [2,]432 [3,]543 [4,]654 [5,]765 [6,]876 [7,]987 [8,]

Re: [R] Idioms for a timeseries operation - moving window

2008-03-02 Thread Gabor Grothendieck
Time to upgrade to the latest version of zoo. You can also read it online here: http://cran.r-project.org/web/packages/zoo On Sun, Mar 2, 2008 at 9:47 AM, Bo Zhou <[EMAIL PROTECTED]> wrote: > Ah! Gabor is the man. > > BTW I don't have zoo-faq. (I actually googled it out) > > This is what I have:

Re: [R] Idioms for a timeseries operation - moving window

2008-03-02 Thread Bo Zhou
Ah! Gabor is the man. BTW I don't have zoo-faq. (I actually googled it out) This is what I have: Vignettes in package ‘zoo’: zoo-quickrefzoo Quick Reference (source, pdf) zoo zoo: An S3 Class and Methods for Indexed Totally Ord

Re: [R] Idioms for a timeseries operation - moving window

2008-03-02 Thread Gabor Grothendieck
See ?embed and from zoo see: ?rollapply ?rollmean ?rollmax There is a function coded in C in the caTools package for speed. On Sun, Mar 2, 2008 at 9:24 AM, Bo Zhou <[EMAIL PROTECTED]> wrote: > > Hi Guys, > > Need your wisdom on this. > > Say I have a time series (in zoo format) like this > > >

  1   2   >