[R] Problems with F77_CALL and Visual C++

I am working on a DLL and I use Visual C++. In the Visual c++ code I use F77_CALL(dcopy)(.); Visual C++ gives the error 'daxpy_' identifier not found. Do you know hot to solve this problem? Thanks in adavnce [[alternative HTML version deleted]] _

[R] Create *png or *.jpeg plot under Linux

I have an example where i plot an image and load it dynamically in an image analysis tool to create a film (dynamic graph). This works well under Windows but i can't plot to *.png or *.jpeg files under Linux (Ubuntu) with the same command. I'm using: png(file,width=200,height=200) or jpeg(fil

Re: [R] Error: X11 cannot allocate additional graphics colours.

On Mon, 15 Oct 2007, michael watson (IAH-C) wrote: Thank you both. I use R on linux both remotely (over CGI) and I log in using Exceed from windows. The problem occurs over both. The latter may explain the font problem: it is a crucial fact. The X server is running on the Windows box, not

Re: [R] Image analysis and image questions

Thank for the answer but i'm searching for built in (matrix operations) functions (filters) or functions for pattern recognition to identify special regions or shapes (in my case species, regions, boundary's). I can then transfer the results back and mark the image in ImageJ. In addition i found o

Re: [R] Wildcards

On Mon, 15 Oct 2007, jim holtman wrote: > for (i in list.files(pattern=".*\\.R$")) source(i) Using Sys.glob() is perhaps more natural for most users: for(file in Sys.glob("*.R")) source(file) It is referenced from ?list.files these days. > On 10/15/07, subura <[EMAIL PROTECTED]> wrote: >> >>

Re: [R] survreg's algorithm

Did you look at the C source code? There are 4 different variants (survregN.c, where N <- 2:5) , depending on whether the distribution is built-in or not, and penalized likelihood is being used or not. They all look like NR to me, but I confess I haven't read the code in extreme detail. It is well

Re: [R] partitioning data [SEC=UNCLASSIFIED]

Hi Stephen, Check the help for predict.glm(). The argument for passing new data is actually 'newdata', as in: > pred = predict(glm.model, newdata=form[150001:20,-1], > type="response") Cheers Joe -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of

Re: [R] Need some help

> -Original Message- > From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf > Of azzza > Sent: Monday, October 15, 2007 6:06 PM > To: r-help@r-project.org > Subject: Re: [R] Need some help > > > > Thanks Jholtman. > However, the plot didnt come out the way I envisone dit to be. On

[R] Bootstrapping Contrasts for Repeated Measures ANOVA

I have executed a Repeated Measures ANOVA with one DV (latency) and one within subject factor (acoustic condtion: 3 levels) by bootstrapping my sampling distribution of F from the empirical sample distribution. I chose to resample because the sample distribution deviates from normality a lot. The

Re: [R] Linear regression and slope from 1

Dan-157 wrote: > > > Dear R-Users, > > I am new to R, so please excuse the ignorance. > > I have data: > x (2.14, 2.41, 1.09, 0.17, 8.18) > y (3.81, 5.13, 0.63, 0.75, 6.35) > > I would like to use simple linear regression and test 2 things: > 1) slope of line of best fit is statistically di

Re: [R] survreg's algorithm

Gad Abraham wrote: > Hi, > > I'm using survreg() from the survival package for parametric survival > regression (modelling inter-arrival times of patients to a waiting list > as exponentially distributed, with various regressors such as queue size > and season). > > Does anyone know which algo

Re: [R] Need help ploting time series(2)

Here are a couple of ways that you can do it: x <- expand.grid(YEAR=c(2003,2004), MONTH=1:12, DAY=1, STATE=LETTERS[1:4]) x$SALES <- runif(nrow(x), 10, 100) # add date for plotting x$date <- ISOdate(x$YEAR, x$MONTH, x$DAY) # sort into date order for plotting x <- x[order(x$date),] # you can use lat

[R] survreg's algorithm

Hi, I'm using survreg() from the survival package for parametric survival regression (modelling inter-arrival times of patients to a waiting list as exponentially distributed, with various regressors such as queue size and season). Does anyone know which algorithm survreg() uses for this? Tha

