[R] summary in functions

2011-02-16 Thread Sam Steingold
summary() in functions seems to print nothing. str() does print something. why? -- Sam Steingold (http://sds.podval.org/) on CentOS release 5.3 (Final) http://mideasttruth.com http://truepeace.org http://iris.org.il http://pmw.org.il http://jihadwatch.org http://ffii.org A year spent in artificial

[R] help with spatstat ppp transformation data

2011-02-16 Thread Ana María Prieto
Hello I am currently studying spatial distribution of lizards in Spain and Germany. For this I am using the spatstat package of R to make a Kolmogorov-Smirnov test of CSR of unmarked points data. I have problems transforming the data to ppp class, and I would like to know if you can help me. I h

Re: [R] A Math question

2011-02-16 Thread rex.dwyer
If y'all want to discuss this more, do it somewhere else, please. This has little to do with R except that both depend on Peano's Axioms. -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of jlu...@ria.buffalo.edu Sent: Tuesday, February

Re: [R] caret::train() and ctree()

2011-02-16 Thread Max Kuhn
Andrew, ctree only tunes over mincriterion and ctree2 tunes over maxdepth (while fixing mincriterion = 0). Seeing both listed as the function is being executed is a bug. I'll setup checks to make sure that the columns specified in tuneGrid are actually the tuning parameters that are used. Max O

Re: [R] leap years in temporal series command ts

2011-02-16 Thread Peter Zorglub
Ok thanks. I got now a temporal series of 2557 elements thanks to the below code: > seriezoo=zooreg(serie, start=as.Date("2002-01-01")) > seriezooTS=as.ts(seriezoo) The structure of seriezooTS is then the following one (str command): > str(seriezooTS) Time-Series [1:2557] from 11

[R] Function scope issues

2011-02-16 Thread Sébastien Bihorel
Dear R-users, I have some questions about the scope of functions, which I would like to illustrate with the following example: ### First example require(Hmisc) combine.levels <- function(x,y) browser() whatever <- function(x,y) combine.levels(x,y) x <- 1:5 names(x) <- LETTERS[x] y <- 6: 10 name

[R] confused by lapply

2011-02-16 Thread Sam Steingold
Description: 'lapply' returns a list of the same length as 'X', each element of which is the result of applying 'FUN' to the corresponding element of 'X'. I expect that when I do > lapply(vec,f) f would be called _once_ for each component of vec. this is not what I see: parse.num

[R] incomplete final line

2011-02-16 Thread mipplor
Hi, I work like this: > data<‐read.table('E:/my documents/r/consumer.xls',header=TRUE) Warning message: In read.table("E:/my documents/r/consumer.xls", header = TRUE) : incomplete final line found by readTableHeader on 'E:/my documents/r/consumer.xls' could someone shoot the trouble for me ?

Re: [R] string parsing

2011-02-16 Thread Sam Steingold
> * David Winsemius [2011-02-16 13:33:32 -0500]: > >> parse.num <- function (s) { >> as.numeric(gsub("M$","e6",gsub("B$","e9",s))); } > > data[1] <- parse.num( data[[1]] ) # as.numeric and gsub are vectorized because parse.num turned out to not be as simple as that: I need to handle "N/A" specia

[R] predict not complete?

2011-02-16 Thread cosmokramer
take this code library(dynlm) test=dynlm(diff(log(y))~ lag(x, -1) + lag(z, -1)) fc=predict(test) fc when I run this and look at "fc", the values there only go as long as the sample (of y) goes, but it should be one further, right? since the other values are lagged. Is predict the wrong functi

[R] create a data frame with the given column names

2011-02-16 Thread Sam Steingold
how do I create a data frame with the given column names _NOT KNOWN IN ADVANCE_? i.e., I have a vector of strings for names and I want to get an _EMPTY_ data frame with these column names. is it at all possible? -- Sam Steingold (http://sds.podval.org/) on CentOS release 5.3 (Final) http://openvo

Re: [R] Cannot allocate memory block

2011-02-16 Thread poisontonic
Uwe Ligges-3 wrote: > > If the available space got too fragmented, there is not single 3.8 block > of memory available any more > Is there anything I can do to prevent this? I've restarted and rerun the whole thing straight up, and still the error...? Ben -- View this message in context

[R] vector splicing

2011-02-16 Thread Sam Steingold
the following seems to work for numeric but not for times: parse.vec <- function (vec, parser = as.numeric, vmode = "numeric", na.string = "N/A") { bad <- (vec == na.string); ret <- vector(mode = vmode, length = length(vec)); if (some(bad)) ret[bad] <- NA; if (so

[R] calling pairs of variables into a function

2011-02-16 Thread squamous
Hi, I have imported three text files into R using read.table. Their variables are called d, e and f. I want to run a function on all the possible combinations of these three files. The only way I know how to do that is like this: bigfunction(d,e) bigfunction(d,f) bigfunction(e,d) bigfunction(

[R] 3 questions about the poisson regression of contingency table

2011-02-16 Thread Lao Meng
Hi all: I have 3 questions about the poisson regression of contingency table. Q1¡¢How to understand the "independent poisson process"as many books or paper mentioned? For instance: Table1 --- treat canernon-cancersum -

Re: [R] leap years in temporal series command ts

2011-02-16 Thread Gabor Grothendieck
On Wed, Feb 16, 2011 at 11:37 AM, Peter Zorglub wrote: > Ok thanks. I got now a temporal series of 2557 elements thanks to the below > code: > >> seriezoo=zooreg(serie, start=as.Date("2002-01-01")) >> seriezooTS=as.ts(seriezoo) > > The structure of seriezooTS is then the following one (str command

Re: [R] Timeseries Data Plotted as Monthly Boxplots

2011-02-16 Thread Thomas Adams
Katrina, What I have done, if I understand what you are after, was to create a list for each month of data - in order. Then, create a boxplot - in order - by month/year. I do this for our ensemble streamflow forecasts. The key us to create the list of values by month. Regards, Tom Sent from m

Re: [R] summary in functions

2011-02-16 Thread Simon Blomberg
On 17/02/11 09:44, Sam Steingold wrote: summary() in functions seems to print nothing. str() does print something. why? summary() returns the summary information as its value. If you want to see this value from inside a function, use print(summary()). The reason you see the summary informat

Re: [R] logsumexp function in R

2011-02-16 Thread Petr Savicky
On Wed, Feb 16, 2011 at 04:14:38PM -0500, Brian Tsai wrote: > Hi, > > I was wondering if anyone has implemented a numerically stable function to > compute log(sum(exp(x))) ? Try the following log(sum(exp(x - max(x + max(x) If this is not sufficient, consider using CRAN package Rmpfr with

Re: [R] incomplete final line

2011-02-16 Thread Joshua Wiley
Hi, On Wed, Feb 16, 2011 at 10:11 AM, mipplor wrote: > Hi, > > I work like this: > >> data<‐read.table('E:/my documents/r/consumer.xls',header=TRUE) > Warning message: > In read.table("E:/my documents/r/consumer.xls", header = TRUE) : >  incomplete final line found by readTableHeader on 'E:/my >

Re: [R] package ChemometricsWithR

2011-02-16 Thread Prof Brian Ripley
On Wed, 16 Feb 2011, deborahsills wrote: Did you find out where one can get the package ChemometricsWithR. All we *know* is that it is said to be with the book, but it has been submitted to CRAN. If accepted, it will appear there. Thank you! Deborah PLEASE do read the posting guide ht

Re: [R] Flipping axis 1 in CCA plot - vegan

2011-02-16 Thread Jari Oksanen
esgarbanzo uw.edu> writes: > > > Hello, > I've successfully run my CCA and the ordination is not displayed quite > right. Is it possible to flip the x-axis (ie., make it a mirror image of > itself)? > Howdy, I don't know wha tyou mean when you say that the ordination is not displayed "right".

<    1   2