[R] Is there any "month" object like "LETTERS" ?

2009-09-11 Thread megh
There is an object "LETTERS" which displays all letters from "a" to "z". Is there any similar object whicg displays the "months" as well in chronological order? like "jan", "feb",...,"dec" Thanks, -- View this message in context: http://www.nabble.com/Is-there-any-%22month%22-object-lik

Re: [R] Best R text editors?

2009-09-11 Thread Patrick Connolly
On Fri, 11-Sep-2009 at 06:12AM +0200, Johannes Huesing wrote: |> Martin Maechler [Wed, Sep 02, 2009 at 09:17:42AM CEST]: |> > > "PaCo" == p connolly |> > > on Wed, 02 Sep 2009 12:19:31 +1200 writes: |> > |> > PaCo> On Mon, 31-Aug-2009 at 08:25PM +1000, Jim Lemon wrote: |> >

[R] Constructing variables conditional on two indicators

2009-09-11 Thread Reynaerts, Jo
Dear R users I have a data frame that looks like this: m j X1 X2 1 1 1 2 1 ... 1 J 2 1 2 2 2 ... 2 J . . . M 1 M 2 M ... M J The data frame essentially looks like a

Re: [R] Constructing variables conditional on two indicators

2009-09-11 Thread Dimitris Rizopoulos
have a look at ave() and/or tapply(). I hope it helps. Best, Dimitris Reynaerts, Jo wrote: Dear R users I have a data frame that looks like this: m j X1 X2 1 1 1 2 1 ... 1 J 2 1 2 2 2 ... 2 J

Re: [R] Linux R version: best?

2009-09-11 Thread Patrick Connolly
On Thu, 10-Sep-2009 at 04:47PM -0400, S. Few wrote: |> For my Redhat 5.2 Linux box, which version of R would be most stable? Redhat 5.2 is *very* old. That was around about the time of Windows98, IIRC. Or do you mean RHEL 5.2? |> |> I am doing forecasting, statistics, etc. |> |> |> Thanks!

Re: [R] executing rscript from VB

2009-09-11 Thread bartjoosen
Maybe I'm missing something, but how about using the shell function? Bart Duncan Murdoch-2 wrote: > > On 9/10/2009 9:25 AM, H Rao wrote: >> Hi, >> I am looking to execute an R script from VB as below. >> The script runs fine but the redirection doesnt seem to happen. The >> redirection opera

[R] how to do this?

2009-09-11 Thread Luca Braglia
Hello R-users I have a situation like this x=c(1,3,2) y=data.frame("a"=1:3, "b"=4:6, "c"=7:9)*2 So we have > t(t(x)) [,1] [1,]1 [2,]3 [3,]2 And > y a b c 1 2 8 14 2 4 10 16 3 6 12 18 I would like to obtain a vector with number taken from the data.frame: x is needed as

Re: [R] sppolot: fill below minimum legend value

2009-09-11 Thread Paul Hiemstra
emorway wrote: In the plot below, there are some grid cells that have values below 10, which is the lowest "cut" value I have specified. Is there a way, without adjusting the number of cuts, to tell R to fill in those cells with the lowest possible color (in this case greeen)? There is a white

[R] ipred

2009-09-11 Thread MERAL YAY
Hello, I have a question about "ipred" package. I am working on a data set which contains 1000 individual swho came from 588 regions. I calculated misclassification error rate with cross validation and its smooted version "bootstrap 632" for different sample sizes in LDF. I am doubting whether er

[R] call Fortran from R

2009-09-11 Thread Giacomo Santini
Dear R users, I have to call fortran program from within R (R 2.8.1 on ubuntu 8.10 machine). Suppose I have a fortran code like this (this is only a toy model, my working model is far more complex, but input/output is similar) DOUBLE PRECISION FUNCTION model(times, alfa, beta) DOU

[R] (no subject)

2009-09-11 Thread MERAL YAY
Hello, I have a question about "ipred" package. I am working on a data set which contains 1000 individual swho came from 588 regions. I calculated misclassification error rate with cross validation and its smooted version "bootstrap 632" for different sample sizes in LDF. I am doubting whether er

[R] How to compare the result of GLM and GAM

2009-09-11 Thread Dilli Prasad Rijal
Dear R users I have basic knowledge of R and unaware of many more. I am confused about how can I compare the result of two different models? Have count data of plant species and want to correlate with altitude, rri, pH, moisture, temperature etc. I have made some models using GLM of different ord

[R] R - box design-scatter plot für means/regres sion/lme?

2009-09-11 Thread Karrer Stefanie
Dear All! It's now weeks that I'm going crazy with R, and as I'm a new user I now ask for help (also because I still have only a few days to finish..)... So shortly I describe you my Experiment in which I was looking for the decomposition of herbivore dung under different treatments: I made a b

