[R] difference between lrm's Model L.R. and anova's Chi-Square

I am running lrm() with a single factor. I then run anova() on the fitted model to obtain a p-value associated with having that factor in the model. I am noticing that the Model L.R. in the lrm results is almost the same as the Chi-Square in the anova results, but not quite; the latter value is

Re: [R] difference between lrm's Model L.R. and anova's Chi-Square

Quoting Frank E Harrell Jr [EMAIL PROTECTED]: anova (anova.Design) computes Wald statistics. When the log-likelihood is very quadratic, these statistics will be very close to log-likelihood ratio chi-square statistics. In general LR chi-square tests are better; we use Wald tests for speed.

Re: [R] Poisson regression in R

glmstat wrote: I have these questions: (1) Use Poisson regression to estimate the main effects of car, age, and dist (each treated as categorical and modelled using indicator variables) and interaction terms. (2) It was determined by one study that all the interactions were unimportant and

Re: [R] Could not install aroma.affymetrix

Hi My, On Sat, Mar 1, 2008 at 8:04 PM, My Coyne [EMAIL PROTECTED] wrote: I don't know if this is the correct forum to ask the following question; however, when I search the aroma.affymetrix discussion group, it suggested that I should posted the question to r-help. Here it goes. I think

Re: [R] Newbie: Incorrect number of dimensions

Here is the structure str(all_differ) num [1:6280, 1:8] 2.22e-16 2.22e-16 2.22e-16 2.22e-16 2.22e-16 ... - attr(*, dimnames)=List of 2 ..$ : NULL ..$ : chr [1:8] rawp Bonferroni Holm Hochberg ... I tried, but it still giving me the same error message Incorrect Dimension

Re: [R] Newbie: Incorrect number of dimensions

You have a 'matrix' not a 'list', so I think you want something like this: all_differ[all_differ[, 'BY'] = 0.01, 'BY'], this will return the values where column is 'BY' in all rows which is less or equal to 0.01. On 02/03/2008, Keizer_71 [EMAIL PROTECTED] wrote: Here is the structure

Re: [R] jpeg and margin text

Thomas Schwander wrote: Hi guys, I use R 2.6.2 and Windows XP. I’ve got the following question: I wrote a lot of text into the margin of a plot. When I use the “jpeg”-function, only a little part of the margin text is displayed in the final jpeg? Any ideas to change it?

Re: [R] an efficient pairwise matrix cell's comparison function

Does this do what you want? A - matrix(sample(0:2, 25, TRUE), ncol=5) B - matrix(1:25, ncol=5) C - ifelse(A == 0, 0, B) A [,1] [,2] [,3] [,4] [,5] [1,]11121 [2,]10110 [3,]00102 [4,]01200 [5,]121

Re: [R] an efficient pairwise matrix cell's comparison function

C - B C[A==0] - 0 would be somewhat more efficient. On Sun, 2 Mar 2008, jim holtman wrote: Does this do what you want? A - matrix(sample(0:2, 25, TRUE), ncol=5) B - matrix(1:25, ncol=5) C - ifelse(A == 0, 0, B) A [,1] [,2] [,3] [,4] [,5] [1,]11121 [2,]10

[R] [OT] normal (as in Guassian)

