Re: [R] help with sockets in R

2010-09-21 Thread Prof Brian Ripley
readLines() is for a text-mode connection; readChar() is for a binary-mode connection. Given that you asked for possible re-encoding by the 'encoding' argument, you cannot safely mix them (text-mode access is re-encoded, binary-mode is not). However, we don't know if re-encoding was active in

Re: [R] Web forum - should I make one?

2010-09-21 Thread Vojtěch Zeisek
Hello Dne St 22. září 2010 00:30:26 Yihui Xie napsal(a): > Ahh... I have almost forgotten the forum links I put in the R Wiki > four years ago! Thanks for reminding me! > > As many Chinese users know, the Chinese R forum at http://cos.name/cn/ > is pretty mature now. From my experience in maintai

Re: [R] Web forum - should I make one?

2010-09-21 Thread Yihui Xie
Ahh... I have almost forgotten the forum links I put in the R Wiki four years ago! Thanks for reminding me! As many Chinese users know, the Chinese R forum at http://cos.name/cn/ is pretty mature now. From my experience in maintaining this forum, I can fully understand Vojtěch's proposal, and in f

Re: [R] cwhmisc package error

2010-09-21 Thread Dennis Murphy
Hi: It's not hard to roll your own here. For first character capitalization, mycap <- function(str) { x <- substring(toupper(str), 1, 1) y <- substring(tolower(str), 2, nchar(str)) paste(x, y, sep = '') } It's not especially fast but it works if the objective is to capitalize only the fir

[R] randomForest - partialPlot - Reg

2010-09-21 Thread Vijayan Padmanabhan
Dear R Group I had an observation that in some cases, when I use the randomForest model to create partialPlot in R using the package "randomForest" the y-axis displays values that are more than -1! It is a classification problem that i was trying to address. Any insights as to how the y axis can

Re: [R] How to fix error: 'x' and 'y' lengths differ

2010-09-21 Thread David Winsemius
On Sep 21, 2010, at 11:27 PM, klsk89 wrote: Hi, im a student so still very new to R. Hope someone could help me out here =) They are 3 slug control products, bustaslug, product X and Y. Im ask to explore the data by plot() and tapply(). But when i try to plot or use the tapply command,

Re: [R] How to fix error: 'x' and 'y' lengths differ

2010-09-21 Thread Joshua Wiley
Hi, The problem has to do with names. You have a dataset called "slugs", inside this dataset is a column also called "slugs". Now, you attach() the slugs dataset, which adds its contents to the search path, but it adds it *after* the global environment (your workspace) where the dataset "slugs"

[R] How to fix error: 'x' and 'y' lengths differ

2010-09-21 Thread klsk89
Hi, im a student so still very new to R. Hope someone could help me out here =) They are 3 slug control products, bustaslug, product X and Y. Im ask to explore the data by plot() and tapply(). But when i try to plot or use the tapply command, it tells me that the x and y lengths differ. so im thin

Re: [R] cwhmisc package error

2010-09-21 Thread David Winsemius
On Sep 21, 2010, at 10:11 PM, Adrienne Wootten wrote: R-listers, I'm working on a project where I need to get the lowercase of the county variable in a dataset alike to the example below (only difference is that the full dataset has all the states and counties in the southeast United St

[R] cwhmisc package error

2010-09-21 Thread Adrienne Wootten
R-listers, I'm working on a project where I need to get the lowercase of the county variable in a dataset alike to the example below (only difference is that the full dataset has all the states and counties in the southeast United States). I keep getting this strange error with the lowerize funct

Re: [R] reshape is re-ordering my variables