Re: [R] Need some help

Thanks Jholtman. However, the plot didnt come out the way I envisone dit to be. On the Y axis, i should have sign changes in 1000 tosses, the range being from negative to postitive, and a straight horizontal line across y=0. The X-axis should have the toss number, range 0-1000 I'm glad u mentio

[R] Need help ploting time series(2)

hi: Yesterday I post a message about hoy to plot a time series, but someone told me to post more information about the file so here it is: the file was read using read.table and the name is list. When I use str(list) it tells the following variables: YEAR int: 2003,2003,20032004 MONTH int:1,1,

[R] Linear regression and slope from 1

Dear R-Users, I am new to R, so please excuse the ignorance. I have data: x (2.14, 2.41, 1.09, 0.17, 8.18) y (3.81, 5.13, 0.63, 0.75, 6.35) I would like to use simple linear regression and test 2 things: 1) slope of line of best fit is statistically different from 1 2) y-intercept is statically

[R] Distance matrix in SpDep-package

Hello everybody, I would like to use the SpDep-package (especially the Local Moran index analysis and the Getis-Ord statistics) in R for analysing my data. However, I don't have x-y coordinates, but my data is in a distance matrix format. Is it possible to use the SpDep package with predefined

Re: [R] Error: X11 cannot allocate additional graphics colours.

On Mon, 15 Oct 2007, michael watson (IAH-C) wrote: > Thanks for the response... > > My confusion about plot stems from the fact I am plotting 82 points with 82 > colours, so surely all colours get plotted? > [snip] No. That was the point of this bit: > > hist( mat, breaks=n ) > Or if tha

Re: [R] Make playwith a default graphic device

On 10/15/07, Felix Andrews <[EMAIL PROTECTED]> wrote: > My previous suggestion was inconsistent with the Trellis/Lattice idea > of creating a trellis object without necessarily creating a plot. And > it also interfered with attempts to plot to a file device. So here is > a better solution, based on

Re: [R] Error: X11 cannot allocate additional graphics colours.

Thank you both. I use R on linux both remotely (over CGI) and I log in using Exceed from windows. The problem occurs over both. I don't run Xvfb (though I can do as I have done in the past). I have root priveleges so I will look and see what the settings are, though I don't think my pet sy

Re: [R] Error: X11 cannot allocate additional graphics colours.

(This could be fixed, it has not happened to me in a long time, but I will mention it mainly because it is not something you are likely to think of.) It used to be that the X colours might be defined by the first application that needed them, so if the systems administrator happened to start u

Re: [R] some question about partial prediction in survival

For numerical accuracy, the coxph routine centers each covariate before doing the computation. All of the downstream results (predict, survfit, etc) use this centered data. Terry Therneau __ R-help@r-project.org mailing list https://stat.e

Re: [R] coef se in lme

?vcov -Original Message- From: [EMAIL PROTECTED] on behalf of Irene Mantzouni Sent: Mon 10/15/2007 3:20 PM To: [EMAIL PROTECTED] Subject: [R] coef se in lme Hi all! How is it possible to estimate standard errors for coef obtained from lme? Is there sth like se.coef() for lmer or what

Re: [R] coef se in lme

Let's say that I can estimate se for ranef by group and the se for fixef (modifying the se.fixef and se.ranef functions of "arm" package for lmer). Since, I need the se for coef=ranef+fixef can I estimate it based on the SEs of ranef and fixef? Thank you! Irene __

Re: [R] Error: X11 cannot allocate additional graphics colours.

This is not a X server version issue. There are Linux system (not R) display settings that will dictate the number of simultaneous colors that can be displayed. This will be dependent upon the display resolution defined and the amount of video RAM on the graphics card. The higher the display resol

[R] help with simple goodness of fit test

Hello, I've read the other posts with regard to "chisq.test" and "goodness of fit" and am still missing something. 1. I create a simple vector of randomly generated lognormal values with mean=0 and sd=1; >d1 <- rlnorm(100,meanlog=0,sdlog=1); 2. I also create a vector of probabilities that are expe

Re: [R] Get data from matrix

