Re: [R] Converting the graphics window to a data matrix

See library(grid) ?grid.cap and this page for background: http://developer.r-project.org/Raster/raster-RFC.html On Thu, May 19, 2011 at 2:19 PM, Beutel, Terry S < terry.beu...@deedi.qld.gov.au> wrote: > I'm wondering if anyone is aware of a way to take what is visible in the > R graphics win

Re: [R] Add a vector to the values in a specific dimension of an array

It depends on how big A is and how much memory you have. Here is one lazy way. A <- aperm(aperm(A, c(2,1,3,4)) + x, c(2,1,3,4)) Bill Venables -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of huron Sent: Thursday, 19 May 2011 1:47

Re: [R] Converting the graphics window to a data matrix

On Thu, May 19, 2011 at 5:19 AM, Beutel, Terry S wrote: > I'm wondering if anyone is aware of a way to take what is visible in the > R graphics window, pixelate it, and express that pixellation in a > numeric matrix. For example, I am using the following command to create > a  black (present) and

Re: [R] Help, please

On Wed, 18 May 2011, Julio César Flores Castro wrote: Hi, I am using R 2.10.1 and I have a doubt. Do you know how many cases can R handle? At least millions. I want to use the library npmc but if I have more than 4,500 cases I get an error message. If I use less than 4500 cases I don?t ha

Re: [R] Add a vector to the values in a specific dimension of an array

Hi: Reordering the dimensions, then doing a vectorized addition, then reordering (back) again is faster, it seems. > m <- 20; n <- 30; p <- 40; q <- 30 > a <- NA > length(a) <- m * n * p * q > dim(a) <- c(m, n, p, q) > x <- 1:n > a[1:m,,1:p,1:q] <- 0 > b <- a > > # Approach 1 > system.time({ +

[R] Converting the graphics window to a data matrix

I'm wondering if anyone is aware of a way to take what is visible in the R graphics window, pixelate it, and express that pixellation in a numeric matrix. For example, I am using the following command to create a black (present) and white (absent) spatial simulation of grass tussock distribution.

[R] Add a vector to the values in a specific dimension of an array

Hello, A simple question, although I can't find an answer via my google/forum search: I have a 4-dimensional array; call it A[1:M,1:N,1:P,1:Q]. I have a vector x that is N by 1. I would like to "quickly" add x to the 2nd dimension of A; in other words, I want a quicker way of doing the followi

Re: [R] Help, please

Hi, On Wed, May 18, 2011 at 6:29 PM, Julio César Flores Castro wrote: > Hi, > > I am using R 2.10.1 and I have a doubt. As a general rule of thumb, it's usually best to be using the latest version of R (which is no R 2.13.0) -- 2.10.1 came out in December, 2009, so ... try to upgrade if you can.

Re: [R] Covariable Logistic Regression In R

In my humble opinion, the statistical methods should be mastered before the software. Frank Anamaria Crisan wrote: > > I specifically wanted to do an adjusted covariable analysis and I want > incorporate genetic information in that. The language tossed around is > just > to do an 'adjusted analys

Re: [R] Finding row in matrix with minimum column

Hi Worik, See ?which.min x <- matrix(c(7, 7, 9, 7, 7, 9, 2, 9), ncol = 2, byrow = FALSE) which.min(x[,2]) x[which.min(x[,2]), ] HTH, Jorge On Wed, May 18, 2011 at 10:38 PM, Worik R <> wrote: > Friends > > If I have a matrix such as... > > [,1] [,2] > [1,]77 > [2,]79 > [3,]

[R] Finding row in matrix with minimum column