Re: [R] call Fortran from R

2009-09-11 Thread Duncan Murdoch
On 11/09/2009 4:39 AM, Giacomo Santini wrote: Dear R users, I have to call fortran program from within R (R 2.8.1 on ubuntu 8.10 machine). Suppose I have a fortran code like this (this is only a toy model, my working model is far more complex, but input/output is similar) DOUBLE PREC

[R] R: how to do this?

2009-09-11 Thread Luca Braglia
> Da: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] > Per conto di Luca Braglia > Inviato: venerdì 11 settembre 2009 10.17 > A: r-help@r-project.org > Oggetto: [R] how to do this? > in this case the ouput should be > > 2 > 16 > 12 > Obviously this id.match = rep(0,

Re: [R] call Fortran from R

2009-09-11 Thread Giacomo Santini
Thanks! now it works ! Giacomo Duncan Murdoch wrote: On 11/09/2009 4:39 AM, Giacomo Santini wrote: Dear R users, I have to call fortran program from within R (R 2.8.1 on ubuntu 8.10 machine). Suppose I have a fortran code like this (this is only a toy model, my working model is far more c

Re: [R] how to do this?

2009-09-11 Thread Linlin Yan
Try this: > y[matrix(c(seq_along(x), x), ncol = 2)] [1] 2 16 12 On Fri, Sep 11, 2009 at 4:17 PM, Luca Braglia wrote: > Hello R-users > > I have a situation like this > > x=c(1,3,2) > > y=data.frame("a"=1:3, "b"=4:6, "c"=7:9)*2 > > So we have > >> t(t(x)) >     [,1] > [1,]    1 > [2,]    3 > [3,]

[R] Multilevel models with sampling weights at both levels

2009-09-11 Thread David Kaplan
Greetings, Is there a package in R that will run multilevel models (e.g. students nested in schools) where sampling weights can be employed at both levels? Thanks in advance. David -- === David Kaplan, Ph.D. Professor Department of

[R] R: how to do this?

2009-09-11 Thread Luca Braglia
> Da: Linlin Yan [mailto:yanlinli...@gmail.com] > Inviato: venerdì 11 settembre 2009 11.19 > A: Luca Braglia > Cc: r-help@r-project.org > Oggetto: Re: [R] how to do this? > > Try this: > > y[matrix(c(seq_along(x), x), ncol = 2)] > [1] 2 16 12 Very interesting! So I can give as index to a data

Re: [R] Is there any "month" object like "LETTERS" ?

2009-09-11 Thread Barry Rowlingson
On Fri, Sep 11, 2009 at 8:13 AM, megh wrote: > > There is an object "LETTERS" which displays all letters from "a" to "z". Is > there any similar object whicg displays the "months" as well in > chronological order? like "jan", "feb",...,"dec" You could construct a vector of the first of t

Re: [R] How to do rotation for polygon?

2009-09-11 Thread baptiste auguie
Hi, See if this helps, polygon.regular <- # return a matrix of xy coordinates for a regular polygon centered about (0,0) function (sides = 5) { n <- sides if (n < 3) n <- 3 if (n > 8) n <- 50 th <- pi * 2/n costh <- cos(th) sinth <- sin(th)

Re: [R] Best R text editors?

2009-09-11 Thread Jim Lemon
On 09/11/2009 05:15 PM, Patrick Connolly wrote: ... |> and in previous versions, you could always do M-x cua-mode for |> the same effect. Talk about a well-hidden function mostly directed |> at beginners ... Perhaps the thinking was that by the time they find it, they'll already have noticed

Re: [R] Simple time series questions

2009-09-11 Thread DKOD
This script worked for me. Be sure to put in your correct link. link <- "C:\\R_Home\\Charts & Graphs Blog\\R_Chart_Doc\\text_data.csv" testdata<- read.table(link, head = T, sep = ",",na.strings = "na") test_date = as.Date(testdata$Date,"%d-%m-%y") plot(test_date, testdata$Model, type="l

Re: [R] Best R text editors?

2009-09-11 Thread Ted Harding
On 11-Sep-09 10:41:21, Jim Lemon wrote: > On 09/11/2009 05:15 PM, Patrick Connolly wrote: >> ... >> |> and in previous versions, you could always do M-x cua-mode for >> |> the same effect. Talk about a well-hidden function mostly directed >> |> at beginners ... >> >> Perhaps the thinking was tha

Re: [R] Best R text editors?

2009-09-11 Thread Duncan Murdoch
On 11/09/2009 6:53 AM, (Ted Harding) wrote: On 11-Sep-09 10:41:21, Jim Lemon wrote: On 09/11/2009 05:15 PM, Patrick Connolly wrote: ... |> and in previous versions, you could always do M-x cua-mode for |> the same effect. Talk about a well-hidden function mostly directed |> at beginners ...

Re: [R] Is there any "month" object like "LETTERS" ?

2009-09-11 Thread Gabor Grothendieck
See ?Constants On Fri, Sep 11, 2009 at 3:13 AM, megh wrote: > > There is an object "LETTERS" which displays all letters from "a" to "z". Is > there any similar object whicg displays the "months" as well in > chronological order? like "jan", "feb",...,"dec" > > Thanks, > -- > View this mes

Re: [R] Color index in image function

2009-09-11 Thread Bernardo Rangel Tura
On Thu, 2009-09-10 at 05:27 -0700, FMH wrote: > Thank you for the scripts, but the label and the values in the x and y-axis > suddently dissapear even the 'axis' function is used as stated in the command > below. Could you help on this? > > axis(1, at = seq(100, 800, by = 100)) > axis(2, at = se

Re: [R] how to do this?

2009-09-11 Thread Henrique Dallazuanna
Try this: y[cbind(seq(length(x)), x)] On Fri, Sep 11, 2009 at 5:17 AM, Luca Braglia wrote: > Hello R-users > > I have a situation like this > > x=c(1,3,2) > > y=data.frame("a"=1:3, "b"=4:6, "c"=7:9)*2 > > So we have > > > t(t(x)) > [,1] > [1,]1 > [2,]3 > [3,]2 > > And > > > y >

Re: [R] Issue displaying DATES on a plot with two ordinates

2009-09-11 Thread Jim Lemon
Hi Clair, Good one for pointing out that dates on the x-axis are messy (and for bringing my mistake in the coloring of the ordinates to my attention). Attached is a revised version of the function that will be in the next version of plotrix. Try this: twoord.plot(x.Left, y.Left, x.Right, y.Ri

[R] What determines the unit of POSIXct differences?

2009-09-11 Thread Heinz Tuechler
Dear All, what determines if a difference between POSIXct objects gets expressed in days or seconds? In the following example, it's sometimes seconds, sometimes days. as.POSIXct('2009-09-01') - as.POSIXct(NA) Time difference of NA secs c(as.POSIXct('2009-09-01'), as.POSIXct(NA)) - c(as.POSI

Re: [R] Merge data frames but prefer values in one

2009-09-11 Thread Henrique Dallazuanna
Maybe: do.call(rbind, lapply(with(xy <- rbind(x, y), split(xy, list(a, b), drop = TRUE)), tail, 1)) On Fri, Sep 11, 2009 at 3:45 AM, jo wrote: > Thanks for the post-processing ideas. But is there any way to do that > in one step? > > On Thu, Sep 10, 2009 at 7:20 PM, Henrique Dallazuanna > wrot

[R] Graph visualization

2009-09-11 Thread Arber Ngjela
Hello, I am working with graph and adjacency matrix, the package 'graph' seems to be appropriate for this. An example in the package > mat <- rbind(c(0, 0, 1, 1), + c(0, 0, 1, 1), + c(1, 1, 0, 1), + c(1, 1, 1, 0)) > rownames(mat) <- colnames(mat) <- letters

Re: [R] Plot area reduction

2009-09-11 Thread jim holtman
you can use the 'mar' parameter of 'par' to setup the margins of the graph and then use the height/width option when opening the graphic device to control the plotting area. On Thu, Sep 10, 2009 at 6:10 AM, rajesh j wrote: > Hi, > > I need my plot to occupy a thin strip-like area but the plot are

[R] Searching 1st number within a ch. string

2009-09-11 Thread RON70
Hi all, Is there any process to conduct a search for a particular digit or letter in a ch. string? For example I want to make a search where 1st numeric figure appears in the string "asd123". Here the answer should be "4" Thanks, -- View this message in context: http://www.nabble.com/Searching

Re: [R] Searching 1st number within a ch. string

2009-09-11 Thread Henrique Dallazuanna
Try this: head(unlist(gregexpr("[0-9]", "asd123")), 1) On Fri, Sep 11, 2009 at 10:02 AM, RON70 wrote: > > Hi all, > > Is there any process to conduct a search for a particular digit or letter > in > a ch. string? For example I want to make a search where 1st numeric figure > appears in the stri

Re: [R] Simple time series questions

2009-09-11 Thread DKOD
Try this script. I converted test_date to numeric decimal year link <- "C:\\R_Home\\Charts & Graphs Blog\\R_Chart_Doc\\text_data.csv" testdata<- read.table(link, head = T, sep = ",",na.strings = "na") test_date = as.Date(testdata$Date,"%d-%m-%y") # Convert dates to decimal year my_yr <-

Re: [R] What determines the unit of POSIXct differences?

2009-09-11 Thread jim holtman
'-' calls 'difftime' which, if you don't specify the units, makes the following assumptions in the code: > difftime function (time1, time2, tz = "", units = c("auto", "secs", "mins", "hours", "days", "weeks")) { time1 <- as.POSIXct(time1, tz = tz) time2 <- as.POSIXct(time2, tz = tz)

Re: [R] Searching 1st number within a ch. string

2009-09-11 Thread Gabor Grothendieck
Try this: > regexpr("[0-9]", "asd123") [1] 4 attr(,"match.length") [1] 1 On Fri, Sep 11, 2009 at 9:02 AM, RON70 wrote: > > Hi all, > > Is there any process to conduct a search for a particular digit or letter in > a ch. string? For example I want to make a search where 1st numeric figure > appe

Re: [R] What determines the unit of POSIXct differences?

2009-09-11 Thread Heinz Tuechler
Jim - Thank you very much for this explanation and the hint to use difftime. Heinz At 15:09 11.09.2009, jim holtman wrote: '-' calls 'difftime' which, if you don't specify the units, makes the following assumptions in the code: > difftime function (time1, time2, tz = "", units = c("auto", "sec

[R] : How wo read stability VAR plot?

2009-09-11 Thread Arif Chandra
I have made program code for Vector Auto Regressive in terms of completing my undergraduate program using R. I have an important question related to my project. If I have: data(Canada) var.2c <- VAR(Canada, p = 2, type = "const") var.2c.stabil <- stability(var.2c, type = "O

[R] Sorting

2009-09-11 Thread megh
I have following object : > date2 [,1] [,2] [1,] "apr" "1992" [2,] "aug" "1992" [3,] "dec" "1992" [4,] "feb" "1992" [5,] "jan" "1992" [6,] "jul" "1992" [7,] "jun" "1992" [8,] "mar" "1992" [9,] "may" "1992" [10,] "nov" "1992" [11,] "oct" "1992" [12,] "sep" "1992" [1

Re: [R] Sorting

2009-09-11 Thread Henrique Dallazuanna
Try this: data2[order(data2[,2], match(data2[,1], tolower(month.abb))),] On Fri, Sep 11, 2009 at 10:46 AM, megh wrote: > > I have following object : > > > date2 > [,1] [,2] > [1,] "apr" "1992" > [2,] "aug" "1992" > [3,] "dec" "1992" > [4,] "feb" "1992" > [5,] "jan" "1992" > [6,] "j

[R] Modify functions in base packages (needed for completion for proto objects)

2009-09-11 Thread Vitalie S.
Hello everyone, I am trying to implement completion for proto objects. Proto extends environment in a hierarchical way. Thus completion should list all the names in all it's parent environments. For "normal" classes defining names.class would do the job, but completion for "environment"

Re: [R] Sorting

2009-09-11 Thread jim holtman
You might want to add a date column it you are planning to operate on the data with dates: > a # your data X..1. X..2. 1apr 1992 2aug 1992 3dec 1992 4feb 1992 5jan 1992 6jul 1992 7jun 1992 8mar 1992 9may 1992 10 nov 1992 11 oct 1992 12 sep

Re: [R] Simple time series questions

2009-09-11 Thread Gabor Grothendieck
Here it is using zoo and classic graphics. 1. Just paste this into your R session: library(zoo) URL <- "http://www.nabble.com/file/p25398419/test%2Bchart%2Bdata.csv"; z <- read.zoo(URL, header = TRUE, format = "%d-%m-%y", sep = ",") cols <- c("green", "red", "blue") plot(z, screen = 1, col = c

Re: [R] Modify functions in base packages (needed for completion for proto objects)

2009-09-11 Thread Gabor Grothendieck
See ?assignInNamespace On Fri, Sep 11, 2009 at 10:02 AM, Vitalie S. wrote: > > Hello everyone, > > I am trying to implement completion for proto objects. Proto extends > environment in a hierarchical way. Thus completion should list all the names > in all it's parent environments. > > For "normal

Re: [R] Graph visualization

2009-09-11 Thread Steve Lianoglou
Hi, On Sep 11, 2009, at 8:22 AM, Arber Ngjela wrote: Hello, I am working with graph and adjacency matrix, the package 'graph' seems to be appropriate for this. An example in the package mat <- rbind(c(0, 0, 1, 1), + c(0, 0, 1, 1), + c(1, 1, 0, 1), + c

Re: [R] Best R text editors?

2009-09-11 Thread Clint Bowman
On Fri, 11 Sep 2009, Duncan Murdoch wrote: On 11/09/2009 6:53 AM, (Ted Harding) wrote: On 11-Sep-09 10:41:21, Jim Lemon wrote: > On 09/11/2009 05:15 PM, Patrick Connolly wrote: > > ... > > | > and in previous versions, you could always do M-x cua-mode for > > | > the same effect. Talk ab

[R] voice recognition software

2009-09-11 Thread volinsky
Hello all. Has anyone out there had experience using R with voice recognition software? I got this query from a student with nerve damage in her hand who wants to program in R. Thanks, Chris __ R-help@r-project.org mailing list https://stat.ethz

Re: [R] Mapping factors to a new set of factors

2009-09-11 Thread james.arnold
Thanks Jorge, That is what I was looking for. Cheers, James From: Jorge Ivan Velez [mailto:jorgeivanve...@gmail.com] Sent: 08 September 2009 18:41 To: Arnold, James Cc: r-help@r-project.org Subject: Re: [R] Mapping factors to a new set of factors Hi James, Take a look at the "recode" function

Re: [R] Mantel test least square line

2009-09-11 Thread swertie
Thank you it is what I was looking for. -- View this message in context: http://www.nabble.com/Mantel-test-least-square-line-tp25235402p25401329.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org mailing list https:

[R] Reading gzip data from a non zero file offset

2009-09-11 Thread Gregory Jefferis
Dear R users, I have a file that contains a text header followed by a gzipped data chunk. I thought I would be able to read this by doing something like this. con<-file('myfile','rb') header<-readLines(con,10) gzf<-gzcon(con) data<-readBin(gzf,"int",n=1e7) But what I find is that gzcon resets

Re: [R] Best R text editors?

2009-09-11 Thread Ted Harding
On 11-Sep-09 14:16:44, Clint Bowman wrote: > On Fri, 11 Sep 2009, Duncan Murdoch wrote: > >> On 11/09/2009 6:53 AM, (Ted Harding) wrote: >>> On 11-Sep-09 10:41:21, Jim Lemon wrote: >>> > On 09/11/2009 05:15 PM, Patrick Connolly wrote: >>> > > ... >>> > > | > and in previous versions, you coul

[R] format (?) problems with data imported from postgres

2009-09-11 Thread Tomas Lanczos
Good day, I read some data from a PostgreSQL database by a following script: library(Rdbi) library(RdbiPgSQL) # conn becomes an object which contains the DB connection: conn <- dbConnect(PgSQL(), host="localhost", dbname="BVS", user="postgres", password = "***") query_duj_kal <- dbSendQuery(

[R] RPostgreSQL package and libpq.dll file

2009-09-11 Thread Lore M
Dear all, I'd like to use the package RPostgreSQL. I'm using R version 2.8.1 and I've download the last version RPostgreSQL. When I load the package, I get something like "the file LIBPQ.DLL is missing". Do you have any idea about what I'm suppose to do ? Thanks everyone. __

Re: [R] Mapping factors to a new set of factors

2009-09-11 Thread james.arnold
Thanks a lot Phil, Recode is exactly what I was looking for. I managed to get my old function working using sapply, but the performance was horrendously slow! One other thing was that the lvls vector can only seem to be set within the global scope of R, and local variables within a function do

Re: [R] Mapping factors to a new set of factors

2009-09-11 Thread james.arnold
Thanks Jim, I managed to get my function to work by putting it within an sapply. Unfortunately the performance was terrible, so I've gone with the suggestion of other posters to use recode in the car package. Cheers, James -Original Message- From: jim holtman [mailto:jholt...@gmail.com

Re: [R] Best R text editors?

2009-09-11 Thread Ramon Diaz-Uriarte
Ted, I think I share your feelings about mice (e.g., that is why I use window managers where the mouse is not needed or is actually discouraged) but ... On Fri, Sep 11, 2009 at 4:46 PM, Ted Harding wrote: > On 11-Sep-09 14:16:44, Clint Bowman wrote: >> On Fri, 11 Sep 2009, Duncan Murdoch wrote: >

[R] State Space models in R

2009-09-11 Thread Giovanni Petris
Hello everybody, I am writing a review paper about State Space models in R, and I would like to cover as many packages as I reasonably can. So far I am familiar with the following tools to deal with SS models: * StructTS, Kalman* (in stats) * packages dse[1-2] * package sspir * package dlm I

Re: [R] Is there any "month" object like "LETTERS" ?

2009-09-11 Thread Don MacQueen
On my R 2.9.2 installation, ?Constants just says "No documentation ...". I don't know why Gabor and I would have different results. However, ?LETTERS does provide the requested information. -Don At 7:30 AM -0400 9/11/09, Gabor Grothendieck wrote: See ?Constants On Fri, Sep 11, 2009 at 3:13 A

Re: [R] RPostgreSQL package and libpq.dll file

2009-09-11 Thread Joe Conway
Lore M wrote: > Dear all, I'd like to use the package RPostgreSQL. I'm using R > version 2.8.1 and I've download the last version RPostgreSQL. When I > load the package, I get something like "the file LIBPQ.DLL is > missing". Do you have any idea about what I'm suppose to do ? Thanks > everyone. l

Re: [R] Is there any "month" object like "LETTERS" ?

2009-09-11 Thread William Dunlap
> -Original Message- > From: r-help-boun...@r-project.org > [mailto:r-help-boun...@r-project.org] On Behalf Of Don MacQueen > Sent: Friday, September 11, 2009 8:16 AM > To: Gabor Grothendieck; megh > Cc: r-help@r-project.org > Subject: Re: [R] Is there any "month" object like "LETTERS" ? >

Re: [R] format (?) problems with data imported from postgres

2009-09-11 Thread Tomas Lanczos
Being more specific: On Fri, 2009-09-11 at 17:03 +0200, Tomas Lanczos wrote: > Good day, > > I read some data from a PostgreSQL database by a following script: > > library(Rdbi) > library(RdbiPgSQL) > # conn becomes an object which contains the DB connection: > conn <- dbConnect(PgSQL(), host="l

Re: [R] R 2.9.2 memory max - object vector size

2009-09-11 Thread gug
At the risk of stating the obvious: - rm(.) # clears specific objects out of memory as soon as they're no longer needed in the routine. - sapply(ls(), function(x) object.size(get(x))) #lists all objects with the memory each is using. - rm(list=ls()) #clears out all objects, e.g. be

Re: [R] Simple time series questions

2009-09-11 Thread gug
Thanks - that works great. Do you have any suggestions about the grid() problem - i.e. that the vertical gridlines do not line up with the x-axis tickmarks (which are years)? I can't see on what basis the vertical gridlines are being positioned, but it doesn't look good that they are not lined u

[R] Rmetrics timeDate - business days between dates

2009-09-11 Thread Leonardo Miceli
Hi One of the most important calculation in applied finance is the number of days between dates. That kind of calculus become annoying when a specific calendar must be used. That is the case for the business days calculus. The package timeDate has a function isBizday to perform that kind of thing

Re: [R] : How wo read stability VAR plot?

2009-09-11 Thread Pfaff, Bernhard Dr.
> > >I have made program code for Vector Auto Regressive in terms >of completing my undergraduate program using R. I have an important >question related to my project. >If I have: >data(Canada) >var.2c <- VAR(Canada, p = 2, type = "const") >var.2c.stabil <- stability(var.2c,

Re: [R] Graph visualization

2009-09-11 Thread Martin Morgan
Arber Ngjela wrote: > Hello, > I am working with graph and adjacency matrix, the package 'graph' seems to be > appropriate for this. > An example in the package >> mat <- rbind(c(0, 0, 1, 1), > + c(0, 0, 1, 1), > + c(1, 1, 0, 1), > + c(1, 1, 1, 0)) >> rowna

[R] Slow enumeration of vectors with restrictions

2009-09-11 Thread Bryan Keller
I have a vector T and a scalar m such that 1 <= m <= sum(T). I also have a vector U that satisfies three properties... 1. length(T)==length(U) 2. sum(U)==m 3. for (i in 1:length(T)) U[i]<=T[i] The function "nextu" (given below) calculates the "next" vector subject to the same restrictions.

Re: [R] State Space models in R

2009-09-11 Thread Gabor Grothendieck
You might try: RSiteSearch("state space") and try the search engine at rseek.org On Fri, Sep 11, 2009 at 11:04 AM, Giovanni Petris wrote: > > Hello everybody, > > I am writing a review paper about State Space models in R, and I would > like to cover as many packages as I reasonably can. > > So

Re: [R] State Space models in R

2009-09-11 Thread spencerg
Hi, Giovanni: 1. I've used primarily "sspir" and "dlm". There is apparently a new release of "sspir", which I have not used. When I last tried "sspir", it did NOT have a forecast function, while "dlm" did. The two packages have functions with the same name but incompatible code. It

Re: [R] Is there any "month" object like "LETTERS" ?

2009-09-11 Thread Steve Murray
month.abb should do the trick _ View your other email accounts from your Hotmail inbox. Add them now. __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEA

[R] Simple time series questions

2009-09-11 Thread gug
I'm sure this is a really simple problem, but I've spent hours digging and I keep running into roadblocks. I'm trying to get a simple chart with three time series. Similar to the attached example http://www.nabble.com/file/p25398419/Excel%2Bchart%2Bexample.pdf Excel+chart+example.pdf , somethin

Re: [R] Completion for proto objects

2009-09-11 Thread Vitalie S.
On Fri, 11 Sep 2009 16:07:45 +0200, Gabor Grothendieck wrote: See ?assignInNamespace thanks for the hint. It works to some extent. Functions are redefined in "base" but utils:::specialCompletions seems not to preform properly: ls.proto <- function(x){ # function to complete .local

Re: [R] transposing a distance matrix in R

2009-09-11 Thread Steve Lianoglou
Hi Jeannine, I'm just forwarding this Q&A back to the r-help list, you'll get more eyes on it and other people might have better solutions. Answers inline: On Sep 11, 2009, at 12:25 PM, Jeannine Cavender-Bares wrote: Dear Steve, Greetings! You helped me earlier this summer with an R quest

Re: [R] Completion for proto objects

2009-09-11 Thread Gabor Grothendieck
You might check the bugs section on the proto home page and if any of those seem to apply try the devel version of proto (or the workaround mentioned) as all known bugs are fixed in the devel version. If that does not help let me know offline and I will try to help you. See: http://r-proto.google

[R] how to determine if a variable is already set?

2009-09-11 Thread carol white
Hi, It might be a primitive question but how it is possible to determine if a variable is initialized in an environment? Suppose that we start a R session and wants to run a script which use the variable i. Which function could evaluate if i is already initialized or not and if not, then ask int

[R] Moving to Mac OS X

2009-09-11 Thread Victor Manuel Garcia Guerrero
Hi all, I have been using R on Windows for a long time, but now I am ready to move to Mac OS X, but I was just wondering if R works better on Mac than on Windows. This is because I had some memory size issues in the past. Another question is if some of you know an R editor for Mac (just like

[R] Working with large matrix

2009-09-11 Thread A Ezhil
Dear All, I have large matrix (46000 x 11250). I would like to do the linear regression for each row. I wrote a simple function that has lm() and used apply(mat,1,func). The issue is that it takes ages to load the file and also to finish the lm. I am using LINUX 64 bit with 32G mem. Is there an

Re: [R] how to determine if a variable is already set?

2009-09-11 Thread Marc Schwartz
On Sep 11, 2009, at 12:15 PM, carol white wrote: Hi, It might be a primitive question but how it is possible to determine if a variable is initialized in an environment? Suppose that we start a R session and wants to run a script which use the variable i. Which function could evaluate if i

[R] how to determine if a variable has already been initialized?

2009-09-11 Thread carol white
Hi, It might be a primitive question but how it is possible to determine if a variable is initialized in an environment? Suppose that we start a R session and wants to run a script which uses the variable i. Which function could evaluate if i exists already and has been initialized or not? The p

Re: [R] transposing a distance matrix in R

2009-09-11 Thread Henrique Dallazuanna
Try this also: subset(as.data.frame.table(m), !is.na(Freq)) On Fri, Sep 11, 2009 at 1:48 PM, Steve Lianoglou < mailinglist.honey...@gmail.com> wrote: > Hi Jeannine, > > I'm just forwarding this Q&A back to the r-help list, you'll get more eyes > on it and other people might have better solutions

Re: [R] AIC and goodness of prediction - was: Re: goodness of "prediction" using a model (lm, glm, gam, brt,

2009-09-11 Thread Kingsford Jones
Hi Corrado, Not being familiar with your research goals or data I can't make recommendations, but I can suggest a couple of places to look for information: Frank Harrell's Regression Modeling Strategies and his Design library available on CRAN, and Hastie et al's The Elements of Statistical Learn

[R] R/S Programmer Employment Opportunity - New York, NY

2009-09-11 Thread Wesley Madison
R / S Programmer - New York Description: Kaplan Test Prep & Admissions is looking for a R / S Programmer to join their research based in New York City. We are looking for highly motivated individuals to work in a customer-focused environment. This is a unique opportunity to develop with a leading

Re: [R] Moving to Mac OS X

2009-09-11 Thread Marc Schwartz
On Sep 11, 2009, at 12:15 PM, Victor Manuel Garcia Guerrero wrote: Hi all, I have been using R on Windows for a long time, but now I am ready to move to Mac OS X, but I was just wondering if R works better on Mac than on Windows. This is because I had some memory size issues in the past.

Re: [R] Moving to Mac OS X

2009-09-11 Thread Steve Lianoglou
Hi, On Sep 11, 2009, at 1:15 PM, Victor Manuel Garcia Guerrero wrote: Hi all, I have been using R on Windows for a long time, but now I am ready to move to Mac OS X, but I was just wondering if R works better on Mac than on Windows. This is because I had some memory size issues in the pa

[R] help with plotting

2009-09-11 Thread Edward Chen
HI all, raw_urine = read.table("Z:\\bruce.9.3.09.sample.stability.analysis\\urine\\mz.spot.sam.dat.new", header = TRUE ) pvalue = read.table("Z:\\bruce.9.3.09.sample.stability.analysis\\urine\\all.urine.features.t.test.result", header = TRUE ) library(compositions) p = function(a,b){ y = pvalue[,a]

[R] Barplot+Table

2009-09-11 Thread Andy Choens
I am trying to automate a report that my company does every couple of years for the state of Maine. In the past we have used SPSS to run the data and then used complicated Excel template to make the tables/graphics which we then imported into Word. Since there are 256 tables/graphics for this re

Re: [R] Working with large matrix

2009-09-11 Thread Steve Lianoglou
Hi, On Sep 11, 2009, at 1:15 PM, A Ezhil wrote: Dear All, I have large matrix (46000 x 11250). I would like to do the linear regression for each row. I wrote a simple function that has lm() and used apply(mat,1,func). The issue is that it takes ages to load the file and also to finish th

Re: [R] help with plotting

2009-09-11 Thread Steve Lianoglou
Hi Edward, On Sep 11, 2009, at 1:47 PM, Edward Chen wrote: HI all, raw_urine = read.table("Z:\\bruce.9.3.09.sample.stability.analysis\\urine\\mz.spot.sam.dat.new ", header = TRUE ) pvalue = read.table("Z:\\bruce.9.3.09.sample.stability.analysis\\urine\\all.urine.features.t.test.result ", he

[R] How to block data across multiple columns?

2009-09-11 Thread A Singh
Dear all, Does anyone have any suggestions on how to block multiple columns of data one at a time in the midst of an analysis, having specified the blocking variable? I am running a random effects model using lmer, and my data set has multiple columns. Individuals in the study are grouped i

Re: [R] Barplot+Table

2009-09-11 Thread Henrique Dallazuanna
Try the textplot function in the gplots package: nf <- layout(matrix(1:2), height = lcm(12)) par(mar = c(2, 8, 2, 2)) barplot(prop.table(example.table, 1), beside = TRUE) textplot(round(prop.table(example.table, 1), 2), halign = 'left', mar = c(0,0,0,0), show.colnames = FALSE, cex = 2,

[R] R help:

2009-09-11 Thread Matthew Fantle
Hi I have written a code to do some averaging of data over uneven intervals. The for loop keeps missing particular depths and I once got an error message reading: *** caught segfault *** address 0xc023, cause 'memory not mapped' Possible actions: 1: abort (with core dump, if enabled) 2

[R] How to Label Certain Lags for a PACF Graph

2009-09-11 Thread Stephanie Cooke
When I use the command for PACF, lags 5, 10, 15, and 20 are labeled. I would like to label lag 1. I would greatly appreciate if someone could tell me how to do this. Below is the command that I am using: pacf(data$R1,main="Series R1 Residuals") [[alternative HTML version deleted]] __

Re: [R] Matrix multiplication and random numbers

2009-09-11 Thread RFish
Hi Sorry I don't seem to have explained what I'm trying to do very clearly. The piece of code below multiplies the two matrices together a number of times based on the value in the matmult(InitialPop,1) term in this case one (year), this gives me the end population for the analysis. InitialPop<

[R] NLME

2009-09-11 Thread Jeff in NC
I am trying to use one of the variance function classes (e.g., varClasses -varConstPower or varPower) When I run the model under default conditions or specify varIdent explicitly I do not have any problems. If I try to use anything else I get error messages of the form "Error in na.fail.defaul

Re: [R] How to Label Certain Lags for a PACF Graph

2009-09-11 Thread Jorge Ivan Velez
Hi Stephanie, Take a look at the "xaxt" argument under ?par as well as ?axis. Here is an example modified from ?pacf : acf(lh, xaxt = "n") axis(1, 0:16, 0:16, cex.axis = 0.8) HTH, Jorge On Fri, Sep 11, 2009 at 12:02 PM, Stephanie Cooke wrote: > When I use the command for PACF, lags 5, 10,

Re: [R] How to Label Certain Lags for a PACF Graph

2009-09-11 Thread Henrique Dallazuanna
Try this: pacf(data$R1, main = "Series R1 Residuals", xaxt = 'n') axis(1, c(1, axTicks(1))) On Fri, Sep 11, 2009 at 1:02 PM, Stephanie Cooke wrote: > When I use the command for PACF, lags 5, 10, 15, and 20 are labeled. I > would > like to label lag 1. I would greatly appreciate if someone could

Re: [R] Moving to Mac OS X

2009-09-11 Thread Noah Silverman
Hi, I'm a daily user of both mac and Linux so wanted to offer some thoughts: 1) R runs great on a Mac. There is a standard install from the cran website that has a nice GUI built into it. You can do things like drag files to the console and it will fill in the path name. 2) I like using B

Re: [R] Matrix multiplication and random numbers

2009-09-11 Thread Henrique Dallazuanna
If I understand you can use replicate: replicate(10, matmult(InitialPop, 1)) On Fri, Sep 11, 2009 at 1:11 PM, RFish wrote: > > Hi > > Sorry I don't seem to have explained what I'm trying to do very clearly. > The piece of code below multiplies the two matrices together a number of > times base

  1   2   >