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

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 spe

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

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 thin

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

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

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

To all, I am undergoing an analysis involving big matrices of about 3x200 which I have to handle in a more efficient way. So I would like some advice to build such efficient function to deliver the following result: - starting with 2 matrices of the same dimension (eg. A and B)

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

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

[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

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 bette

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.

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

[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 my.super

Re: [R] question on lag.zoo

Its a FAQ. x is of class "zoo" but z is of class "zooreg" so one is using zoo.lag and the is using zooreg.lag. See question #6 in the zoo FAQ: vignette("zoo-faq") and also ?lag.zoo On Sun, Mar 2, 2008 at 9:12 AM, Bo Zhou <[EMAIL PROTECTED]> wrote: > > Hi Guys, > > I'm using zoo package now

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 =

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 Ord

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

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

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 outpu

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

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 impor

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

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

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

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

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

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

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

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

[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: http://www.nabble.com/Imputation-Packages-tp15790749p

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

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

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

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

Re: [R] summarizing matrix data

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

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

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 PROTECTE

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 PROT

Re: [R] Imputation Packages

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

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

[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 rotation-sym

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

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

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

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

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,

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:

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 man

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

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

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

[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, Species

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

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

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

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

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 th

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

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

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 q

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

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

[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 error*

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

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

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 opening

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

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 ridicu

Re: [R] Need help to locate my mistake

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

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

[R] reducing RODBC odbcQuery memory use?

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] Need help to locate my mistake

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

Re: [R] Recommended Packages

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] 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.h

Re: [R] emacs and R

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

[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 https://stat.ethz.ch/mailman/listinfo/r-

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 PR

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-

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 ca

[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], Variance=apply(data[,y],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 PN

Re: [R] glm: offset

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.

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