SubMatrix = Matrix[10,] Cheers, Jorge - Original Message - From: pintinho <[EMAIL PROTECTED]> Date: Monday, October 15, 2007 2:17 pm Subject: [R] Get data from matrix > > Hi, > > I have a matrix that has a variable number of columns. I do not > know, a > priori, the number of c

Re: [R] Error: X11 cannot allocate additional graphics colours.

Thanks for the response... My confusion about plot stems from the fact I am plotting 82 points with 82 colours, so surely all colours get plotted? As for updating X, I recently installed the latest version of XFree86 for my version of linux, RHEL 4. As for Brian's e-mail you quoted, I do try a

Re: [R] Strange results converting string to number

On 16/10/2007, at 8:30 AM, pintinho wrote: > > Hi, > > I am getting a strange result while converting a string vector into > numeric > vector: > >> Datas[1] > [1] 37315 > >> as.numeric(Datas[1]) > [1] 2 > > Can anyone help me?? It would seem that ``Datas'' is a ***factor*** and NOT a ``string

[R] Strange results converting string to number

Hi, I am getting a strange result while converting a string vector into numeric vector: > Datas[1] [1] 37315 > as.numeric(Datas[1]) [1] 2 Can anyone help me?? -- View this message in context: http://www.nabble.com/Strange-results-converting-string-to-number-tf4629811.html#a13220089 Sent from

Re: [R] Rmpi building

Dirk Eddelbuettel wrote: > Sean, > > On 15 October 2007 at 12:03, Sean Davis wrote: > | I am trying to build Rmpi on Suse 10.2 linux. While I would like to use > | the RPM version of liblam-7.1.2, I have not been able to do so. It > | seems that Rmpi makes some pretty strong assumptions when try

[R] partitioning data