Friends If I have a matrix such as... [,1] [,2] [1,]77 [2,]79 [3,]9 2 [4,]79 And I want to find the row number that has the minimum value of column 2 (row 3 in this case) how can I do it? Is there a simple way? cheers Worik [[alternative HTML versio

Re: [R] identical function names from 2 packages

Thanks!! Nick On 5/18/11 7:06 PM, Ista Zahn wrote: Hi Nick, See ?"::" Best, Ista On Wed, May 18, 2011 at 10:02 PM, Nick Matzke wrote: Hi, If I load 2 packages that have a function with the same name, how do I tell R to run one or the other? (Instead of having R automatically use the first-

Re: [R] R Style Guide -- Was Post-hoc tests in MASS using glm.nb

On Wed, May 18, 2011 at 6:28 PM, David Scott wrote: >  Another style guide is at: > http://www1.maths.lth.se/help/R/RCC/ > > Listed as a first draft and dated 2005, but still worth a read. Has some > references also. That URL obsolete (I need to have it removed) - a more recent/stable URL is [5]

Re: [R] identical function names from 2 packages

Hi Nick, See ?"::" Best, Ista On Wed, May 18, 2011 at 10:02 PM, Nick Matzke wrote: > Hi, > > If I load 2 packages that have a function with the same name, how do I tell > R to run one or the other? > > (Instead of having R automatically use the first- or last-loaded one, > whichever it is.  (Whi

Re: [R] identical function names from 2 packages

Is it not something like? package1::myfunction package2::myfunction On Wednesday, May 18, 2011 at 9:02 PM, Nick Matzke wrote: > Hi, > > If I load 2 packages that have a function with the same > name, how do I tell R to run one or the other? > > (Instead of having R automatically use the first

Re: [R] Parameter Estimates needed from lmer output

David Cross tcu.edu> writes: > > If I remember correctly, coef(m1) would do it ... > but it has been a while since I last used lmer, and I am > working only from memory. > fixef() for fixed effects only ranef() for random effects only coef() gives a matrix of the predicted coefficients in

[R] identical function names from 2 packages

Hi, If I load 2 packages that have a function with the same name, how do I tell R to run one or the other? (Instead of having R automatically use the first- or last-loaded one, whichever it is. (Which is it, by the way.)) Cheers! Nick -- ==

Re: [R] matrix help (first occurrence of variable in column)

Is this what you were after: > mdat <- matrix(c(1,0,1,1,1,0), nrow = 2, ncol=3, byrow=TRUE, + dimnames = list(c("T1", "T2"), + c("sp.1", "sp.2", "sp.3"))) > > mdat sp.1 sp.2 sp.3 T1101 T2110 > # do 'rle' on each column and see

Re: [R] R Style Guide -- Was Post-hoc tests in MASS using glm.nb

Another style guide is at: http://www1.maths.lth.se/help/R/RCC/ Listed as a first draft and dated 2005, but still worth a read. Has some references also. I think I recall Hadley having a style guide which he requested his students followed, but I didn't like it too much (sorry Hadley) . I

Re: [R] issue with odfWeave running on Windows XP; question about installing packages under Linux

Sorry for the delayed response. An upgrade of the XML package has broken odfWeave; see this thread: https://stat.ethz.ch/pipermail/r-help/2011-May/278063.html That may be your issue. We're working on the problem now. I'll post to R-Packages when we have a working update. If you like, I can se

Re: [R] scales argument in bwplot (lattice)

Hi Peter A little late but catching up see ? combineLimits from the latticeExtra package a very welcome addition in particular when combined with useOuterStrips with multiple conditioning Regards Duncan Duncan Mackay Department of Agronomy and Soil Science University of New England ARMIDAL

Re: [R] R Style Guide -- Was Post-hoc tests in MASS using glm.nb

I used to think like that. However I have recently re-read John Chambers' "Software for Data Analysis" and now I'm starting to see the point. S4 classes and methods do require you to plan your classes and methods well and the do impose a discipline that can seem rigid and unnecessary. But I ha

Re: [R] Parameter Estimates needed from lmer output

If I remember correctly, coef(m1) would do it ... but it has been a while since I last used lmer, and I am working only from memory. Cheers David Cross d.cr...@tcu.edu www.davidcross.us On May 18, 2011, at 6:29 PM, Stephen Peterson wrote: > Hello, > I am looking for some help on how I may b

[R] Ordinal Package

Hi, Does anyone know if Ordinal package can be used for continuous covariate. In CLMM function, can independent variables be continuous variables? I used all dummy variable as the independent variables, the CLMM function works, but when I added continuous variables in the model, it gave an error.

[R] Parameter Estimates needed from lmer output

Hello, I am looking for some help on how I may be able to view estimated values for 3 response variables with 1 fixed and 1 random effect using lmer. My data is proportional cover of three habitat variables (bare ground, grass cover, shrub cover) that was collected during 3 years (1976, 2000, 2010)

[R] Help, please

Hi, I am using R 2.10.1 and I have a doubt. Do you know how many cases can R handle? I want to use the library npmc but if I have more than 4,500 cases I get an error message. If I use less than 4500 cases I don´t have problems with this library. Is there any way to increase the number of cas

Re: [R] Covariable Logistic Regression In R

I specifically wanted to do an adjusted covariable analysis and I want incorporate genetic information in that. The language tossed around is just to do an 'adjusted analysis', but I am not certain if this is common terminology as I am relatively new to biostatistics. Perhaps I did not express myse

Re: [R] R Style Guide -- Was Post-hoc tests in MASS using glm.nb

On Thu, 19 May 2011, Timothy Bates wrote: lol:-) is there a way to get a fortune on start up in R? One possible way is to put something like this into your .Rprofile if(interactive()) { library("fortunes") fortune() } The interactive() condition avoids having it in the output of certain

Re: [R] R Style Guide -- Was Post-hoc tests in MASS using glm.nb

On Wed, May 18, 2011 at 4:44 PM, Timothy Bates wrote: > lol:-) >  is there a way to get a fortune on start up in R? You could require() and fortune() one in your .Rprofile > On 19 May 2011, at 00:33, Rolf Turner wrote: > >> >> On 19/05/11 10:26, bill.venab...@csiro.au wrote: >> >> >>> Most of [

Re: [R] R Style Guide -- Was Post-hoc tests in MASS using glm.nb

On Wed, May 18, 2011 at 4:33 PM, Rolf Turner wrote: > On 19/05/11 10:26, bill.venab...@csiro.au wrote: > >> >> Most of [the Google style guide's] advice is very good (meaning I agree >> with it!) but some is a bit too much (for example, the blanket advice never >> to use S4 classes and methods -

Re: [R] R Style Guide -- Was Post-hoc tests in MASS using glm.nb

lol:-) is there a way to get a fortune on start up in R? On 19 May 2011, at 00:33, Rolf Turner wrote: > > On 19/05/11 10:26, bill.venab...@csiro.au wrote: > > >> Most of [the Google style guide's] advice is very good (meaning I agree with >> it!) but some is a bit too much (for example, the b

Re: [R] R Style Guide -- Was Post-hoc tests in MASS using glm.nb

On 19/05/11 10:26, bill.venab...@csiro.au wrote: Most of [the Google style guide's] advice is very good (meaning I agree with it!) but some is a bit too much (for example, the blanket advice never to use S4 classes and methods - that's just resisting progress, in my view). I must respectf

[R] using hglm to fit a gamma GLMM with nested random effects?

Apologies for continuing to ask about this but . . in my quest to fit a gamma GLMM model to my data (see partial copy of thread below), I'm exploring using hglm today. The question of the day has to do with the errors I'm currently getting from the hglm package. Can hglm handle a model with nested

Re: [R] R Style Guide -- Was Post-hoc tests in MASS using glm.nb

Hi Bert, I think people should know about the Google Sytle Guide for R because, as I said, it represents a thoughtful contribution to the debate. Most of its advice is very good (meaning I agree with it!) but some is a bit too much (for example, the blanket advice never to use S4 classes and m

Re: [R] Plots: I've deleted axes, now to delete space

I've never used an R command without reading the documentation first. I think that would be impractical. I am not an expert at deciphering the documentation though and I post here because I did in fact read the documentation, do extensive google searches, ask friends/collegues and still find no

Re: [R] Covariable Logistic Regression In R

First of all, use the correct terminology for the statistical method. It is not productive to identify 'significant' genes unless you use quite complex methods. This kind of endeavor would take at least 6 statistics courses in order to be successful. Frank Anamaria Crisan wrote: > > Hello, > >

Re: [R] Plots: I've deleted axes, now to delete space

On 05/16/2011 06:18 PM, adele_thomp...@cargill.com wrote: Re-sizing within the dev command works well. I'm not sure why I would need the dev.off(). I have the plot commands run. Then I have the dev.copy2pdf command. Thanks again for your help. Well, you really should get into the habit of rea

Re: [R] Help with Memory Problems (cannot allocate vector of size)

How about reading the posting guide and afterwards searching the list archive: http://r.789695.n4.nabble.com/R-help-f789696.html Searching for "Error: cannot allocate vector of size" will give you hundreds of results as this question is asked VERY frequently! Jannis On 05/18/2011 07:55 P

[R] June ** R ** /S+ Courses: Nationwide Back2back (1) R/S+ Fundamentals and (2) R/S-Plus Advanced Programming. in San Francisco, New York City, Washington DC

XLSolutions has scheduled the first 2011 back2back courses in New York City, San Francisco, Washington DC. Taught by top R/S+ gurus! West Coast <---back2back---> East Coast (1) R/S-PLUS Fundamentals and Programming Techniques http://www.xlsolutions-corp.com/coursedetail.asp?id=30 * San Franc

Re: [R] ggplot geom_boxplot vertical margins

Exactly! Thanks, I couldn't find that anywhere! On Wed, May 18, 2011 at 1:59 PM, Felipe Carrillo wrote: > Is this what you want? You can control how much space you > want to see on the sides of the plot: > > df<-data.frame(x=factor(1:100),y=rnorm(1000)) > ggplot(df,aes(x=x,y=y))+geom_boxplot() +

Re: [R] ggplot geom_boxplot vertical margins

Is this what you want? You can control how much space you want to see on the sides of the plot: df<-data.frame(x=factor(1:100),y=rnorm(1000)) ggplot(df,aes(x=x,y=y))+geom_boxplot() + scale_x_discrete(expand=c(0,0))   Felipe D. Carrillo Supervisory Fishery Biologist Department of the Interior US

[R] ggplot geom_boxplot vertical margins

If you plot: df<-data.frame(x=factor(1:100),y=rnorm(1000)) ggplot(df,aes(x=x,y=y))+geom_boxplot() How do I remove those pesky margins on the sides of the plot area? Or maybe just reduce their size to something more like the spacing of the boxes? Thanks, Justin ___

Re: [R] Convolution confusion:

On May 18, 2011, at 15:47 , Alex Hofmann wrote: > Hi, > > I'm new to R, and I'm a bit confused with the "convolve()" function. > If I do: > x<-c(1, 2, 3) > convolve(x, rev(x), TRUE, "open") > = 9 12 10 4 1 > > But I expected: 3 8 14 8 3 (like in Octave/MATLAB - conv(x, reverse(x)) ) > > 3 2 1

[R] Covariable Logistic Regression In R

Hello, I would like some help figuring out how to run Covariable Logistic Regression in R. I've been searching for a while on how to get this done in R (I have had the luck previously of using a software package that just does it) and I am coming up empty handed. Any experience or insights would b

Re: [R] data network format and grouping analysis

The following works to get an igraph object from a matrix edgelist: dat2 <- matrix(rep(seq(1,5,1), 4), nrow=10, ncol=2) graph.edgelist( dat2 ) I tried with NA's but graph.edgelist did not allow NA's. Wouldn't you just leave those rows out with NA's in them? An NA means there is no edge, right?

Re: [R] Overlaying maps

Using my mind-reading skills, plot.Map(), while deprecated, does provide an add= option. If this doesn't help, you'll need to read the posting guide and provide (a lot) more information. Ray Brownrigg On Thu, 19 May 2011, michael.laviole...@dhhs.state.nh.us wrote: > I'm having difficulty over

Re: [R] Grouped bar plot

On Wed, May 18, 2011 at 2:38 PM, Peter Ehlers wrote: > On 2011-05-18 11:13, jctoll wrote: >> >> Hi, >> >> I am trying to produce a grouped bar plot from a data.frame and I'm >> having difficulties figuring out how to do so.  My data is 500 rows by >> 4 columns and basically looks like so: >> >>> h

Re: [R] Grouped bar plot

On 2011-05-18 11:13, jctoll wrote: Hi, I am trying to produce a grouped bar plot from a data.frame and I'm having difficulties figuring out how to do so. My data is 500 rows by 4 columns and basically looks like so: head(x) V1V2V3V4 1 XOM 0.2317915 0.1610068 1.6

Re: [R] Use paste function to select column of data

On 2011-05-18 12:12, John Poulsen wrote: Hello, I want to build a function to call up a column of a data.frame by the names of the columns. I have column names that are sequentially named (col1, col2, etc.). How do I change a character expression into something that will be understood as a d

Re: [R] Use paste function to select column of data

On May 18, 2011, at 3:12 PM, John Poulsen wrote: Hello, I want to build a function to call up a column of a data.frame by the names of the columns. I have column names that are sequentially named (col1, col2, etc.). How do I change a character expression into something that will be unde

Re: [R] assign $y of predict() function output to variable

On May 18, 2011, at 1:35 PM, Asan Ramzan wrote: Hello R-help Below is the output from the predict() function. How can I assign $y to a variable. Newvar <- predict(function,df2)$y -- David. $x V1 1 36.28 2 34.73 3 33.74 4 69.87 5 58.88 6 89.44 7 43.97 8 41.94 9

[R] Use paste function to select column of data

Hello, I want to build a function to call up a column of a data.frame by the names of the columns. I have column names that are sequentially named (col1, col2, etc.). How do I change a character expression into something that will be understood as a data.frame column. For example: example<

Re: [R] Simple ordering or sorting question

That did it. Thanks!! David === David Kaplan, Ph.D. Professor Department of Educational Psychology University of Wisconsin - Madison Educational Sciences, Room, 1082B 1025 W. Johnson Street Madison, WI 53706 email: dkap...@edu

Re: [R] Simple ordering or sorting question

It looks like your row numbers are "characters" because that is the sort sequence you are getting. Try df <- df[order(as.numeric(rownames(df))), ] On Wed, May 18, 2011 at 2:42 PM, David Kaplan wrote: > Greetings, > > I'm trying to simply reorder a data frame on the row numbers.  So, for > exampl

[R] Simple ordering or sorting question

Greetings, I'm trying to simply reorder a data frame on the row numbers. So, for example, instead of getting 1,2,3,4,5,6,7,8,9,10,11, ... 100 ..., I get instead 1, 10, 100, 101, 102, 103, 104, 105, 106, 107, 108, 109, 11, ... I've tried commands such as df <- df[order(rownames(df)),] and

Re: [R] Dataset Quasi Poisson

1) This mailing list is not for homework. 2) I would recommend reading the introduction to R that comes with every installation of R, since your answer is in there. Alternatively, you could google R and quasi poisson. On Wed, May 18, 2011 at 11:42 AM, ilpoeta84 wrote: > Hello, I'm looking for a

[R] data network format and grouping analysis

Hi everyone, I have a dataset of friendship with this format: ego alter 47461 2 97421 3 14738 1NA 47472NA 974323 14739 21 4748313 97443 5 14740 314 47494NA 97454NA 14741 4NA 47505NA 9746513

[R] Grouped bar plot

Hi, I am trying to produce a grouped bar plot from a data.frame and I'm having difficulties figuring out how to do so. My data is 500 rows by 4 columns and basically looks like so: > head(x) V1V2V3V4 1 XOM 0.2317915 0.1610068 1.6941637 2 AAPL 0.6735488 0.7433611 0.15

Re: [R] hierarchical clustering within a size limit

Hi Peter, Thanks for your help. A second simple question that I cannot solve is the following. labels = cutree(hc, h=500) # members of cluster 1: x[labels==1] # members of cluster 2: x[labels==2] When x is >= 8 the index numbers appear in the output: [['[1]', '180066408', '180066464', '18006646

Re: [R] R-code in R-file documentation

I guess what you want is cat(readLines("file.r"), sep = "\n") Regards, Yihui -- Yihui Xie Phone: 515-294-2465 Web: http://yihui.name Department of Statistics, Iowa State University 2215 Snedecor Hall, Ames, IA On Wed, May 18, 2011 at 4:03 AM, Brian Oney wrote: > Hello List, > I would like to

[R] Help with Memory Problems (cannot allocate vector of size)

While doing pls I found the following problem > BHPLS1 <- plsr(GroupingList ~ PCIList, ncomp = 10, data = PLSdata, jackknife > = >FALSE, validation = "LOO") when not enabling jackknife the command works fine, but when trying to enable jackknife i get the following error. >BHPLS1 <- plsr(Gr

Re: [R] strucchange package Linux help

On Wed, 18 May 2011, Hasan Diwan wrote: When I run the code below on Macintosh and Windows, the plot comes out fine. However, on Linux, the png generated is invalid from R console, and loading strucchange crashes rkward. I can replicate nothing of this. I ran the script both in a plain R 2.13.

[R] text mining problem using TM package

Hi, I’m using R (TM package) for text mining and I’m having problems filtering articles out of my data set by local meta data. Here is the code: *data <- ("C:/… /19970331")* * * * * *rs <- ReutersSource(data , encoding = "UTF-8")* *RC <- VCorpus(DirSource(data), readerControl = list(reade

[R] Email out of R (code)

In case you're using Unix/Linux, have a look at www.r-project.org/doc/Rnews/Rnews_2007-1.pdf (page 30 - 32) Wolfgang . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . Wolfgang Raffelsberger, PhD IGBMC, 1 rue Laurent Fries, 67404 Illkirch Strasbourg, France Tel (+33) 388 6

[R] Overlaying maps

I'm having difficulty overlaying maps when writing to a file graphics device. My command sequence has the structure plot(map1) par(new = T) plot(map2) On the screen device, it works fine. When I attempt something like png(file = "map.png") plot(map1) par(new = T) plot(map2) dev.off() only the

[R] assign $y of predict() function output to variable

Hello R-help Below is the output from the predict() function. How can I assign $y to a variable. >predict(function,df2) $x    V1 1   36.28 2   34.73 3   33.74 4   69.87 5   58.88 6   89.44 7   43.97 8   41.94 9   33.34 10  38.47 11  35.16 12  42.94 13  46.76 14  53.24 15  52.43 16  50.40 17

Re: [R] Smooth contour of a map

It's perfect, thank you! I would like post the final code if someone need help in this subject, but I try to correct a last problem, how can I constrain the contourLines() function to take the corner points of the map in his result ... it does not consider this point like a contour point. Le 1

Re: [R] logistic regression lrm() output

Why is a one unit change in x an interesting range for the purpose of estimating an odds ratio? The default in summary() is the inter-quartile-range odds ratio as clearly stated in the rms documentation. Frank array chip wrote: > > Hi, I am trying to run a simple logistic regression using lrm()

Re: [R] Help with 2-D plot of k-mean clustering analysis

One idea: Pick the three largest clusters, their centers determine a plane. project your data into that plane. albyn On Wed, May 18, 2011 at 06:55:39PM +0200, Claudia Beleites wrote: > Hi Meng, > > > I would like to use R to perform k-means clustering on my data which > >included 33 samples me

Re: [R] Help with 2-D plot of k-mean clustering analysis

Hi Meng, I would like to use R to perform k-means clustering on my data which included 33 samples measured with ~1000 variables. I have already used kmeans package for this analysis, and showed that there are 4 clusters in my data. However, it's really difficult to plot this cluster in 2-D for

Re: [R] Email out of R (code)

I do not know. I was not aware and could hardly find any information on create.post(). From what I have seen at first glance, it seems that create.post() either opens your standard email program or web browser, which the python code does not. Instead it needs the R-library interfacing Python. I als

[R] strucchange package Linux help

When I run the code below on Macintosh and Windows, the plot comes out fine. However, on Linux, the png generated is invalid from R console, and loading strucchange crashes rkward. Is this a known issue on Linux and, if so, is there a workaround? Many thanks! require(strucchange) data("RealInt") bp

[R] logistic regression lrm() output

Hi, I am trying to run a simple logistic regression using lrm() to calculate a odds ratio. I found a confusing output when I use summary() on the fit object which gave some OR that is totally different from simply taking exp(coefficient), see below: > dat<-read.table("dat.txt",sep='\t',header=T

Re: [R] Loop stopping after 1 iteration

Didn't mean to snub you guys, Hugo and David. I didn't see your posts before. Thanks for the advice. -- View this message in context: http://r.789695.n4.nabble.com/Loop-stopping-after-1-iteration-tp3532988p3533217.html Sent from the R help mailing list archive at Nabble.com. _

Re: [R] Email out of R (code)

How does this compare to create.post() ? Kevin On Tue, May 17, 2011 at 3:44 PM, Daniel Malter wrote: > Hi all, > > I thought I would post code to send an email out of R. The code uses > Grothendieck and Bellosta's interface package rJython for executing Python > from R. The code itself provide

[R] Dataset Quasi Poisson

Hello, I'm looking for a dataset for Quasipoisson regression. The result must be significantly different from the classic poisson regression. You can help me? Please It is for my last university exam Thanks a lot -- View this message in context: http://r.789695.n4.nabble.com/Dataset-Quasi-Poiss

Re: [R] Loop stopping after 1 iteration

William, num_obs obviously isn't a vector, therefore length(num_obs) will evaluate to one. Hence your for loop control part will expand to for (i in 1:1) while it should probably read: for (i in 1:num_obs) Best Hugo On Wednesday 18 May 2011 17:18:15 armstrwa wrote: > Hi all

Re: [R] Loop stopping after 1 iteration

On May 18, 2011, at 11:18 AM, armstrwa wrote: Hi all, This is a very basic question, but I just can't figure out why R is handling a loop I'm writing the way it is. Here is the script I have written: grid_2_series<-function(gage_handle,data_type,filename) series_name<-paste(gage_handle,d

Re: [R] Loop stopping after 1 iteration

I knew it would be something simple. Thanks for catching that, Martyn. Billy -- View this message in context: http://r.789695.n4.nabble.com/Loop-stopping-after-1-iteration-tp3532988p3533041.html Sent from the R help mailing list archive at Nabble.com. __

Re: [R] Loop stopping after 1 iteration

Hi, The answer to (2) is that num_obs is a scalar, so length(num_obs) is 1. You probably wanted to do for (i in 1:num_obs) instead. Best wishes Martyn -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of armstrwa Sent: 18 May 2011 1

Re: [R] Help with 2-D plot of k-mean clustering analysis

On Wed, May 18, 2011 at 7:41 AM, Meng Wu wrote: > Hi, all > >  I would like to use R to perform k-means clustering on my data which > included 33 samples measured with ~1000 variables. I have already used > kmeans package for this analysis, and showed that there are 4 clusters in my > data. Howeve

Re: [R] Integral Symbol

Can’t you just embed it in the html as a symbol? ∫ or ∫ I’d have thought you also just put it into straight into the document as a character – ∫– , as long as the html is stored as unicode U+222B http://en.wikipedia.org/wiki/Integral_symbol On 18 May 2011, at 4:04 PM, Javi Hidalgo wrote: >

[R] Loop stopping after 1 iteration

Hi all, This is a very basic question, but I just can't figure out why R is handling a loop I'm writing the way it is. Here is the script I have written: grid_2_series<-function(gage_handle,data_type,filename) series_name<-paste(gage_handle,data_type,sep="_") data_grid<-read.table(file=paste(fi

Re: [R] Help with 2-D plot of k-mean clustering analysis

I wonder if it makes sense to reduce the dimensionality of the variables somehow? David Cross d.cr...@tcu.edu www.davidcross.us On May 18, 2011, at 9:41 AM, Meng Wu wrote: > Hi, all > > I would like to use R to perform k-means clustering on my data which > included 33 samples measured with

Re: [R] Integral Symbol

Thanks! This is what I was looking for. Apparently, it is not supported then it should be as "integral", but in the pdf version appears the integral symbol. Regards, Javier. > Date: Wed, 18 May 2011 09:39:28 -0400 > From: murdoch.dun...@gmail.com > To: havyhida...@hotmail.com > CC: r-help@r-pr

[R] retrieving gbif data

Hi, I am trying to use the "gbif" function in the "dismo" package and it does not seem to work. I get the same error message every time: Error in if (sp) geo <- TRUE : argument is not interpretable as logical It does not matter whether I query for my species of interest or if I copy and paste tho

Re: [R] R example code of Split-plot Manova

Hi, I'm a PhD student in Milan (Italy). I read the OBrienKaiser example in ?Anova in the car package. I think that this is not a Manova split-plot design. I need to know if someone knows a R code for MANOVA split-plot. Is there someone who can help me? Thanks for kindness Riccardo -- View this m

[R] Convolution confusion:

Hi, I'm new to R, and I'm a bit confused with the "convolve()" function. If I do: x<-c(1, 2, 3) convolve(x, rev(x), TRUE, "open") = 9 12 10 4 1 But I expected: 3 8 14 8 3 (like in Octave/MATLAB - conv(x, reverse(x)) ) 3 2 1 x 1 2 3 = 3 2 1 0 6 4 2 0 0 9 6 3 = 3 8 14 8 3 The thing is, t

Re: [R] Change pattern in histograms in ggplot2

There are a number of discussion threads on the google groups ggplot2 page: here are two of them. http://groups.google.com/group/ggplot2/browse_thread/thread/ca546f7f4d636deb/e0763a54b7735c35?lnk=gst&q=fill+pattern#e0763a54b7735c35 http://groups.google.com/group/ggplot2/browse_thread/thread/9a9c

Re: [R] Changing order of facet grid in ggplot2

data$tp <- factor(data$tp, levels = c("TP-ANY","TP-SUB","TP-CLIN")) qplot(life,geom="bar",weight=sym,ylim=c(0,1),legend=F,data=data) + facet_grid(. ~ tp) On Wednesday, May 18, 2011 at 9:14 AM, Christopher Desjardins wrote: Hi I am running the following code: > > sym <- c(sym1,sym2,sym4) > > life

[R] Help with 2-D plot of k-mean clustering analysis

Hi, all I would like to use R to perform k-means clustering on my data which included 33 samples measured with ~1000 variables. I have already used kmeans package for this analysis, and showed that there are 4 clusters in my data. However, it's really difficult to plot this cluster in 2-D format

Re: [R] subsetting a list of dataframes

> -Original Message- > From: r-help-boun...@r-project.org > [mailto:r-help-boun...@r-project.org] On Behalf Of Lara Poplarski > Sent: Tuesday, May 17, 2011 4:14 PM > To: r-help@r-project.org > Subject: Re: [R] subsetting a list of dataframes > > Thank you all, this is exactly what I had

[R] Changing order of facet grid in ggplot2

Hi I am running the following code: sym <- c(sym1,sym2,sym4) lifedxm <- c("O-BD","O-WELL","O-UNI") life <- c(lifedxm,lifedxm,lifedxm) tp <- c("TP-ANY","TP-ANY", "TP-ANY", "TP-SUB", "TP-SUB", "TP-SUB", "TP-CLIN" , "TP-CLIN", "TP-CLIN") data <- data.frame(sym,life,tp) qplot(life,geom="bar",w

Re: [R] R Style Guide -- Was Post-hoc tests in MASS using glm.nb

This is the first time I've seen an R Style Guide. I will admit that I haven't looked for one previously, but nevertheless I still haven't seen one. My code style simply evolved (perhaps, chugged along) by reading posts from other users who post to the r-help community. I regularly program with a

[R] R Style Guide -- Was Post-hoc tests in MASS using glm.nb

Thanks Bill. Do you and others think that a link to this guide (or another)should be included in the Posting Guide and/or R FAQ? -- Bert On Tue, May 17, 2011 at 4:07 PM, wrote: > Amen to all of that, Bert.  Nicely put.  The google style guide (not perfect, > but a thoughtful contribution on th

Re: [R] Integral Symbol

On 18/05/2011 9:09 AM, Javi Hidalgo wrote: Thanks. I was exactly reading the manual Writing R Extensions, on section Mathematics. Where, it informs about basic LaTeX style support. However, It seems like it does not support the LaTeX integral symbol \int, but it does support i.e.: the summation

[R] Change pattern in histograms in ggplot2

Hi, I am wondering if there is a way to change the pattern of the fill in histogram in ggplot2? By default the fill is solid and I'd like to add some sort of pattern to make it more visible that these are different levels of a factor. Thanks! Chris [[alternative HTML version deleted]] __

[R] matrix help (first occurrence of variable in column)

Dear R help, Apologies for the less than informative subject line. I will do my best to describe my problem. Consider the following matrix: mdat <- matrix(c(1,0,1,1,1,0), nrow = 2, ncol=3, byrow=TRUE, dimnames = list(c("T1", "T2"), c("sp.1", "sp.2", "

Re: [R] Smooth contour of a map

You may be looking for the par settings of xaxs="i", yaxs="i", which if you add them to the plot call will prevent the "regular" behavior of adding 4% padding to the axis widths. ?par -- David. On May 18, 2011, at 8:27 AM, Pierre Bruyer wrote: I've pratically resolved my problem (the code

Re: [R] subsetting a list of dataframes

On May 17, 2011, at 7:13 PM, Lara Poplarski wrote: Thank you all, this is exactly what I had in mind, except that I still have to get my head around apply et al. Back to the books for me then! Read the lapply( ...) call as: "For every element in the object named `data`, send that element

  1   2   >