Hi Folks, Apologies to anyone who'd prefer not to see this query on this list; but I'm asking because it is probably the forum where I'm most likely to get a good answer! I'm interested in the provenance of the name normal distribution (for what I'd really prefer to call the Gaussian

Re: [R] Newbie: Incorrect number of dimensions

This is exactly what i am looking for. many thanks Chris On 3/2/08, Henrique Dallazuanna [EMAIL PROTECTED] wrote: You have a 'matrix' not a 'list', so I think you want something like this: all_differ[all_differ[, 'BY'] = 0.01, 'BY'], this will return the values where column is 'BY' in all

[R] question on lag.zoo

Hi Guys, I'm using zoo package now. I found lag is not doing what I assumed. x - zoo(11:21) z - zoo(1:10, yearqtr(seq(1959.25, 1961.5, by = 0.25)), frequency = 4) x 1 2 3 4 5 6 7 8 9 10 11 11 12 13 14 15 16 17 18 19 20 21 lag(x) 1 2 3 4 5 6 7 8 9 10 12 13 14 15 16 17

Re: [R] difference between lrm's Model L.R. and anova's Chi-Square

[EMAIL PROTECTED] wrote: Quoting Frank E Harrell Jr [EMAIL PROTECTED]: anova (anova.Design) computes Wald statistics. When the log-likelihood is very quadratic, these statistics will be very close to log-likelihood ratio chi-square statistics. In general LR chi-square tests are better; we

Re: [R] question on lag.zoo

On Sun, Mar 2, 2008 at 9:28 AM, Gabor Grothendieck [EMAIL PROTECTED] wrote: Its a FAQ. x is of class zoo but z is of class zooreg so one is using That should be lag.zooreg zoo.lag and the is using zooreg.lag. See question #6 in the zoo FAQ: vignette(zoo-faq) and also ?lag.zoo On

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

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

[R] Idioms for a timeseries operation - moving window

Hi Guys, Need your wisdom on this. Say I have a time series (in zoo format) like this x - zoo(11:21) x 1 2 3 4 5 6 7 8 9 10 11 11 12 13 14 15 16 17 18 19 20 21 I want to do a moving window sampling of it. The result can either be a matrix or a dataframe like this

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

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:

[R] Newbie:Export Data into Excel from R

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 = NA)

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

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

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

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

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

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] Make plots with GNUplot. Have anyone tried that?

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 output

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

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: Hi, All i

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

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

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

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] [OT] normal (as in Guassian)

A nice survey of this territory is: http://books.google.com/books?id=TN3_d7ibo30Cpg=PA85lpg=PA85dq=stigler+normal+oxymoronsource=webots=OwGhmnDk3Osig=J7ou_L8-_Mu4L14c3KJAhefrD4Ihl=en I particularly like the phrase: [normal] is in this respect a rare one-word oxymoron. url:

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

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] COMPAR.GEE Output

Charles Willis willis.charlie at 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

[R] summarizing matrix data

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

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

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

Re: [R] Variance Calculation in R

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

[R] regression output to latex

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

[R] Variance Calculation in R

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

[R] Imputation Packages

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:

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

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,]05

Re: [R] Variance Calculation in R

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

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] elegant way to minus on each row of a matrix

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 Thanks in advance.

Re: [R] summarizing matrix data

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

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

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

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:

Re: [R] Variance Calculation in R

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

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

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

Re: [R] summarizing matrix data

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

Re: [R] Variance Calculation in R

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

Re: [R] Imputation Packages

ArunPrasad agurubar at 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

[R] new to latex to pdf

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

[R] About R-Project

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] [OT] normal (as in Guassian)

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

[R] Need help to locate my mistake

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] summarizing matrix data

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)), dim

Re: [R] Variance Calculation in R

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

[R] discrete variable

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

[R] emacs and R

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] emacs and R

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

Re: [R] discrete variable

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, 4, 3, 15,

Re: [R] summarizing matrix data

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

Re: [R] discrete variable

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 many

Re: [R] discrete variable

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 vector

Re: [R] Need help to locate my mistake

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(fuelData)[1]

Re: [R] Need help to locate my mistake

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)[1] xOnes-

[R] listing components of an object

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,

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

Johannes Hüsing johannes at 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

Re: [R] [PS] discrete variable

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

Re: [R] emacs and R

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 R. Can

Re: [R] listing components of an object

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), ## using

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

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 the

Re: [R] emacs and R

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 how

[R] Recommended Packages

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

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

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 query

Re: [R] discrete variable

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

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

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,

[R] Runtime error

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

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

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

Re: [R] Recommended Packages

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

Re: [R] Need help to locate my mistake

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

Re: [R] About R-Project

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.html and

[R] insert vector image into a plot

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

Re: [R] emacs and R

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

Re: [R] Need help to locate my mistake

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

Re: [R] insert vector image into a plot

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

[R] R data Export to Excel

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

[R] glm: offset

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)

Re: [R] R data Export to Excel

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] unable to start device PNG and unable to open connection to X11 display

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

Re: [R] glm: offset

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 (1) a

Re: [R] R data Export to Excel

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 #output variance

Re: [R] R data Export to Excel

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

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] plotting to stdout (while reading from stdin?)

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 works

[R] R function to convert a number to text

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] insert vector image into a plot

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

[R] Constrained regression

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

[R] handling big data set in R

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

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,

Re: [R] Problem with PCA

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

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

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(png::,

Re: [R] R data Export to Excel

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

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

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

[R] Formating a zoo dataset .

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

  1   2   >