I am trying to train on part of my data and test on another part: > glm.model = glm(as.factor(h_finished) ~ . , family=binomial, data=form[1:15,]) > pred = predict(glm.model, data=form[150001:20,-1], type="response") > t = table(pred, form[150001:20,1]) Error in table(pred, form[15000

[R] coef se in lme

Hi all! How is it possible to estimate standard errors for coef obtained from lme? Is there sth like se.coef() for lmer or what is the anaytical solution? Thank you! __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help P

Re: [R] avoiding a loop?

I knew it was simple. Thanks very much. -tom jim holtman <[EMAIL PROTECTED]> wrote: > ?table to count the factors > > > x > [1] "a" "b" "c" "d" "e" > > paste(head(x, -1), tail(x, -1), sep='') > [1] "ab" "bc" "cd" "de" > > > On 10/15/07, Tom Sgouros <[EMAIL PROTECTED]> wrote: > > > > Hi

Re: [R] Wildcards

subura said the following on 10/15/2007 12:04 PM: > Care to explain how i can use a wildcard expression to "source" all files > ending with .R in a subdirectory ? I've tried something like this > 'source(glob2rx("*.R"))' without success. > > Thank you Try R.files <- list.files(my.path, pattern

Re: [R] Wildcards

for (i in list.files(pattern=".*\\.R$")) source(i) On 10/15/07, subura <[EMAIL PROTECTED]> wrote: > > Care to explain how i can use a wildcard expression to "source" all files > ending with .R in a subdirectory ? I've tried something like this > 'source(glob2rx("*.R"))' without success. > > Thank

Re: [R] avoiding a loop?

?table to count the factors > x [1] "a" "b" "c" "d" "e" > paste(head(x, -1), tail(x, -1), sep='') [1] "ab" "bc" "cd" "de" On 10/15/07, Tom Sgouros <[EMAIL PROTECTED]> wrote: > > Hi All: > > I feel like there must be a slick R-native no-loop way to get the counts > for the entries in a factor,

[R] Wildcards

Care to explain how i can use a wildcard expression to "source" all files ending with .R in a subdirectory ? I've tried something like this 'source(glob2rx("*.R"))' without success. Thank you -- View this message in context: http://www.nabble.com/Wildcards-tf4627981.html#a13214214 Sent from the

Re: [R] grouping modalities

> Dear Sirs: > > Is there a way to group multiple responses in variables sets. SPSS > have this feature: this is possible to group a set of variables by > their common categories. I would like to do the same in R. Example: my.df = data.frame(var1=c(1,2,2,1,2,2,1,2), var2=c(2,2,1,1,2,2,1,1),

Re: [R] Get data from matrix

Hi, mat[,10:ncol(mat)] On 15/10/2007, pintinho <[EMAIL PROTECTED]> wrote: > > > Hi, > > I have a matrix that has a variable number of columns. I do not know, a > priori, the number of columns. > > How can I get a sub matrix, for example, from row 10 to the end of the > columns? > > In MatLab I w

Re: [R] Error: X11 cannot allocate additional graphics colours.

You knew this? http://tolstoy.newcastle.edu.au/R/e2/help/06/09/0640.html I cannot replicate your error. I use n <- 1000 on R-2.6.0, and it still works. Only a guess, but maybe your X setup is out of date. Maybe an update would help? As for why axis triggers this, axis uses all the c

Re: [R] Rmpi building

Sean, On 15 October 2007 at 12:03, Sean Davis wrote: | I am trying to build Rmpi on Suse 10.2 linux. While I would like to use | the RPM version of liblam-7.1.2, I have not been able to do so. It | seems that Rmpi makes some pretty strong assumptions when trying to | guess the MPI version. Has

Re: [R] Bad EMF export

I think you mean savePlot: I have never heard of saveCopy. We have seen this before, and the problem was the EMF viewer, not the EMF file. So how is this being viewed? And as ever > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html > and provide commented, minimal, s

[R] avoiding a loop?

Hi All: I feel like there must be a slick R-native no-loop way to get the counts for the entries in a factor, but I'm unable to see how. Right now I'm doing this: hn.ent<-data.frame(rep(0,length(levels(hnf))), row.names=levels(hnf)) colnames(hn.ent)<-c("count") for (lev in levels(hnf)) { hn.e

Re: [R] Get the last 3 chars of a string

I was hoping to avoid using regex except when necessary (u know what they say), but I'm beginning to think that's the way things are done in R. Thanks On 10/15/07, jim holtman <[EMAIL PROTECTED]> wrote: > Looks fine by me. There are lots of other ways of doing it. Does > this look any nicer? >

[R] Get data from matrix

Hi, I have a matrix that has a variable number of columns. I do not know, a priori, the number of columns. How can I get a sub matrix, for example, from row 10 to the end of the columns? In MatLab I would use something like this: SubMatrix = Matrix[10, 1:end] Thanks a lot. -- View this messa

Re: [R] Get the last 3 chars of a string

Looks fine by me. There are lots of other ways of doing it. Does this look any nicer? > x <- c('asdfghk', 'qwerrey') > gsub(".*(...)$", '\\1', x) [1] "ghk" "rey" > On 10/15/07, Sergio Correia <[EMAIL PROTECTED]> wrote: > I want to extract the last 3 letters of a string. > > So far, I've done t

Re: [R] Get the last 3 chars of a string

Here is an alternative: sub(".*(..)$", "\\1", x) and using strapply in gsubfn its even shorter: library(gsubfn) strapply(x, "..$") On 10/15/07, Sergio Correia <[EMAIL PROTECTED]> wrote: > I want to extract the last 3 letters of a string. > > So far, I've done this: > > > symbol = 'XYZ.VX" > > s

[R] Get the last 3 chars of a string

I want to extract the last 3 letters of a string. So far, I've done this: > symbol = 'XYZ.VX" > substr(symbol,nchar(symbol)-2,nchar(symbol)) [1] ".VX" It works, but the code looks UGLY as hell. Am I missing something? Or is this the way it's supposed to be? Thanks, Sergio On 10/15/07, pintinho

Re: [R] ERROR while importing data

Stefan, I found what was wrong. Thank you very much anyway! Stefan Grosse-2 wrote: > > On Monday 15 October 2007 06:43:52 pm pintinho wrote: > pi> I tried lots of methods to import (read.csv, read.table, RODBC, > read.delim) > pi> and the same message appears for all these methods. I thin

[R] residual s.e. by group in mixed model

Hi all! I would like to extract the residual standard error by group in a lme() model... Is there a direct method? Also, a rather statistical question; I need to estimate the standard error of the mean(residuals) in a model..is this the same as the residual s.e.? Thank you!! Irene _

Re: [R] Need some help

I really depends on what you want to plot. You can plot the cumulative count of the sign changes with: x <- sample(c(0,1), 1000, TRUE) plot(cumsum(x), type='l') There are other things you can do. You can get a rough histogram of the length of the run by: > stem(rle(x)$length) The decimal po

Re: [R] ERROR while importing data

On Monday 15 October 2007 06:43:52 pm pintinho wrote: pi> I tried lots of methods to import (read.csv, read.table, RODBC, read.delim) pi> and the same message appears for all these methods. I think it is a bigger pi> problem. pi> pi> Can anyone help me solving this issue? A little bit more in

Re: [R] Bad EMF export

On Mon, 2007-10-15 at 18:21 +0200, Poirier Clement wrote: > > How are you doing this? Via code? > > I just type saveCopy(file, type="emf", etc) after the barplot() > function ; but using windows commands (file, save as) do the same. > > > If so, more than likely you forgot to add: > > > >

[R] ERROR while importing data

Hi everyone, When I try to import data do R, the following message appears: "\U sequences are not supported on Windows". I tried lots of methods to import (read.csv, read.table, RODBC, read.delim) and the same message appears for all these methods. I think it is a bigger problem. Can

[R] String concatenation, File Path Handling to pass to download.file( ) [backslash in DOS paths]

Gabor, Thanks much. Your solution is elegant. My overall scheme is to take present date, and check whether it is a weekend, if not, then create a string based on the date, to concatenate into a url link for download.file( ). The files I need to download have a part which is in the format: mmddyy.

[R] String concatenation, File Path Handling to pass to download.file( ) [backslash in DOS paths]

Gabor, Thanks much. Your solution is elegant. My overall scheme is to take present date, and check whether it is a weekend, if not, then create a string based on the date, to concatenate into a url link for download.file( ). The files I need to download have a part which is in the format: mmddyy.

Re: [R] Bad EMF export

> How are you doing this? Via code? I just type saveCopy(file, type="emf", etc) after the barplot() function ; but using windows commands (file, save as) do the same. > If so, more than likely you forgot to add: > > dev.off() > > after the code that generates the plot. If you don't close

Re: [R] clipping off words inside a vector of strings

On Mon, 2007-10-15 at 12:04 -0400, Gonçalo Ferraz wrote: > Hi, > > I have a vector of strings (class character) with 6 elements (length > 6). I call it 'names'. > > "Graham Chapman" > "John Cleese" > "Terry Gilliam" > "Eric Idle" > "Terry Jones" > "Michael Palin" > > And I want to turn it into

Re: [R] Need some help

Umm, yes, what you had makes a lot of sense. How would I represent that in a plot of the number of sign changes in the to Y axis, and the toss number (from 0 to 1000) in the x-axis? jholtman wrote: > > You might want to check out 'rle'. This will give you the 'lengths' > of runs of the sam

Re: [R] clipping off words inside a vector of strings

Lemon Curry ? sapply( strsplit( monty, " " ), function(x) { paste( substring(x,1,3), collapse = " " ) }) is a way to do it, ... There is probably a better way to do that using the gsubfn package Gonçalo Ferraz wrote: > Hi, > > I have a vector of strings (class character) with 6 elements (leng

[R] clipping off words inside a vector of strings

Hi, I have a vector of strings (class character) with 6 elements (length 6). I call it 'names'. "Graham Chapman" "John Cleese" "Terry Gilliam" "Eric Idle" "Terry Jones" "Michael Palin" And I want to turn it into another vector of strings called 'shortnames' with the same length. The new vect

[R] glmmML vs. lmer - fitting overdispersed Poisson outcome.

Group, I have count data with one observation per subject. I would like to fit a glmm to these data in order to account for overdispersion in the outcome. The lmer() function does not appear to be able to handle data that have only one observation per-cluster id, even though separate variance

[R] Rmpi building

I am trying to build Rmpi on Suse 10.2 linux. While I would like to use the RPM version of liblam-7.1.2, I have not been able to do so. It seems that Rmpi makes some pretty strong assumptions when trying to guess the MPI version. Has anyone successfully built Rmpi against lam-7.1.2 recently usin

Re: [R] Variable which has the maximum value of DF

You can also use arr.ind=TRUE to get the indices of the maximum value: > Lines <- "var1 var2 var3 var4 var5 var6 + + 0 2 1 2 0 0 + 2 3 7 6 0 1 + 1.54 9 9 6 0 + 1.06 1022 3 3 + " > DF <- read.table(text

Re: [R] Need some help

You might want to check out 'rle'. This will give you the 'lengths' of runs of the same value and therefore when the value changes (sign change?) you can see how often: > x <- sample(c(-1,1), 1000, TRUE) > rle(x) Run Length Encoding lengths: int [1:483] 2 2 1 4 3 1 1 1 1 2 ... values : num [1

Re: [R] Bad EMF export

On Mon, 2007-10-15 at 17:36 +0200, Poirier Clement wrote: > Hello dear useRs, > > I'm trying to export a barplot into an emf file. My problem is that > the plot is properly printed into the file, except the bars that do > not appear :( > I've experienced some problems also with simple points p

Re: [R] Bad EMF export

Poirier Clement wrote: > Hello dear useRs, > > I'm trying to export a barplot into an emf file. My problem is that > the plot is properly printed into the file, except the bars that do > not appear :( > I've experienced some problems also with simple points plots, in which > points did not ap

Re: [R] The "condition has length > 1" issue for lists

Your logic test is not correct. Here is what I think you want. See if it makes sense. > x <- list() # create some test data > x[[1]] <- read.table(textConnection("1 2 3 4 + 5 6 7 8 + 4 3 2 1 + 8 7 6 5")) > # create another list element > x[[2]] <- t(x[[1]]) > str(x[[1]]) 'data.frame': 4 obs.

[R] Bad EMF export

Hello dear useRs, I'm trying to export a barplot into an emf file. My problem is that the plot is properly printed into the file, except the bars that do not appear :( I've experienced some problems also with simple points plots, in which points did not appear (same problem). Can you help m

Re: [R] Variable which has the maximum value of DF

Henrique Dallazuanna wrote: > Perhaps, > > names(which.max(sapply(DF, max))) > Nice. I was thinking along the lines of M <- as.matrix(DF) colnames(M)[col(M)[which.max(M)]] > > On 15/10/2007, Lauri Nikkinen <[EMAIL PROTECTED]> wrote: > >> Hi, >> >> Suppose I have a data.frame like this >> >>

Re: [R] some question about partial prediction in survival

On Sun, 14 Oct 2007, coldeyes.Rhelp wrote: > Hi there: > i got a problem to get the prediction from a model recently. for > example if i use a survival analysis to predict the risk. i use the code > like below: i found the the prediction is not equal to (coef * x + coef > * sex) , could someone h

Re: [R] a question on impulse responses

On Sat, Oct 13, 2007 at 03:08:58PM +0300, Martin Ivanov wrote: > > Dear R users, > > I am using the vars package to calculate the impulse response functions and > the forecast error variance decomposition of a VAR model. Unfortunately I do > not know whether these functions assume unit or one s

Re: [R] for loop if else conditional

See ?seq.Date, e.g. now <- Sys.Date() dd <- seq(now - 20, now, by = "day") dd[as.POSIXlt(dd)$wday %% 6 != 0] and have a look at R News 4/1. On 10/15/07, Vishal Belsare <[EMAIL PROTECTED]> wrote: > date <- as.POSIXlt(Sys.time()) #present date > for (i in 1:difftime(as.PO

[R] for loop if else conditional

date <- as.POSIXlt(Sys.time()) #present date for (i in 1:difftime(as.POSIXlt(Sys.Date()),"2007-10-01")) if (date$wday != 0 & date$wday != 6) {print(date);assign("date", (date-86400))} else (assign("date", (date-86400))) I am trying to print dates from present day to a day in the

Re: [R] Need some help

Quite helpful indeed. Greatly appreciated. Another problem I had was trying to simulate an example from my book. Simulating 1000 coin tosses, and finding the frequency of sign changes. So how will we plot this using R? (frequency of sign changes in Y axis) Daniel Nordlund wrote: > >> -Or

[R] for loop if else conditional

date <- as.POSIXlt(Sys.time()) #present date for (i in 1:difftime(as.POSIXlt(Sys.Date()),"2007-10-01")) if (date$wday != 0 & date$wday != 6) {print(date);assign("date", (date-86400))} else (assign("date", (date-86400))) I am trying to print dates from present day to a da

[R] for loop if else conditional

date <- as.POSIXlt(Sys.time()) #present date for (i in 1:difftime(as.POSIXlt(Sys.Date()),"2007-10-01")) if (date$wday != 0 & date$wday != 6) {print(date);assign("date", (date-86400))} else (assign("date", (date-86400))) I am trying to print dates from present day to a day

Re: [R] Data manipulation

Hi [EMAIL PROTECTED] napsal dne 15.10.2007 14:36:59: > 2007/10/15, Stephen Tucker <[EMAIL PROTECTED]>: > > > > Hi Klaus, > > > > I am not exactly sure what you are asking for, but something like this? > > This > > would be option (2) from your list - I don't know that it would be too > > difficu

Re: [R] Fonts do not display properly in R 2.5.1 on Red HatEnterprise Linux 4

OK, tried google, got very, very lost. There are lots of different packages out there. Can anyone tell me where I can download the "Unicode X11 (meta-)fonts" for Red Hat that R needs? -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of michael watson (IAH-C)

Re: [R] Fonts do not display properly in R 2.5.1 on Red Hat Enterprise Linux 4

Thank you Brian, setting the locale using. Sys.setlocale("LC_CTYPE","en_GB") Meant that my test plot command worked fine. Will now install Unicode X11 fonts -Original Message- From: Prof Brian Ripley [mailto:[EMAIL PROTECTED] Sent: 15 October 2007 14:12 To: michael watson (IAH-C) Cc:

[R] Error: X11 cannot allocate additional graphics colours.

Dear All Another one I have touched on before with a much older OS and version. My sessionInfo() is: > sessionInfo() R version 2.5.1 (2007-06-27) i686-redhat-linux-gnu locale: LC_CTYPE=en_GB.UTF-8;LC_NUMERIC=C;LC_TIME=en_GB.UTF-8;LC_COLLATE=en_GB.U TF-8;LC_MONETARY=en_GB.UTF-8;LC_MESSAGES=en_G

Re: [R] Variable which has the maximum value of DF

Perhaps, names(which.max(sapply(DF, max))) On 15/10/2007, Lauri Nikkinen <[EMAIL PROTECTED]> wrote: > > Hi, > > Suppose I have a data.frame like this > > Lines <- "var1 var2 var3 var4 var5 var6 > > 0 2 1 2 0 0 > 2 3 7 6 0 1 > 1.54 9 9

[R] Variable which has the maximum value of DF

Hi, Suppose I have a data.frame like this Lines <- "var1 var2 var3 var4 var5 var6 0 2 1 2 0 0 2 3 7 6 0 1 1.54 9 9 6 0 1.06 1022 3 3 " DF <- read.table(textConnection(Lines), skip=1) names(DF) <-

Re: [R] Need help ploting time series

In the zoo package plot.zoo and xyplot.zoo can do this. library(zoo) vignette("zoo") vignette("zoo-quickref") ?plot.zoo ?xyplot.zoo On 10/14/07, gsmkb86 <[EMAIL PROTECTED]> wrote: > > Hi: > I been having a lot of trouble trying to plot multiple time series on the > same plot. What I want to do is

Re: [R] Fonts do not display properly in R 2.5.1 on Red Hat Enterprise Linux 4

This looks like what happens when you use a UTF-8 locale and don't have Unicode X11 (meta-)fonts installed. Try running in LC_CTYPE=en_GB: if that works you will both have a workaround and know where to look for a solution. On Mon, 15 Oct 2007, michael watson (IAH-C) wrote: > Dear All > > I p

[R] Fonts do not display properly in R 2.5.1 on Red Hat Enterprise Linux 4

Dear All I posted a similar question quite some time ago, but that was on an old OS and an old version of R. This time I have RHEL 4, which is still supported as an OS, and R 2.5.1 which is not *that* old. My sessionInfo() gives: > sessionInfo() R version 2.5.1 (2007-06-27) i686-redhat-linux-g

Re: [R] Data manipulation

2007/10/15, Stephen Tucker <[EMAIL PROTECTED]>: > > Hi Klaus, > > I am not exactly sure what you are asking for, but something like this? > This > would be option (2) from your list - I don't know that it would be too > difficult in R that you would want to use another tool. > > filt <- function(x)

Re: [R] cumulative frequency plots for factors

Perhaps: par(mfrow=c(3,1)) lapply(tapply(df$time, df$events, ecdf), plot) On 15/10/2007, Dieter Vanderelst <[EMAIL PROTECTED]> wrote: > > Dear list, > > I have a data frame with a number of events (factor) and the times at > which they occurred (continuous variable): > > event time > A 10 > A 12

[R] [R-pkgs] new package 'nnls'

A new package 'nnls' is now available on CRAN. The package provides an R interface to the Lawson-Hanson NNLS algorithm for non-negative least squares that solves the least squares problem A x = b with the constraint x >= 0. The Lawson-Hanson NNLS algorithm was published in Lawson CL, Hanson RJ (

[R] cumulative frequency plots for factors

Dear list, I have a data frame with a number of events (factor) and the times at which they occurred (continuous variable): event time A 10 A 12 B 15 A 17 C 13 ... Is it possible in R to make a plot against time of the cumulative frequency of occurrence of each event? This would be, a raising

[R] The "condition has length > 1" issue for lists

I have the following code: list1 <- list() for (i in list.files(pattern="filename1")){ x <- read.table(i) list1[[i]] <- x } list2 <- list() for (i in list.files(pattern="filename2*")){ x <- read.table(i) list2[[i]] <- x } anslist <- vector('list', length(list1)) for(i i

[R] Need help ploting time series

Hi: I been having a lot of trouble trying to plot multiple time series on the same plot. What I want to do is the following: I have a table with sales per day on different states and what i would like to do is plot time series for all the diferents states in the same plot (divided in small squares

Re: [R] Download old R-2.5.1 version?

For Windows: http://cran.r-project.org/bin/windows/base/old/ For Linux: http://cran.r-project.org/bin/linux/ On 15/10/2007, Hans-Peter <[EMAIL PROTECTED]> wrote: > > Is it possible to download somewhere the former 2.5.1 windows bin version > of > R? (Need it for testing). > > Thanks, > Hans-Peter

[R] Download old R-2.5.1 version?

Is it possible to download somewhere the former 2.5.1 windows bin version of R? (Need it for testing). Thanks, Hans-Peter [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLE

[R] iplots question

Are there parameters to set position and size of the windows created by the commands in iplots? I could not find them. __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project

[R] Odp: Save variable with loop part

Hi [EMAIL PROTECTED] napsal dne 15.10.2007 11:12:57: > Hello > > does anyone know how to save in a variable with a loop part. > > Im reading several csv files with read.table and would like to save with a loop: > > for (i in 0:9){ > } > > should give read0, read1, read2, read3, read4, read5,

Re: [R] how to generate and evaluate a design using Algdesign

Thanks Steve, Now it is a bit clear to me, thanks very much. - Original Message - From: "S Ellison" <[EMAIL PROTECTED]> To: <[EMAIL PROTECTED]>; "sun" <[EMAIL PROTECTED]> Sent: Friday, October 12, 2007 12:59 PM Subject: Re: [R] how to generate and evaluate a design using Algdesign >

Re: [R] oanda and yahoo get.hist.quote

Hello Alexander I doubt that such an analyis is very useful as the data is not sampled synchronously (equity close in the US for ^gspc and even that is not always at the same time, some average price from Oanda data). Also fx data from others sources as suggested in another mail on this list wo

[R] Save variable with loop part

Hello does anyone know how to save in a variable with a loop part. Im reading several csv files with read.table and would like to save with a loop: for (i in 0:9){ } should give read0, read1, read2, read3, read4, read5, read6, read7, read8, read9 st. like read(i)<-read.table("myfile.csv".)

  1   2   >