2010-09-21 Thread Dennis Murphy
Hi: Reshaping multiple variables is nontrivial. Try the following (untested): reshape(rcw, idvar = 'ICU', varying = list(c(paste('Q6.RC', 1:4, sep = '.'), c(paste('Q6.FT.RC', 1:4, 'years', sep = '.'), c(paste('Q6.FT.RC', 1:4, 'months',

Re: [R] Package for calculating bandwidths

2010-09-21 Thread Dennis Murphy
Try this: # install.packages('sos')[if necessary] library(sos) findFn('bandwidth cross validation kernel smoothing') found 24 matches; retrieving 2 pages The ks package appears to be a place to start... HTH, Dennis On Tue, Sep 21, 2010 at 9:31 AM, Brocker84 wrote: > > No idea? > -- > Vie

Re: [R] Error in eval(expr, envir, enclos)

2010-09-21 Thread Dennis Murphy
If this is the same error from your initial post, try the following: (i) ls()# look for an object named Regression (ii) if ((i) is TRUE) str(Regession) Assuming you're getting the same error you reported in your initial post, either the data o

[R] unsubscribe [SEC=UNCLASSIFIED]

2010-09-21 Thread Story, Paul
Paul Story Environmental Officer Australian Plague Locust Commission Department of Agriculture, Fisheries and Forestry GPO Box 858 Canberra ACT 2601 Phone: (02) 6272-5176 Fax (02) 6272-5074 Mobile 0429-469306 Mail to: paul.st...@daff.gov.au mailto:paul.st...@daff.gov.au> or pg...

[R] help with sockets in R

2010-09-21 Thread Christopher Bare
Hi R gurus, I'm trying to use a ReSTful web service from within R. Specifically, I need to make HTTP PUT requests. I'm able to make the request and that goes well enough, but I'm having trouble properly consuming the HTTP response. The problem comes in when I'm trying to parse out the response bo

Re: [R] Size of the legend

2010-09-21 Thread Dennis Murphy
Hi: It is not good graphics for a legend to overlap data points, so let's see how one can get said legend away from the data. One approach in base graphics is to extend the y-axis limits and to place the legend in the open area. It would probably be better to split the legend into two halves and

[R] How to convert ARMA process to infinite AR?

2010-09-21 Thread Juan Pablo Calle
Hi, I need a function to convert an ARMA process to an "infinite" AR process. I know that in the Stats package exist one to convert an ARMA to MA, but i can't find one function to convert to AR. There exists one? THANKS __ R-help@r-project.org maili

Re: [R] Doing operations by grouping variable

2010-09-21 Thread Seth W Bigelow
Aah, that is the sort of truly elegant solution I have been seeking. And it's wrapped up in a nice programming shortcut to boot (i.e., the within statement). I retract anything I may have said about tapply being clunky. Many thanks --Seth Dr. Seth W. Bigelow Biologist, USDA-FS Pacific Southwe

Re: [R] How to convert a character into a filename?

2010-09-21 Thread Nordlund, Dan (DSHS/RDA)
> -Original Message- > From: r-help-boun...@r-project.org [mailto:r-help-boun...@r- > project.org] On Behalf Of Tucson August > Sent: Tuesday, September 21, 2010 3:15 PM > To: r-help@r-project.org > Subject: [R] How to convert a character into a filename? > > Dear list, > > How to convert

Re: [R] How to convert a character into a filename?

2010-09-21 Thread Tucson August
On Tue, Sep 21, 2010 at 3:39 PM, David Winsemius wrote: > > On Sep 21, 2010, at 6:14 PM, Tucson August wrote: > > Dear list, >> >> How to convert a character to a filename? >> > > Filenames are character mode so that shouldn't be a problem. > > > such as: >> >> x <- "height" # "height" here is

Re: [R] group means of multi-way table?

2010-09-21 Thread David Winsemius
On Sep 21, 2010, at 6:52 PM, Simon Kiss wrote: hello, can someone tell me how to generate the means for a data frame that looks like this? My data frame has many more variables, but I won't bother you with those; these are the one's that I'm interested in. Needless to say, z is the variabl

Re: [R] How to convert a character into a filename?

2010-09-21 Thread Tucson August
that works! Thanks a lot!!! On Tue, Sep 21, 2010 at 3:29 PM, stephen sefick wrote: > file <- foo.data > > filename <- paste("plant", x, sep="") > out <- paste("C:/", "plant/", filename, ".csv",sep="") > > write.csv(file, out) > > On Tue, Sep 21, 2010 at 5:14 PM, Tucson August > wrote: > > D

[R] group means of multi-way table?

2010-09-21 Thread Simon Kiss
hello, can someone tell me how to generate the means for a data frame that looks like this? My data frame has many more variables, but I won't bother you with those; these are the one's that I'm interested in. Needless to say, z is the variable in which I'm interested. I'd like to find out th

Re: [R] How to convert a character into a filename?

2010-09-21 Thread David Winsemius
On Sep 21, 2010, at 6:14 PM, Tucson August wrote: Dear list, How to convert a character to a filename? Filenames are character mode so that shouldn't be a problem. such as: x <- "height" # "height" here is actually a name of a colume in a data frame So you need to learn how to refe

Re: [R] How to convert a character into a filename?

2010-09-21 Thread stephen sefick
file <- foo.data filename <- paste("plant", x, sep="") out <- paste("C:/", "plant/", filename, ".csv",sep="") write.csv(file, out) On Tue, Sep 21, 2010 at 5:14 PM, Tucson August wrote: > Dear list, > > How to convert a character to a filename? > such as: > > x <- "height"  # "height" here is ac

Re: [R] lmer() vs. lme() gave different variance component estimates

2010-09-21 Thread Douglas Bates
On Tue, Sep 21, 2010 at 4:21 PM, Peter Dalgaard wrote: > On 09/21/2010 09:02 PM, Douglas Bates wrote: >> On Tue, Sep 21, 2010 at 1:39 PM, Douglas Bates wrote: >>> I haven't had the time to keep up with this discussion, or many of the >>> other discussions on the R-SIG-Mixed-Models email list.  I

Re: [R] Doing operations by grouping variable

2010-09-21 Thread Bill.Venables
You left out the subscript. Why not just do d <- within(data.frame(group = rep(1:5, each = 5), variable = rnorm(25)), scaled <- variable/tapply(variable, group, max)[group]) and be done with it? (Warning: if you replace the second '<-' above by '=', it will not work. It is NOT true t

[R] How to convert a character into a filename?

2010-09-21 Thread Tucson August
Dear list, How to convert a character to a filename? such as: x <- "height" # "height" here is actually a name of a colume in a data frame filename <- paste("plant,x") write.csv (data, file="C:/plant/filename.csv) # having trouble with this statement, how to 'write' the filename here? All I w

[R] Fwd: mode across lists of matrices

2010-09-21 Thread Gregory Ryslik
Hi Everyone, Looks like I found the solution, the mode function is now changed to: mode <- function(x){ cell.mode<-as.numeric(names(which.max(table(x if(length(cell.mode)==0){ cell.mode <- NA } cell.mode } Then you can ta

Re: [R] Need help for EM algorithm ASAP !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

2010-09-21 Thread Carl Witthoft
The problem is clear: you failed to follow proper intarweb protocol in your Subject: line. It should read [R] Need help for EM algorithm ASAP !!!111 ONE ONE ONE LOL __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PL

Re: [R] Prime Factorization

2010-09-21 Thread David Winsemius
On Sep 21, 2010, at 5:31 PM, Hans W Borchers wrote: David Winsemius comcast.net> writes: A further citation that answers the question I raised (and inaccurately predicted no value) regarding prime.sieve : http://finzi.psych.upenn.edu/R/Rhelp02/archive/49773.html This was found with Barons

Re: [R] Prime Factorization

2010-09-21 Thread Hans W Borchers
David Winsemius comcast.net> writes: > > A further citation that answers the question I raised (and > inaccurately predicted no value) regarding prime.sieve : > > http://finzi.psych.upenn.edu/R/Rhelp02/archive/49773.html > > This was found with Barons search facility set for rhelp postings:

Re: [R] removed data is still there!

2010-09-21 Thread Greg Snow
This comes up every now and then. The fact is that the behavior of R in not throwing away information unless explicitly told to, is a feature, and one that I don't want to see go away. Yes in your example doing a table or plot based on iris1$Species gives meaningless results, but anything you

Re: [R] lmer() vs. lme() gave different variance component estimates

2010-09-21 Thread Peter Dalgaard
On 09/21/2010 09:02 PM, Douglas Bates wrote: > On Tue, Sep 21, 2010 at 1:39 PM, Douglas Bates wrote: >> I haven't had the time to keep up with this discussion, or many of the >> other discussions on the R-SIG-Mixed-Models email list. I swamped >> with other duties at present. It's not like I don

Re: [R] multiplying values in data frame by corresponding value in the first column

2010-09-21 Thread Martin Tomko
Thank you Marc, that was an elegant solution, works perfectly! Martin On 9/21/2010 10:24 PM, Marc Schwartz wrote: DF[, -c(1:2)]<- DF[, -c(1:2)] * DF$Group __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do re

Re: [R] Prime Factorization

2010-09-21 Thread David Winsemius
A further citation that answers the question I raised (and inaccurately predicted no value) regarding prime.sieve : http://finzi.psych.upenn.edu/R/Rhelp02/archive/49773.html This was found with Barons search facility set for rhelp postings: http://search.r-project.org/cgi-bin/namazu.cgi?query=

Re: [R] Indexing sublists inside lists.

2010-09-21 Thread Dennis Murphy
Hi: Is this what you were looking for? frd <- rep('fred', 5) bigfred <- lapply(frd, get) On Tue, Sep 21, 2010 at 5:04 AM, Alaios wrote: > I would like to thank you very much for your reply. > Actually I would like to ask you if there is > a small list called fred: > fred <- list(happy = 1:10,

Re: [R] bivariate vector numerical integration with infinite range

2010-09-21 Thread baptiste auguie
Got it, thanks! baptiste On 21 September 2010 22:38, Hans W Borchers wrote: > baptiste auguie googlemail.com> writes: > >> >> Thanks. I am having trouble getting adaptIntegrate to work with a >> multivalued integrand though, and cannot find a working example. >> Anyone had better luck with it?

Re: [R] R-help Digest, Vol 91, Issue 21

2010-09-21 Thread Marc Schwartz
On Sep 21, 2010, at 2:47 PM, kurt_h...@nps.gov wrote: > All > Is there a script in R equivalent to the "if then" transforms one can > perform in Systat? For example, I want to create a "Treatment" column > coded either 1 or 2 for twelve field sites in a large data set. Ideally, > I'd be able

Re: [R] R-help Digest, Vol 91, Issue 21

2010-09-21 Thread Peter Dalgaard
On 09/21/2010 09:47 PM, kurt_h...@nps.gov wrote: > All > Is there a script in R equivalent to the "if then" transforms one can > perform in Systat? For example, I want to create a "Treatment" column > coded either 1 or 2 for twelve field sites in a large data set. Ideally, > I'd be able to t

Re: [R] R-help Digest, Vol 91, Issue 21

2010-09-21 Thread David Winsemius
On Sep 21, 2010, at 3:47 PM, kurt_h...@nps.gov wrote: All Is there a script in R equivalent to the "if then" transforms one can perform in Systat? For example, I want to create a "Treatment" column coded either 1 or 2 for twelve field sites in a large data set. Ideally, I'd be able

Re: [R] bivariate vector numerical integration with infinite range

2010-09-21 Thread Hans W Borchers
baptiste auguie googlemail.com> writes: > > Thanks. I am having trouble getting adaptIntegrate to work with a > multivalued integrand though, and cannot find a working example. > Anyone had better luck with it? The function to be integrated needs a vector as input: f <- function(x) {

Re: [R] multiplying values in data frame by corresponding value in the first column

2010-09-21 Thread Marc Schwartz
On Sep 21, 2010, at 3:00 PM, Martin Tomko wrote: > I am sure there is a simple solution to this... I have a column in a data > frame specifying a grouping (1, -1) for my observations, and need to mutliply > each observation in all the other columns of the data frame by the > corresponding value

Re: [R] Prime Factorization

2010-09-21 Thread David Winsemius
On Sep 21, 2010, at 3:44 PM, Spencer Graves wrote: Hello: It's hard to RTFM if you don't know which FMTR. I only found one with the "obvious search strategy" with my favorite search tools, which may nevertheless be unfamiliar to Cliff Clive. RSiteSearch('prime factorization') # pro

Re: [R] Lattice xyplot and groups

2010-09-21 Thread Axel
> Does using > >     df = df[order(df$type,df$set,df$x),] > > before calling xyplot fix the problem? Thank you very much for your suggestion. It does make the lines open, but unfortunately it mixes different id together. With the new ordering the plot looks something like this: (x=-10, y_id1) ->

Re: [R] bivariate vector numerical integration with infinite range

2010-09-21 Thread baptiste auguie
Thanks. I am having trouble getting adaptIntegrate to work with a multivalued integrand though, and cannot find a working example. Anyone had better luck with it? library(cubature) > > f <- function(x, y) { + res <- 1 / (sqrt(x)*(1+x)) + c(res, res/2, 2*res) + } > > adaptIntegrate(f, lowerLimi

Re: [R] mode across lists of matrices

2010-09-21 Thread Gregory Ryslik
Hey, Looks like I spoke to soon last time. It's almost there but not quite. Now If there are more NA's than any other number, it returns NA. That's not the correct behavior though as I want to return NA only if all the numbers are NA. For example, in the situation below, I want it to return 5 f

[R] multiplying values in data frame by corresponding value in the first column

2010-09-21 Thread Martin Tomko
I am sure there is a simple solution to this... I have a column in a data frame specifying a grouping (1, -1) for my observations, and need to mutliply each observation in all the other columns of the data frame by the corresponding value in the given column. I played with apply, and saw some s

Re: [R] R-help Digest, Vol 91, Issue 21

2010-09-21 Thread Kurt_Helf
All Is there a script in R equivalent to the "if then" transforms one can perform in Systat? For example, I want to create a "Treatment" column coded either 1 or 2 for twelve field sites in a large data set. Ideally, I'd be able to tell R to code sites a-f as 1 and sites g-l as 2. Cheers Kur

Re: [R] package gbm, predict.gbm with offset

2010-09-21 Thread Ridgeway, Greg
You've got it right. Just multiply the predicted values by the exposure time or number of cases and that will get you what you want. Note that, depending on the scale of the predictions, you might want to do all the predictions and +log(offset) on the log scale and then exponentiate it. Should

Re: [R] Prime Factorization

2010-09-21 Thread Spencer Graves
Hello: It's hard to RTFM if you don't know which FMTR. I only found one with the "obvious search strategy" with my favorite search tools, which may nevertheless be unfamiliar to Cliff Clive. RSiteSearch('prime factorization') # produced 3 matches, the first of which is the "primeFa

Re: [R] Prime Factorization

2010-09-21 Thread Hans W Borchers
David Winsemius comcast.net> writes: > > On Sep 21, 2010, at 11:33 AM, Cliff Clive wrote: > > > > > Hi everyone, I have a very quick question: > > > > Is there a ready-made function in R or any R packages to find the > > prime > > factorization of an integer? > > Yes. At least two. The obvio

Re: [R] mode across lists of matrices

2010-09-21 Thread Gregory Ryslik
Perfect! Thank you! On Sep 21, 2010, at 3:27 PM, Henrique Dallazuanna wrote: > Use this function: > > mode <- function(x, ...) > as.numeric(names(which.max(table(x, useNA = 'ifany' > > On Tue, Sep 21, 2010 at 4:18 PM, Gregory Ryslik wrote: > Hi, > > Much better! That seems to work gr

Re: [R] mode across lists of matrices

2010-09-21 Thread Henrique Dallazuanna
Use this function: mode <- function(x, ...) as.numeric(names(which.max(table(x, useNA = 'ifany' On Tue, Sep 21, 2010 at 4:18 PM, Gregory Ryslik wrote: > Hi, > > Much better! That seems to work great. The only time it doesn't work is > that if all the elements are NA for a specific eleme

Re: [R] mode across lists of matrices

2010-09-21 Thread Gregory Ryslik
Hi, Much better! That seems to work great. The only time it doesn't work is that if all the elements are NA for a specific element. Then it gives me a numeric(0) for that cell rather than NA. On Sep 21, 2010, at 2:30 PM, Henrique Dallazuanna wrote: > Try this: > > apply(array(unlist(mymats),

Re: [R] missing package tensorA

2010-09-21 Thread Uwe Ligges
Hi, can you please try the prerelease R-2.12.0 alpha? It comes with both 32- and 64-bit support. tensorA does not pass checks under R-2.11.1 patched 64-bit these days and is currently not available. We will stop support for R-2.11.1 64-bit (not 32-bit!) shortly after the R-2.12.0 release anywa

Re: [R] Lattice xyplot and groups

2010-09-21 Thread Phil Spector
Does using df = df[order(df$type,df$set,df$x),] before calling xyplot fix the problem? - Phil Spector Statistical Computing Facility Department of Statistics

Re: [R] diagnosing download.file() problems

2010-09-21 Thread steven mosher
That is what I feared. I know other people on slow connections have done this without issue ( at least they didnt report an issue ) I had a similar issue with geonames.org who at least published their terms of service (requests per second or sumptin like that.. so could program around it) I'll hu

Re: [R] Trouble installing pwr package

2010-09-21 Thread Gavin Simpson
On Tue, 2010-09-21 at 12:04 -0500, Erik Iverson wrote: > > Erik Iverson wrote: > > > > > > Brian J Mingus wrote: > >> Hi all, > >> > >> I'm having trouble getting access to the pwr. This is on Ubuntu Lucid > >> Lynx, > >> 64 bit. > >> > >> I'm installing pwr via packages.install('pwr'), and loa

Re: [R] lmer() vs. lme() gave different variance component estimates

2010-09-21 Thread Douglas Bates
On Tue, Sep 21, 2010 at 1:39 PM, Douglas Bates wrote: > I haven't had the time to keep up with this discussion, or many of the > other discussions on the R-SIG-Mixed-Models email list.  I swamped > with other duties at present. > > It is important to remember that the nlme and lme4 packages take a

[R] reshape is re-ordering my variables

2010-09-21 Thread Kevin E. Thorpe
Is it an undocumented (at least I missed it if it's documented) feature of the reshape function to do numeric variables followed by character? I ask because that seems to be the case below. > str(rcw) 'data.frame': 23 obs. of 21 variables: $ ICU : int 1 18 17 9 22 19 6 16 25 26

Re: [R] Package for calculating bandwidths

2010-09-21 Thread David Winsemius
On Sep 21, 2010, at 12:31 PM, Brocker84 wrote: No idea? Search? Perhaps with: "adaptive bandwidth cross-validation" AND (for your brain, not the search engine) -- PLEASE do read the posting guide http://www.R-project.org/posting-guide.html --- David Winsemius, MD West Ha

Re: [R] Lists with NULL entries

2010-09-21 Thread Greg Snow
Why not just do: lst2 <- c(list(NULL), lst) -- Gregory (Greg) L. Snow Ph.D. Statistical Data Center Intermountain Healthcare greg.s...@imail.org 801.408.8111 > -Original Message- > From: r-help-boun...@r-project.org [mailto:r-help-boun...@r- > project.org] On Behalf Of Peter Langfeld

[R] Lattice xyplot and groups

2010-09-21 Thread Axel
Hi, I'm trying to plot many (x, y) data files using the xyplot function from the lattice package. Each file can be classified by set name (s1, s2,...) and data type (A, B, ...). Each data set contains a different number of files. If the data is grouped by type or set and visualized as line plot wi

Re: [R] Please Help_Error:cannot allocate vector of size 400.4Mb

2010-09-21 Thread Peter Langfelder
On Tue, Sep 21, 2010 at 9:50 AM, qcshare wrote: > > Hello, everyone, > When I run R, I met: > "error:cannot allocate vector of size 400Mb", My data is large. > What should I do? > Thanks, everyone. > How big is the RAM in your computer? There are a few things you can try: 1. Before running the

Re: [R] lmer() vs. lme() gave different variance component estimates

2010-09-21 Thread Douglas Bates
I haven't had the time to keep up with this discussion, or many of the other discussions on the R-SIG-Mixed-Models email list. I swamped with other duties at present. It is important to remember that the nlme and lme4 packages take a model specification and provide code to evaluate the deviance.

Re: [R] mode across lists of matrices

2010-09-21 Thread Henrique Dallazuanna
Try this: apply(array(unlist(mymats), dim = c(dim(mymats[[1]]), length(mymats))), 1:2, mode) The error was in c(length(mymats), dim(mymats[[1]])) On Tue, Sep 21, 2010 at 3:19 PM, Gregory Ryslik wrote: > Ack, apologies for the previous email. What I meant to say is that the > first element is

[R] Regarding Contour maps

2010-09-21 Thread syama adhikarla
Hello, I would like to know how to go about creating contour maps (Isofrequency maps) in R using Kriging procedure. Thank You Regards Syama [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailma

[R] missing package tensorA

2010-09-21 Thread Gabriela Cendoya
Hi: I was trying to download the package MCMCglmm and it give me this message: Aviso: dependency ‘tensorA’ is not available probando la URL 'http://www.vps.fmvz.usp.br/CRAN/bin/windows64/contrib/2.11/MCMCglmm_2.06.zip' Content type 'application/zip' length 8988896 bytes (8.6 Mb) URL abierta downlo

Re: [R] Colorramp in Maptools, how to choose min and max values for the fg= argument

2010-09-21 Thread Jannis
Hi Katrin, your problem is not related to maptools, rather to the base function "colorRampPalette". In my opinion the easiest way to get what you want is to use colorRamp: Vector=c(2,3,4,5,6) color.ramp <- colorRamp(c('white','black')) map.colors <- color.ramp(Vector/10) The 10 refers to the

Re: [R] mode across lists of matrices

2010-09-21 Thread Gregory Ryslik
Ack, apologies for the previous email. What I meant to say is that the first element is calculated incorrectly (1,1) should be 2, instead it is 3. I've been staring at the code for two long. I've copied it in again for convenience. mymats <- vector('list', 4) set.seed(246) # Generate a list of

[R] Please Help_Error:cannot allocate vector of size 400.4Mb

2010-09-21 Thread qcshare
Hello, everyone, When I run R, I met: "error:cannot allocate vector of size 400Mb", My data is large. What should I do? Thanks, everyone. -- View this message in context: http://r.789695.n4.nabble.com/Please-Help-Error-cannot-allocate-vector-of-size-400-4Mb-tp2549004p2549004.html Sent from the

[R] package gbm, predict.gbm with offset

2010-09-21 Thread Markus Loecher
Dear all, the help file for predict.gbm states that "The predictions from gbm do not include the offset term. The user may add the value of the offset to the predicted value if desired." I am just not sure how exactly, especially for a Poisson model, where I believe the offset is multiplicative ?

Re: [R] HOW to create image like this?

2010-09-21 Thread zcrself
I like your answer! -- View this message in context: http://r.789695.n4.nabble.com/HOW-to-create-image-like-this-tp2548152p2548945.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org mailing list https://stat.ethz.ch/

Re: [R] Prime Factorization

2010-09-21 Thread David Winsemius
On Sep 21, 2010, at 11:33 AM, Cliff Clive wrote: Hi everyone, I have a very quick question: Is there a ready-made function in R or any R packages to find the prime factorization of an integer? Yes. At least two. The obvious search strategy with your favrite search tool should work well

Re: [R] function help?

2010-09-21 Thread Cormac Long
Hi there Petr, Apologies for only replying to your post now - I hope the code included below helps you out. An alternative function documentation approach would be the one I took when faced with a similar problem - I wrote some functions to allow me to embed a retrievable function specification i

Re: [R] mode across lists of matrices

2010-09-21 Thread Gregory Ryslik
Hi, That doesn't work as seen in this example. I get a 4x3 matrix as the answer which is definitely not right! But I think we're getting closer! Cheers, G for(i in 1:4) mymats[[i]] <- matrix(sample(1:9), nrow = 3) mymats[[1]][1,1]<-3 mymats[[1]][1,2]<-3 mymats[[1]][1,3]<-1 mymats[[2]][2,1]<-

Re: [R] Package for calculating bandwidths

2010-09-21 Thread Brocker84
No idea? -- View this message in context: http://r.789695.n4.nabble.com/Package-for-calculating-bandwidths-tp2548091p2548976.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org mailing list https://stat.ethz.ch/mailma

Re: [R] Survival curve mean adjusted for covariate: NEED TO DO IN NEXT 2 HOURS, PLEASE HELP

2010-09-21 Thread mark.fisher123
Thanks for the help everyone - problem solved in time! -- View this message in context: http://r.789695.n4.nabble.com/Survival-curve-mean-adjusted-for-covariate-NEED-TO-DO-IN-NEXT-2-HOURS-PLEASE-HELP-tp2548484p2548854.html Sent from the R help mailing list archive at Nabble.com. ___

Re: [R] Error in eval(expr, envir, enclos)

2010-09-21 Thread uttara_n
Hello Joshua, It still gives me the same error. Uttara -- View this message in context: http://r.789695.n4.nabble.com/Error-in-eval-expr-envir-enclos-tp2547917p2548861.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project

[R] Prime Factorization

2010-09-21 Thread Cliff Clive
Hi everyone, I have a very quick question: Is there a ready-made function in R or any R packages to find the prime factorization of an integer? -- View this message in context: http://r.789695.n4.nabble.com/Prime-Factorization-tp2548877p2548877.html Sent from the R help mailing list archive at

Re: [R] Web forum - should I make one?

2010-09-21 Thread chockri adnen
Hello, it's good idea, i'm with you in this project. I know very well Joomla, VB forum, IPB forum, PHP forum, SMF forum... Best regards, Adnen CHOCKRI Department of Quantitative Methods , University of Economics and Management, Sfax, Tunisia Community of the Ubuntu GNU/Linux 2010/9/21 Vojtě

Re: [R] Web forum - should I make one?

2010-09-21 Thread Vojtěch Zeisek
Hello Dne Út 21. září 2010 11:55:02 Marc Schwartz napsal(a): > Before creating yet another R discussion venue, you may wish to look at: > > http://stackoverflow.com/questions/tagged/r My idea was about something like this, but more complex, with categories and tools. And only for R. :-) This

Re: [R] Size of the legend

2010-09-21 Thread David Winsemius
On Sep 21, 2010, at 1:28 PM, Mestat wrote: Hi Denis, Check it out my code... This is not my real data... I would like to manage the size of the legend... Set the legend smaller than it is, because on my real data, the legend is over the values... Thanks for the help... x<-runif(100) y<-run

Re: [R] Web forum - should I make one?

2010-09-21 Thread Erik Iverson
Everyone is free to create as they please. My opinion is that between this list and stackexchange (as Tal pointed out), that there isn't a market for such a thing. If you want a web front-end to this list, such things exist, like nabble. Vojtěch Zeisek wrote: Hello, this might be little off-to

Re: [R] puzzle with integrate over infinite range

2010-09-21 Thread Duncan Murdoch
On 21/09/2010 1:29 PM, baptiste auguie wrote: I see, thank you. I'm still worried by the very dramatic error I obtained just from shifting so slightly the support of the integrand, it took me a while to figure what happened even with this basic example (I knew the integral couldn't be so small!

Re: [R] Web forum - should I make one?

2010-09-21 Thread Tal Galili
And of course there is also: http://stats.stackexchange.com/questions/tagged/r Cheers, T Contact Details:--- Contact me: tal.gal...@gmail.com | 972-52-7275845 Read me: www.talgalili.com (Hebrew) | www.biostatistics.co.il (Hebre

Re: [R] puzzle with integrate over infinite range

2010-09-21 Thread baptiste auguie
I see, thank you. I'm still worried by the very dramatic error I obtained just from shifting so slightly the support of the integrand, it took me a while to figure what happened even with this basic example (I knew the integral couldn't be so small!). For a general integration in [0, infty), ther

Re: [R] Size of the legend

2010-09-21 Thread Mestat
Hi Denis, Check it out my code... This is not my real data... I would like to manage the size of the legend... Set the legend smaller than it is, because on my real data, the legend is over the values... Thanks for the help... x<-runif(100) y<-runif(100) color<-c(rep(1,15),rep(2,20),rep(3,15),rep

Re: [R] Trouble installing pwr package

2010-09-21 Thread David Winsemius
On Sep 21, 2010, at 1:01 PM, Erik Iverson wrote: Brian J Mingus wrote: Hi all, I'm having trouble getting access to the pwr. This is on Ubuntu Lucid Lynx, 64 bit. I'm installing pwr via packages.install('pwr'), and loading it via library(pwr), both of which appear successful. Perhaps th

Re: [R] Trouble installing pwr package

2010-09-21 Thread Phil Spector
Brian - After you load a package with the library command, you can see what objects are available by looking at position 2 in the search path: library(pwr) objects(pos=2) [1] "cohen.ES" "ES.h" "ES.w1" "ES.w2" [5] "pwr.2p2n.test" "pwr.2p.test""pwr.anova.test

Re: [R] puzzle with integrate over infinite range

2010-09-21 Thread Thomas Lumley
On Tue, 21 Sep 2010, baptiste Auguié wrote: Thanks, I'll do that too from now on. It strikes me that in a case such as this one it may be safer to use a truncated, finite interval around the region where the integrand is non-zero, rather than following the advice of ?integrate to use Inf as in

Re: [R] Trouble installing pwr package

2010-09-21 Thread Brian J Mingus
On Tue, Sep 21, 2010 at 11:04 AM, Erik Iverson wrote: > > > Erik Iverson wrote: > >> >> >> Brian J Mingus wrote: >> >>> Hi all, >>> >>> I'm having trouble getting access to the pwr. This is on Ubuntu Lucid >>> Lynx, >>> 64 bit. >>> >>> I'm installing pwr via packages.install('pwr'), and loading i

Re: [R] Trouble installing pwr package

2010-09-21 Thread Erik Iverson
Erik Iverson wrote: Brian J Mingus wrote: Hi all, I'm having trouble getting access to the pwr. This is on Ubuntu Lucid Lynx, 64 bit. I'm installing pwr via packages.install('pwr'), and loading it via library(pwr), both of which appear successful. Perhaps the actual output would help h

Re: [R] partial dbRDA or CCA with two distance objects in Vegan.

2010-09-21 Thread Gavin Simpson
On Wed, 2010-09-22 at 02:34 +1000, Nevil Amos wrote: > Thanks for that, just to confirm I therefore need to use: > #if > #geogdist is a geographic distance matrix > #gen_dist is a genetic distance matrix > #env_var are environmental variables > > mypcnm<-pcnm(geogdist) > > > mydbRDA<-capscale(g

Re: [R] Trouble installing pwr package

2010-09-21 Thread Erik Iverson
Brian J Mingus wrote: Hi all, I'm having trouble getting access to the pwr. This is on Ubuntu Lucid Lynx, 64 bit. I'm installing pwr via packages.install('pwr'), and loading it via library(pwr), both of which appear successful. Perhaps the actual output would help here. Strangely, I neve

[R] Trouble installing pwr package

2010-09-21 Thread Brian J Mingus
Hi all, I'm having trouble getting access to the pwr. This is on Ubuntu Lucid Lynx, 64 bit. I'm installing pwr via packages.install('pwr'), and loading it via library(pwr), both of which appear successful. Strangely, I never get access to the pwr object in R. I tried installing it to /usr/local

Re: [R] Web forum - should I make one?

2010-09-21 Thread Marc Schwartz
Before creating yet another R discussion venue, you may wish to look at: http://stackoverflow.com/questions/tagged/r The R-* e-mail lists are the "official" venues and you can read/post via e-mail. There are also other means of interacting with the e-mail lists using Gmane and Nabble. Howeve

Re: [R] Can ucminf be installed in 64 bit R and one more question?

2010-09-21 Thread Duncan Murdoch
On 21/09/2010 11:43 AM, Hey Sky wrote: Hey, Duncan thanks for your reply. I am not sure which version i have installed but I downloaded it from http://cran.skazkaforyou.com/. when I check the R installed, it says 2.11.1. I do not know I answered your question or not. if not, where I can find

Re: [R] Doing operations by grouping variable

2010-09-21 Thread Seth W Bigelow
Bill, I'd never heard of ave( ) before, but it works like a charm for the purpose outlined below Thanks!! Dr. Seth W. Bigelow Biologist, USDA-FS Pacific Southwest Research Station 1731 Research Park Drive, Davis California "William Dunlap" 09/21/2010 08:52 AM To "Seth W Bigelow" , "R-

  1   2   3   >