[R] Coefficient of determination in a regression model with AR(1) residuals

Dear R-users, I used lm() to fit a standard linear regression model to a given data set, which led to a coefficient of determination (R^2) of about 0.96. After checking the residuals I realized that they follow an autoregressive process (AR) of order 1 (and therefore contradicting the i.i

Re: [R] printing to the console in color

If this is RGui on Windows, no. It supports only two colours of text (which you can set in the preferences dialiog box). More generally, you send the appopriate escapes to your console to change colour, and the console/terminal documentation will tell you what those are. On Wed, 23 Apr 2008,

Re: [R] Compiling msm on Fedora Core Linux

The problem seems to be in the header paths: > gcc -m64 -std=gnu99 -I/usr/include/R -I/usr/include/R > -I/usr/local/include Is R.h not in /usr/include/R? My guess is that you installed an R rpm and not the R-devel rpm. But if you installed from rpms, this is a Fedora support issue -- this comm

[R] plotting two functions

i wanna compare functions to be simple , let's say i want x^2 and x^5 in same plot ( it's not the case but if i get it i'll understand for others ) how i do it? x<-seq(-10, 10, l=100) > plot(x^2) and? tks [[alternative HTML version deleted]]

Re: [R] problem with "which"

G'day Melanie, On Wed, 23 Apr 2008 17:46:56 -1000 Melanie Abecassis <[EMAIL PROTECTED]> wrote: > This doesn't seem to happen with integers. > Am I missing something ?? Yes, FAQ 7.31: http://cran.r-project.org/doc/FAQ/R-FAQ.html#Why-doesn_0027t-R-think-these-numbers-are-equal_003f > Is there a

Re: [R] problem with "which"

Melanie Abecassis <[EMAIL PROTECTED]> wrote in news:[EMAIL PROTECTED]: > Hi, > I'm having trouble with the "which" or the "seq" function, I'm not > sure. Here's an example : > > > > lat=seq(1,2,by=0.1) > > lat > [1] 1.0 1.1 1.2 1.3 1.4 1.5 1.6 1.7 1.8 1.9 2.0 > > > which(lat==1) > [1] 1 > > w

[R] problem with "which"

Hi, I'm having trouble with the "which" or the "seq" function, I'm not sure. Here's an example : > lat=seq(1,2,by=0.1) > lat [1] 1.0 1.1 1.2 1.3 1.4 1.5 1.6 1.7 1.8 1.9 2.0 > which(lat==1) [1] 1 > which(lat==1.1) [1] 2 > which(lat==1.2) [1] 3 > which(lat==1.3) [1] 4 > which(lat==1.4) [1]

Re: [R] plotting two functions

On Thu, 24 Apr 2008, Manoel Santos wrote: > i wanna compare functions > to be simple , let's say i want x^2 and x^5 in same plot ( it's not the > case but if i get it i'll understand for others ) > how i do it? See ?curve example( curve ) HTH, Chuck > > x<-seq(-10, 10, l=1

Re: [R] plotting two functions

Hi Manoel, Try this: x<-seq(-10, 10, l=100) plot(x^2,type="l",col=2,ylim=range(x^2,x^5),ylab='f(x)') points(x^5,type="l",col=4) legend('topleft',c(expression(x^2),expression(x^5)),col=c(2,4),lty=1) See ?plot HTH, Jorge On Wed, Apr 23, 2008 at 11:23 PM, Manoel Santos <[EMAIL PROTECTED]> wrot

Re: [R] plotting two functions

by the way, i want both in the same figure, not side by side or others.. tks [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-pro

Re: [R] Import_from_outputfile

You can issue a read for just the headers and then skip first 3 lines so as to issue a read for just the body. fn <- "myfile.dat" DF <- read.table(fn, skip = 3, na.strings = "--", col.names = read.table(fn, nrows = 1, as.is = TRUE)) On Wed, Apr 23, 2008 at 9:29 AM, Ole Roessler <[EMAIL PROTE

Re: [R] Import_from_outputfile

One way is to first read in the data with readLines. This will create a character vector of the input lines. You can then remove the lines you want and then use textConnect to read in the modified data: x <- readLines(file) x <- x[-c(2,3)] # delete lines two and three input <- read.table(textCo

Re: [R] Replecing empty values with a letter (continuation)

Judith I think you'll find the issue is that read.csv creates a data.frame and your subsetting syntax dat[is.na(dat)] is assuming it to be a matrix. If this is the case, you have two options: Coerce dat into a matrix with as.matrix(dat); Use apply() to replace the missing data for each column in

Re: [R] Replecing empty values with a letter (continuation)

Judith Flores <[EMAIL PROTECTED]> wrote in news:[EMAIL PROTECTED]: > Sorry, I hit the send botton by accident. Here is my > code so far: > > dat<-read.csv('myfile.csv', na.strings="") > dat[is.na(dat]<-'N' > >But it doesn't replace the for the letter 'N'. > > What happens when you just t

Re: [R] Can I get rid of this for loop using apply?

Sure, I am creating a partial dependence plot (reference Friedman's stochastic gradient paper from, I want to say, 2001). The idea is to find the relationship between one of the predictors, say x1, and y by creating the following plot: take a random sample of actual data points, hold other predict

Re: [R] Can I get rid of this for loop using apply?

Thanks for the help. That explains why my time testing showed no difference. Is there any way to speed up the program? It is unbearably slow if I increase the number of loops. Mike On Wed, Apr 23, 2008 at 6:23 PM, hadley wickham <[EMAIL PROTECTED]> wrote: > On Wed, Apr 23, 2008 at 4:23 PM, M

Re: [R] Can I get rid of this for loop using apply?

On Wed, Apr 23, 2008 at 7:31 PM, Mike Dugas <[EMAIL PROTECTED]> wrote: > Thanks for the help. That explains why my time testing showed no > difference. Is there any way to speed up the program? It is unbearably > slow if I increase the number of loops. Could you explain exactly what you're tryi

Re: [R] select rows from data based on a vector of char strings

I'm sorry about that. Second line should be attach(yourdata) Cheers, Jorge On Wed, Apr 23, 2008 at 7:48 PM, Jorge Ivan Velez <[EMAIL PROTECTED]> wrote: > > Try this, > > x="label freq1 freq2 > news 54 35 > fun 37 21 > milk19 7 > food 3 3" > > yourdata

Re: [R] select rows from data based on a vector of char strings

u... One of the S language's great strengths is the variety and flexibility of it's data structures. It pays to familiarize yourself with them. So if x is already a dataframe ("dataset" is a rather unhelpful description), then row.names(x) <- x$label x[c("fun","food"),] will extract the row

[R] mysterious (to me) automatic loading of packages when R starts up

Hello, I have recently noticed that in one of my work areas, a number of packages are automatically loaded without my explicitly requesting them to be loaded (see below). This only happens in one particular workspace (located in a folder). I suspect (without any real evidence) that this may be rela

Re: [R] select rows from data based on a vector of char strings

Try this, x="label freq1 freq2 news 54 35 fun 37 21 milk19 7 food 3 3" yourdata=read.table(textConnection(x),header=TRUE) attached(yourdata) flist<-c("fun","food") yourdata[label %in% flist,] I hope this helps, Jorge On Wed, Apr 23, 2008 at 3:13 AM,

Re: [R] select rows from data based on a vector of char strings

Hi, You are right the == doesn't work, but there's a workaround using regular expressions: flist<-"fun|food" grep(flist, data$label) will give you the vector [2 4] which are the numbers of the rows of interest! Dirkheld wrote: > > Hi, > > I have loaded a dataset in R : > data = > > label

[R] filled contour plots (link)

I have moved the figure to the link below. thank you Be a better friend, newshound, and __ R-help@r-project.org mailing list https://stat.ethz.ch/mailma

[R] (no subject)

Hello, I have been trying to apply some of the different ways suggested in the past to replace empty values in a vector for a letter, without success. Actually, I have a data frame that might contain empty values, I need to replace all those empty values for the letter 'N'. This was my l

[R] Replecing empty values with a letter (continuation)

Sorry, I hit the send botton by accident. Here is my code so far: dat<-read.csv('myfile.csv', na.strings="") dat[is.na(dat]<-'N' But it doesn't replace the for the letter 'N'. I am using R v 2/7/0, running it under Windows XP. Thank you, Judith __

[R] help on coxph.wtest

Hi, i need to use pspline. In this pspline function coxph.wtest was used. When I try to make some change to this function by pulling out the pspline function, it turns out R gave me an error msg, saying coxph.wtest cannot be found. Even if i dont change anything in pspline and just rename it and ru

[R] printing to the console in color

Hi, Is there a function that would allow me to print things in different colors to the CONSOLE, say one line in blue and another line in green? Right now, the print() function only prints in navy blue... Thanks! -- View this message in context: http://www.nabble.com/printing-to-the-console-in

Re: [R] combining two (or more) tables by creating another dimension

if(!require(abind)) { install.packages("abind") library(abind) } t3 <- as.table(abind(t1, t2, rev.along = 0)) ### why must they be tables? Bill Venables CSIRO Laboratories PO Box 120, Cleveland, 4163 AUSTRALIA Office Phone (email preferred): +61 7 3826 7251 Fax (if absolutely

Re: [R] Can I get rid of this for loop using apply?

On Wed, Apr 23, 2008 at 4:23 PM, Mike Dugas <[EMAIL PROTECTED]> wrote: > The answer to my post is yes (which I just figured out). > Switching from for to apply isn't going to speed up your code. If you carefully read the source code of apply, you'll see the guts of the work is done by: for (i i

Re: [R] help on coxph.wtest

> Hi > i need to use pspline. In this pspline function coxph.wtest was used. When > I try to make some change to this function by pulling out the pspline > function, it turns out R gave me an error msg, saying coxph.wtest cannot be > found. Even if i dont change anything in pspline and just rename

[R] filled contour plots

The see the figure I refer to on my email, please use this link: http://us.f13.yahoofs.com/bc/480fa972_411e/bc/My+Documents/filled_contour_plot.pdf?bfy26DIBonZ1vybc thanks >hello everyone, >I was wondering if anybody can help me solve 2 > problems related to the function filled.con

[R] Can I get rid of this for loop using apply?

The answer to my post is yes (which I just figured out). Solution: #super small version of R code for pd plot using apply a <- rbind(c(0:)*(max(m$x1)-min(m$x1))/ + min(m$x1),c(0:)*0-9) b <- matrix(rep(c(0:)*(max(m$x1)-min(m$x1))/ + min(m$x1), nrow(m)), nrow(m), 1112, b

[R] 64 bit Linux/Ubuntu memory limit problem?

Dear list, I've recently installed R on a 64 bit machine with 8 GB of RAM. I set this computer up as a dual-boot system, with windows XP 64 and Ubuntu 7.10. I downloaded the Linux 64 bit version of R and installed it. I'm trying to run rather large Random forest models and was running into

[R] Compiling msm on Fedora Core Linux

Greetings all, I'm trying to install the msm package and it is failing on compilation. The problem seems to be the analyticp component? Any advice on how to get it to work? error message is below. I'm running R version 2.6.2 (2008-02-08) x86_64-redhat-linux-gnu on a dell precision 690 with Fe

[R] filled contour plots

sorry. the figure did not go through. sending again as PDF. hope it will work this time. Thanks Be a better friend, newshound, and __ R-help@r-project.

[R] filled contour plots

hello everyone, I was wondering if anybody can help me solve 2 problems related to the function filled.contour. I am entering the following R command: filled.contour(xx,yy,P1, nlevels=20,color=cm.colors, plot.axes={ contour(xx,yy,P1,add=T,col="grey", nlevels=20, draw

[R] BB - a new package for solving nonlinear system of equations and for optimization with simple constraints

Hi, We (Paul Gilbert and I) have just released a new R package on CRAN called "BB" (stands for Barzilai-Borwein) that provides functionality for solving large-scale (and small-scale) nonlinear system of equations. Until now, R didn't have any functionality for solving nonlinear systems. We ho

Re: [R] Multidimensional contingency tables

It seems to me that a combination of ftable and xtabs works fine: > prob1<- data.frame(victim=c(rep('white',4),rep('black',4)), + perp=c(rep('white',2),rep('black',2),rep('white',2),rep('black',2)), + death=rep(c('yes','no'),4), count=c(19,132,11,52,0,9,6,97)) > prob1 victim perp death count 1

[R] help on coxph.wtest

Hi, i need to use pspline. In this pspline function coxph.wtest was used. When I try to make some change to this function by pulling out the pspline function, it turns out R gave me an error msg, saying coxph.wtest cannot be found. Even if i dont change anything in pspline and just rename it and ru

Re: [R] ccf and covariance

On Wed, 23 Apr 2008, Bob Farmer wrote: > Thanks to Prof. Ripley and Phil Spector for pointing out that the > autocorrelation functions must use a "nontraditional" definition of > the covariance, involving a denominator of n (instead of n-1) in order > to satisfy an assumption of second-order stati

[R] combining two (or more) tables by creating another dimension

Dear R community, I wish to combine two tables in one by adding an additional dimension: e.g.: > t1<-as.table(matrix(rnorm(40),nrow=4,ncol=10));rownames(t1)<-c("rowone","rowtwo","rowthree","rowfour") > t1 A B C D E F G

Re: [R] Feed list of vectors to vioplot()?

Off-list it was pointed out to me that the trivial solution to this would look like this: myfunction2 <- function(x, ...){ mylist <- list(x, ...) # plenty of "lapply" stuff do.call(vioplot,mylist) } Thanks for everybodies patience, Joh Johannes Graumann wrote: > > > Thanks for the very

Re: [R] ccf and covariance

Thanks to Prof. Ripley and Phil Spector for pointing out that the autocorrelation functions must use a "nontraditional" definition of the covariance, involving a denominator of n (instead of n-1) in order to satisfy an assumption of second-order stationarity in the (unbiased) covariance estimators

[R] Minimise a parameter of a given function f, with f > 0

Hi, I need to find the minimum value of the parameter, s, such that the function f(s,t) > 0 (where -Inf < t < Inf) I've looked into optim, constrOptim and others but they don't seem to do this. Des anyone have some suggestions? Thanks in advance. [[alternative HTML version deleted]] __

[R] Can I get rid of this for loop using apply?

Hey all, The code below creates a partial dependence plot for the variable x1 in the linear model y ~ x1 + x1^2 + x2. I have noticed that the for loop in the code takes a long time to run if the size of the data is increased. Is there a way to change the for loop into an apply statement? The tr

Re: [R] Feed list of vectors to vioplot()?

Thanks for the very appropriate scolding. Here's my example (based on "?vioplot"): mu<-2 si<-0.6 bimodal<-c(rnorm(1000,-mu,si),rnorm(1000,mu,si)) uniform<-runif(2000,-4,4) normal<-rnorm(2000,0,3) # Working just fine myfunction1 <- function(x, ...){vioplot(x,...)} myfunction1(bimodal,uniform,norm

Re: [R] ccf and covariance

On Wed, 23 Apr 2008, Bob Farmer wrote: > Hi. > It's my understanding that a cross-correlation function of vectors x > and y at lag zero is equivalent to their correlation (or covariance, > depending on how the ccf is defined). The ratio of your values is > MASS::fractions(282568.5/259021) [1] 12

[R] significant variables in GPLS ?

Hello, I am using the gpls package for modelling vegetation classes. My problem is that I now want to know which input variables are significant for the modelling of the classes to recalculate the equation again with just the selected variables. I think I can analyse the significance of the vari

[R] ccf and covariance

Hi. It's my understanding that a cross-correlation function of vectors x and y at lag zero is equivalent to their correlation (or covariance, depending on how the ccf is defined). If this is true, could somebody please explain why I get an inconsistent result between cov() and ccf(type = "covarianc

Re: [R] help needed: Plotting step by step.

Maybe like this? > proc.1 <- rexp(200,1) > proc.2 <- rexp(200,2) > plot( ecdf( proc.1 ), xlim=range( proc.1, proc.2 ) ) > plot( ecdf( proc.2 ), add=T, col.points='red' ) > See ?ecdf ?plot.stepfun HTH, Chuck On Wed, 23 Apr 2008, Greg Snow wrote: > Can you show us the code yo

Re: [R] Feed list of vectors to vioplot()?

PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. It is hard to provide a solution if we do not understand the problem to be solved. Sample data would be helpful along with an understanding of what you w

Re: [R] pdf() and histogram() in function call

Read FAQ 7.22 -- Gregory (Greg) L. Snow Ph.D. Statistical Data Center Intermountain Healthcare [EMAIL PROTECTED] (801) 408-8111 > -Original Message- > From: [EMAIL PROTECTED] > [mailto:[EMAIL PROTECTED] On Behalf Of qian z > Sent: Wednesday, April 23, 2008 10:49 AM > To: r-help@r-pr

[R] Course*** R/S-Plus Advanced Programming ****by XLSolutions Corp / May 2008 in San Francisco

Announcing R/Splus Advanced Programming course: May 29-30, 2008 in San Francisco. *** half of the seats are already taken *** Please email for earlybird rates: Payments due after the class. R/Splus Advanced Programming Course Outline: Day 1 - Overview of R/S fundamentals: Syntax a

Re: [R] help needed: Plotting step by step.

Can you show us the code you used for the 5 to 10 points? (either generate some random data, or use a sample dataset). -- Gregory (Greg) L. Snow Ph.D. Statistical Data Center Intermountain Healthcare [EMAIL PROTECTED] (801) 408-8111 > -Original Message- > From: [EMAIL PROTECTED] > [

Re: [R] g(x,y) = f(x,y) - e(x)- e(y)?

It is not clear exactly what you want to do by subtracting a vector from a matrix, but look at the "sweep" function, it may do what you want. -- Gregory (Greg) L. Snow Ph.D. Statistical Data Center Intermountain Healthcare [EMAIL PROTECTED] (801) 408-8111 > -Original Message- > From

[R] pdf() and histogram() in function call

Here is a function I wrote. It runs no problem, but generate empty pdf files. I can't find what is the problem. create.pdf<- function(x, dir) { dir.create(dir, showWarnings = FALSE) plist<- c("a", "b" , "c", "d") for(j in plist) { filedir<- paste(dir,

[R] nested time series data with measurement error

Dear R helpers, I am trying to fit a model with the main objective of assessing differences, rather than predicting. The treatment was applied to half of the subjects after 9 months of measurement. Then 9 more months of data were collected. The variable "month" has measurement error due to the com

Re: [R] Density estimation

Look at the logspline package. It has a different way of estimating densities that allows for limits to be specified (i.e. probability is 0 beyond the point(s) you specify). Hope this helps, -- Gregory (Greg) L. Snow Ph.D. Statistical Data Center Intermountain Healthcare [EMAIL PROTECTED] (801)

Re: [R] Overall p-value from a factor in a coxph fit

Kåre Edvardsen <[EMAIL PROTECTED]> wrote in news:[EMAIL PROTECTED]: > Prof. Paul, Prof. Frank. > > Thank you very much for helping me out. The Design package did the > trick. > > Here is how the anova table looks like without using the Design > package: > >> anova(Fit1) > Analysis of Devianc

Re: [R] Equivalent of intervals() in lmer

Dear Friends, (1) There may be a solution for those (e.g., experimental psychologists) who are *not at all* interested in generalizing the absolute level of the response variable (say, reaction time, rt) to other subjects, but *only* to generalize the effect of the manipulated variables wi

Re: [R] Comparing kendall's tau values?

Ashton, Gail si.edu> writes: > > I have 3 variables relating to the successful introductions of species > to 95 different areas: introduction frequency; number of successes pre > 1906; number of successes post 1906 > > The data are not normal, nor homo-skedatic, so I am using non-parametric > s

[R] help needed: Plotting step by step.

Hello, I have generated 2 Poisson processes and want to plot them on a single graph in a step by step manner in order to be able to compare them. I tried plot and biplot but it does not help, I could connect two points by hand for point graph if they were 5 or 10 I have more than 200 such point to

Re: [R] Modeling presence only data in R

Actually, the ENFA (function enfa, in the package adehabitat) is not really suitable for modelling environment suitability maps (though it is still useful to discover the patterns of the ecological niche with presence-only data). However, the package adehabitat also contains other functions all

Re: [R] Modeling presence only data in R

Dear Alejandro, Thank you very much for your advice. In fact I use GARP and Openmodeller, and they - up to I know - also generate pseudo-absence when modelling. I am searching some R solution to compare the results with OpenModeller outputs. Could you send-me a PDF (or the complete reference) of L

Re: [R] g(x,y) = f(x,y) - e(x)- e(y)?

> Thanks Phipp very much for your help. I had meant, given that I'd > computed the matrix f[x,y] and the vector e[x], how to take the > difference. What is confusing is how to subtract a vector from a > matrix. I don't want the recycling rule. That function sounds like its describing how to sub

Re: [R] Modeling presence only data in R

Ecological Niche Factor Analysis (ENFA) is implemented in the adehabitat package available on CRAN. Cheers, Alejandro González wrote: > Dear Milton, you have 2 packages for modelling species distribution in > R: grasper and Biomod (by Wilfried Thuiller), but they are all > presence-absence m

Re: [R] Feed list of vectors to vioplot()?

Johannes Graumann wrote: > Hi, > > I have a list of vectors and am trying to coerce them into something that > vioplot will take as groups of data to be plotted independently. Can > someone nudge me into the right direction? > > Thanks, Joh > > __ > R

Re: [R] Modeling presence only data in R

Dear Milton, you have 2 packages for modelling species distribution in R: grasper and Biomod (by Wilfried Thuiller), but they are all presence-absence models, so you must generate pseudoabsences for each species, following recommendations in papers (see Lobo 2007). On the other side, there are

[R] Luis Miguel Delgado Gomez/BBK está ausente d e la oficina.

Estaré ausente de la oficina desde el 22/04/2008 y no volveré hasta el 11/05/2008. Responderé a su mensaje cuando regrese. [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help

Re: [R] Writing list object to a file

"Arun Kumar Saha" <[EMAIL PROTECTED]> wrote: > Hi all, > > I am wondering how to write a 'list' object to a file. I already gone > through some threads like > http://mail.python.org/pipermail/python-list/2001-April/080639.html, however > could not trace out any reliable solution. I tried followin

Re: [R] g(x,y) = f(x,y) - e(x)- e(y)?

Le mer. 23 avr. à 06:59, William Simpson a écrit : > Thanks Phipp very much for your help. I had meant, given that I'd > computed the matrix f[x,y] and the vector e[x], how to take the > difference. What is confusing is how to subtract a vector from a > matrix. I don't want the recycling rule. We

Re: [R] lmer model building--include random effects?

On Apr 23, 2008, at 8:56 AM, Douglas Bates wrote: > On 4/22/08, Ista Zahn <[EMAIL PROTECTED]> wrote: >> Hello, >> This is a follow up question to my previous one >> http://tolstoy.newcastle.edu.au/R/e4/help/08/02/3600.html > >> I am attempting to model relationship satisfaction (MAT) scores >> (

Re: [R] ROracle error at step 1

Hello and thanks for the quick reply. I got it here, http://stat.bell-labs.com/RS-DBI/download/index.html As the following site http://cran.r-project.org/web/packages/ROracle/NEWS mentions that the changes to the newest version are trivial. Although I take your point that pack

[R] Import_from_outputfile

Hello, I want to do some analyses with data originating from a hydrological model that are saved in multiple output files. I want to import the output-files in R without a treatment like preproceesing every output file. The output file is more or less a matrix, but with ad

Re: [R] ROracle error at step 1

Where did you get your ROracle package from? > Built: R 2.2.1; i386-pc-mingw32; 2006-02-13 18:05:59; windows indicates it was too old. All methods-using packages need to have been installed recently, and definitely since 2.4.0. (It is unsafe to use one that has not been installed for the same

Re: [R] lmer model building--include random effects?

On 4/22/08, Ista Zahn <[EMAIL PROTECTED]> wrote: > Hello, > This is a follow up question to my previous one > http://tolstoy.newcastle.edu.au/R/e4/help/08/02/3600.html > I am attempting to model relationship satisfaction (MAT) scores > (measurements at 5 time points), using participant (spouse

Re: [R] Writing list object to a file

?save On Wed, Apr 23, 2008 at 2:51 AM, Arun Kumar Saha <[EMAIL PROTECTED]> wrote: > Hi all, > > I am wondering how to write a 'list' object to a file. I already gone > through some threads like > http://mail.python.org/pipermail/python-list/2001-April/080639.html, however > could not trace out any

Re: [R] how to read in multiple files with unequal number of columns

Is this what you want? I am assuming that you will read the dataframes into a list and then process them like below: > # put dataframe in a list -- would have read them in via a list > x <- list(d, my.fake.data) > # determine maximum number of columns and then pad out the short one > # also use t

[R] ROracle error at step 1

Hi I Can't connect to the Oracle database, any tips? Has anybody actually got ROracle up and running on windows? >> unable to find an inherited method for function "dbConnect", for signature "OraDriver" I can happily connect to the same database through RODBC. Oracle client, version 9.2

Re: [R] Ubuntu vs. Windows

On 4/22/08, Prof Brian Ripley <[EMAIL PROTECTED]> wrote: > On Tue, 22 Apr 2008, Peter Dalgaard wrote: > > > Doran, Harold wrote: > >> Dear List: > >> > >> I am very much a unix neophyte, but recently had a Ubuntu box installed > >> in my office. I commonly use Windows XP with 3 GB RAM on my ma

Re: [R] Time arithmetic

library(zoo) ?yearmon library(chron) ?month.day.year ?strptime See R News 4/1. On Wed, Apr 23, 2008 at 6:13 AM, Worik R <[EMAIL PROTECTED]> wrote: > I am a bit worried I am reinventing the wheel. Isn't there a calendar > system in R? > > I have written a function to add months to a date and re

Re: [R] Problems with Windows RGUI R-2.7.0 when P ATH includes 'æ'

Your email was not properly encoded by the time it reached me (see the ? in crucial places), so I may be missing something. I am guessing you are in a Danish locale with codepage CP1252, and fortunately the subject line appears to have worked. I think I have a reproduction test, in which case

[R] Density estimation

Hi, I am analysing a dataset containing genetic distances within and between species. I want to show a overlap of the distributions of the intra- and interspecific values; on a second graph I use a cut-off value to determine these boundaries. As the dataset contains >30 000 values, I would like

Re: [R] g(x,y) = f(x,y) - e(x)- e(y)?

Thanks Phipp very much for your help. I had meant, given that I'd computed the matrix f[x,y] and the vector e[x], how to take the difference. What is confusing is how to subtract a vector from a matrix. I don't want the recycling rule. Cheers Bill On Tue, Apr 22, 2008 at 9:53 AM, Philipp Pagel <[

[R] [R-pkgs] new package multipol

Hello List please find a new package, multipol, recently uploaded to CRAN. This package generalizes the polynom package (which handles univariate polynomials) to the multivariate case. A short article discussing the package will appear in the next issue of Rnews, Insha'Allah enjoy -- Rob

[R] R appears not to reallocate memory on linux

Hi, I have a question concerning the memory management of R 2.6.1 on a 32-bit linux with 3.2 GB RAM. Especially, I found the current memory size reported from gc() to be different from the amount reported by e.g. top. The 'max used' number of gc() seems more to match the actually used memory. s

[R] Time arithmetic

I am a bit worried I am reinventing the wheel. Isn't there a calendar system in R? I have written a function to add months to a date and return the number of days resulting. I am newish to R so I am hoping there is a package that can do this sort of date arithmetic for me... Worik DaysInMonth

Re: [R] optimization setup

The error comes from the way you specify the parameters: At least not as elegant, but it works: function4 <- function(x){ theta1 <- x[1] theta2 <- x[2] theta3 <- x[3] function3(theta1,theta2,theta3) } fit<-optim(par=c(1, 1.2, .2), fn=function4) Bart threshol

Re: [R] how to read in multiple files with unequal number of columns

Thank you John. It was useful to know about this package. I tried merge_all and I got this error: Error in .subset2(x, i, exact = exact) : subscript out of bounds It could be due to the way my data is and I will try the other solutions suggested by the other kind souls on this list. Best wish

[R] Problems with Windows RGUI R-2.7.0 when PA TH includes 'æ'

I have a subdirectory named 'Fælles Filer' in my PATH. When I start the RGUI, I receive the following error/warning messages: Error in Sys.setenv(PATH = PATH) : invalid input in wtransChar and Warning message: package "methods" in options("defaultPackages") was not found The problems disappear

Re: [R] find directories

Prof Brian Ripley a écrit : > On Wed, 23 Apr 2008, mel wrote: > >> Hello, >> Function dir() is very useful for finding files. >> However, is there a analog function for finding directories >> based eg on a pattern ? > > No. > >> Apologies if I've missed something obvious. > > list.dirs <- funct

Re: [R] find directories

On Wed, 23 Apr 2008, mel wrote: > Hello, > Function dir() is very useful for finding files. > However, is there a analog function for finding directories > based eg on a pattern ? No. > Apologies if I've missed something obvious. list.dirs <- function(...) { x <- dir(...) x[file_test(

Re: [R] Regression inclusion of variable, effect on coefficients

At 20:06 21/04/2008, Thiemo Fetzer wrote: >Hello :) > >I am happy to hear that I am not necessarily asking stupid questions. > >The thing is, that I have data on x1 and x4 for the whole sample. However, >theoretically, it is clear that the informational content of x1 is not as >high as of x4. x4 pr

Re: [R] logit GLM without intercept

At 16:56 21/04/2008, Robert Junker wrote: >Dear Statisticians, > > > >I would like to analyse my data with a GLM with binomial error distribution >and logit link function. The point is that I want a model fitted without >intercept, i.e. the fitted curve should start at y=0.5 for x=0. > >I tried it

Re: [R] Documentation General Comments

Good morning, Firstly I'd like to say that I'm a huge fan of R and I think it's great system. Part of the problem in searching for information is knowing what buzzwords / keywords to use. I was recently caught out like this as I didn't see my problem as a cumulative sum (keyword=cumsum) only as

Re: [R] select rows from data based on a vector of char strings

On 4/23/2008 3:13 AM, Dirkheld wrote: > Hi, > > I have loaded a dataset in R : > data = > > label freq1 freq2 > news 54 35 > fun 37 21 > milk19 7 > food 3 3 > etc > > And I have a vector > flist<-c("fun","food") > > Now I want to use the vector 'fli

[R] select rows from data based on a vector of char strings

Hi, I have loaded a dataset in R : data = label freq1 freq2 news 54 35 fun 37 21 milk19 7 food 3 3 etc And I have a vector flist<-c("fun","food") Now I want to use the vector 'flist' for selecting these values from 'data' so that I get the followin

Re: [R] Writing list object to a file

Previous question that I asked was originated from this problem : suppose I have following time series : library(zoo) date1 = seq(as.Date("01/01/01", format = "%m/%d/%y"), as.Date("12/31/08", format = "%m/%d/%y"), by = 1) len1 = length(date1); data1 = zoo(matrix(rnorm(len1, mean=0, sd=0.5), nrow

Re: [R] problem with data for metaMDS analysis please help

On Tue, 2008-04-22 at 20:19 -0400, stephen sefick wrote: > am at my wit's end. I am not sure what is wrong with this data matrix. It > is sparse because it is a matrix of species, but I have looked at the row > totals and column totals and they are positive. > rmetaMDS(x.d) > Error in if (autotra

[R] find directories

Hello, Function dir() is very useful for finding files. However, is there a analog function for finding directories based eg on a pattern ? Apologies if I've missed something obvious. Thanks Vincent __ R-help@r-project.org mailing list https://stat.ethz.

  1   2   >