[R] alternatives to mosiac plots and ?3D? mosiac plots

2009-11-27 Thread DispersionMap
i had some great success with mosiacplots... http://stat.ethz.ch/R-manual/R-patched/library/graphics/html/mosaicplot.html http://stat.ethz.ch/R-manual/R-patched/library/graphics/html/mosaicplot.html ...does anyone know of any other hidden gems out there that allow you to examine the relationshi

Re: [R] Debugging the "unexpected string constant" error in the zp2ssg2.R

2009-11-27 Thread Jason Rupert
Evidently it just needed a little more tweaking... http://n2.nabble.com/Need-help-solving-the-unexpected-string-constant-error-td4077984.html#a4078425 Hopefully all this will work out... - Original Message From: Jason Rupert To: R-help@r-project.org Sent: Fri, November 27, 2009 6:1

Re: [R] [R-es] R en Windows 7 a 64 bits

2009-11-27 Thread Jorge Ivan Velez
My apologies to all of you. I sent this email to the wrong email list :) Best, Jorge 2009/11/27 Jorge Ivan Velez <> > Buenas noches César, > > Tal y como se estipula en la sección 8.1 de [1], la R 32-bits funciona en > Windows 64-bits aunque con las limitaciones de RAM que todos conocemos (sol

Re: [R] How to compute Rolling analysis of Standard Deviation using ZOO package?

2009-11-27 Thread Saji Ren
yeah, my mistake is that i did not make the data "c" into a zoo class. now the problem solved. thank you everybody 2009/11/27 Gabor Grothendieck > We will need a minimal reproducible example (as per last line on every > message to r-help) to answer as trying it with made up data seems to > work

Re: [R] my failing understanding ...

2009-11-27 Thread David Winsemius
On Nov 27, 2009, at 9:39 PM, David Winsemius wrote: I get a different result, undoubtedly because I have different packages loaded: > showMethods(plot) Function "plot": > sessionInfo() R version 2.10.0 Patched (2009-10-29 r50258) x86_64-apple-darwin9.8.0 locale: [1] en_US.UTF-8/en_US.UTF

Re: [R] my failing understanding ...

2009-11-27 Thread David Winsemius
I get a different result, undoubtedly because I have different packages loaded: > showMethods(plot) Function "plot": > sessionInfo() R version 2.10.0 Patched (2009-10-29 r50258) x86_64-apple-darwin9.8.0 locale: [1] en_US.UTF-8/en_US.UTF-8/C/C/en_US.UTF-8/en_US.UTF-8 attached base packages

Re: [R] How to combine matrix with sapply function?

2009-11-27 Thread Ben Bolker
Peng Yu gmail.com> writes: > > I want an 'apply' function that can give me results by bind the result > from each function call of g. Could somebody let me know how to do it? > > g<-function(x){ [snip] > } > > x=1:3 > > sapply(x,g) do.call(rbind,sapply(x,g))?

Re: [R] Export kde object as shapefile

2009-11-27 Thread David Winsemius
I'm not a programmer either, if by "programmer" you mean one who took a computer science degree and is employed with some title that implies a creator of code. I am a user of software toward particular purposes. The purposes need to be defined in sufficiently specific manner to guide that

Re: [R] [R-es] R en Windows 7 a 64 bits

2009-11-27 Thread Jorge Ivan Velez
Buenas noches César, Tal y como se estipula en la sección 8.1 de [1], la R 32-bits funciona en Windows 64-bits aunque con las limitaciones de RAM que todos conocemos (solo puedes usar hasta 4 GB de RAM en el mejor de los casos). Para instalarlo sólo descarga el ejecutable usual de CRAN, haz doble

[R] How to combine matrix with sapply function?

2009-11-27 Thread Peng Yu
I want an 'apply' function that can give me results by bind the result from each function call of g. Could somebody let me know how to do it? g<-function(x){ if(x==1){ rbind(c(1,3,8)) } else if(x==2) { rbind(c(4,7,2), c(3,7,3)) #if I use line, sapply() below gives me a list #rbind(

[R] Debugging the "unexpected string constant" error in the zp2ssg2.R

2009-11-27 Thread Jason Rupert
I posted the full extent of this question to Nabble, but essentially I'm having difficulty determining the route cause of the "unexpected string constant" error in the zp2ssg2.R (which I'm trying to port over from Octave). Here is my current Nabble posting: http://n2.nabble.com/Need-help-solvi

[R] Questions about use of multinomial for discrimination.

2009-11-27 Thread Ricardo Antunes
Dear All, I am looking at discriminating among several individuals based on a few variable sets (I think some variables do not make sense unless they are entered together, so I "force" them into the models together, hence datasets). I have done so with linear discriminant analysis (LDA) using

[R] Orthogonal regression in R

2009-11-27 Thread Georg D. Blind
Dear R-Users! The most recent entries on orthogonal regression that I have found on this list date back to 2005. Is there, by now, any package available that allows for a multivariate orthogonal regression in R? Thanks and best regards. Georg georg.blind at gmx.net _

[R] my failing understanding ...

2009-11-27 Thread Kjetil Halvorsen
The following I do not understand, but then I did'nt really use S4 methods ... > showMethods(plot) Function: plot (package graphics) x="ANY" x="lmList.confint" x="merMCMC" (inherited from: x="ANY") > plot(x=moda0MCMC) Error in as.double(y) : cannot coerce type 'S4' to vector of type 'double

Re: [R] Graphic Device - View/get all graphics

2009-11-27 Thread Sarah Goslee
I'm not entirely sure I understand the problem, but one possible solution would be to use some device like postscript that saves the graphics to a file. (Or some other device appropriate for your intended use and OS.) You could still save the data to a list if you wanted - just make it part of wha

Re: [R] Export kde object as shapefile

2009-11-27 Thread T.D.Rudolph
On Thu, Nov 26, 2009 at 6:56 PM, David Winsemius wrote: http://old.nabble.com/Export-kde-object-as-shapefile-tp26532782p26545381.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org mailing list https://stat.ethz.ch/ma

Re: [R] Build of XML package failed

2009-11-27 Thread Duncan Temple Lang
Hi Luis. You can change the two lines PROBLEM buf WARN; to the one line warning(buf); That should compile. If not, please show us the compilation command for DocParse.c, i.e. all the arguments to the compiler, just above the error messages. D. Luis Tito de Morais wrote: > Hi list,

[R] Iteration idioms & laziness

2009-11-27 Thread Alexander Søndergaard
Hi all, I'm new to R. Having a functional background, I was wondering what's the idiomatic way to iterate. It seems that for loops are the default given there's no tail-call optimization. I'm curious to know whether there is a way to transform the following toy snippet into something that doesn't

Re: [R] using reshape to do ANOVA mixed models

2009-11-27 Thread Peter Dalgaard
Or Duek wrote: Hi, I just started with R and I found that there are many options to rearrange the data to do mixed models. I want to use the reshape function. I have 2 between subject variables and one within. I was able to change the data structure but still - the result of the aov functions are

[R] Build of XML package failed

2009-11-27 Thread Luis Tito de Morais
Hi list, It may be a FAQ, but I searched the web and Uni of Newcastle Maths and Stats and R mailing list archive on this issue but was unable to find a solution. I would appreciate any pointer to help me solving this. I am using R version 2.10.0 (2009-10-26) on linux mandriva 2010.0 I tried to i

Re: [R] Graphic Device - View/get all graphics

2009-11-27 Thread MarcioRibeiro
Hello, After some research, I didn't resolve my doubt yet... Maybe I didn't explain very well my question... So, I have a bootstrap simulation, but suppose just a simple part of my function: teste <-function(x){ rw_mean_app<-rnorm(x, mean = 10, sd = 2) rwy_mean_app<-rnorm(x, mean = 10, sd = 2) rw

Re: [R] Replace Values in Matrix

2009-11-27 Thread smu
matritz[is.na(matritz)] <- 0 On Fri, Nov 27, 2009 at 04:15:45PM -0200, Romildo Martins wrote: > Hello, > > how to replace the "NA" by number zero? > > > matrizt > [,1] [,2] [,3][,4] > [1,] 1.000NA NA NA > [2,] 0

[R] using reshape to do ANOVA mixed models

2009-11-27 Thread Or Duek
Hi, I just started with R and I found that there are many options to rearrange the data to do mixed models. I want to use the reshape function. I have 2 between subject variables and one within. I was able to change the data structure but still - the result of the aov functions are calculating ever

Re: [R] Modular inverses

2009-11-27 Thread William Dunlap
> -Original Message- > From: r-help-boun...@r-project.org > [mailto:r-help-boun...@r-project.org] On Behalf Of SJ Robson-Davis > Sent: Friday, November 27, 2009 8:52 AM > To: r-help@r-project.org > Subject: [R] Modular inverses > > I want to find the inverse of an integer k mod p (prime.)

Re: [R] Replace Values in Matrix

2009-11-27 Thread baptiste auguie
Hi, Try this, matrizt[is.na(matrizt)] <- 0 HTH, baptiste 2009/11/27 Romildo Martins : > Hello, > > how to replace the "NA" by number zero? > >> matrizt >                 [,1]           [,2]             [,3]        [,4] > [1,] 1.000            NA             NA         NA > [2,] 0.6717685 0

[R] Replace Values in Matrix

2009-11-27 Thread Romildo Martins
Hello, how to replace the "NA" by number zero? > matrizt [,1] [,2] [,3][,4] [1,] 1.000NA NA NA [2,] 0.6717685 0.1453253 NA NA [3,] 0.3971276 0.1493241 0.14532526 NA [4,] 0.1493241 0.1453

Re: [R] problem with "dynformula" from "plm" package [RE-POST]

2009-11-27 Thread Owen Powell
Hi David, Thank you for the response. I forgot to mention that I'd already tried what (I think) you propose (adding ".form" to the end of the "lag", "log" and "diff") and I still see the same results (posted below). Specifically, I still see no lags, logs or diffs in my model. Any other ideas?

Re: [R] problem with "dynformula" from "plm" package [RE-POST]

2009-11-27 Thread David Winsemius
On Nov 27, 2009, at 12:04 PM, Owen Powell wrote: Hi David, Thank you for the response. I forgot to mention that I'd already tried what (I think) you propose (adding ".form" to the end of the "lag", "log" and "diff") and I still see the same results (posted below). Specifically, I still see no

[R] Modular inverses

2009-11-27 Thread SJ Robson-Davis
I want to find the inverse of an integer k mod p (prime.) Is there a function that can do this for me? I know i could simply write (k^(p-2)) %% p, but i need to do this for large primes (above 100) and this gives the warning message: Warning message: probable complete loss of accuracy in modulus so

Re: [R] problem with "dynformula" from "plm" package [RE-POST]

2009-11-27 Thread David Winsemius
You might also note that the authors hint on the help page that one might want to use the formula() operation on the result of dynformula. Following that path would have gotten us to a more successful conclusion. > grun.fe <- plm(formula = formula(a),data=EmplUK) > grun.fe Model Formula: l

Re: [R] problem with "dynformula" from "plm" package [RE-POST]

2009-11-27 Thread Owen Powell
Great! That did the trick, thanks David. To summarize for the list, to get "dynformula" to work (example): R> a = dynformula(emp~wage+capital,log.form=list(capital=FALSE,TRUE),lag.form=list(emp=2,c(2,3)),diff.form=list(FALSE,capital=TRUE)) it might be necessary to use: R> plm(formula = formula

[R] R 2.10 Line Type Issue...

2009-11-27 Thread David B. Thompson, Ph.D., P.E., D.WRE, CFM
Morning folks (at least here on the western side of the U.S.)... This morning I constructed a contour plot of some bivariate distributions I'm working with. When I attempted to add a second contour to the plot using a dashed line (lty=2), R immediately went off to la-la land, requiring a f

Re: [R] Best way to preallocate numeric NA array?

2009-11-27 Thread William Dunlap
> -Original Message- > From: r-help-boun...@r-project.org > [mailto:r-help-boun...@r-project.org] On Behalf Of Rob Steele > Sent: Thursday, November 26, 2009 8:03 AM > To: r-h...@stat.math.ethz.ch > Subject: [R] Best way to preallocate numeric NA array? > > These are the ways that occur t

Re: [R] Simple Function doesn't work?

2009-11-27 Thread Colin Millar
Erm... Maybe the sequence bit wont work ... A bit hasty there And also the length check should be if (length(x) != length(y) & !(length(x) == 1 | length(y) == 1)) stop ("inputs not compatible") # or something I missed out the not! -Original Message- From: r-help-boun...@r-project.org

Re: [R] Simple Function doesn't work?

2009-11-27 Thread Colin Millar
Hi, You would also make your code more efficient and possible more readable by doing ReturnsGrid <- function(x, y, m) { x + (seq.int(m) - 1) * (y - x) / m } (xx <- ReturnsGrid(0, 9, 3)) #[1] 0 3 6 And if you want to supply vector x and y you could do something like (there are probably better

Re: [R] questions on the ff package

2009-11-27 Thread Jens Oehlschlägel
> I wonder how efficiently it is to do the following command on a frequent > basis. > nrow(matFF) <- nrow(matFF)+1 Obviously there is overhead (closing file, enlarging file, openeing file). I recommend you measure yourself whether this is acceptable for you. > no large file copying is needed ea

Re: [R] Save workspace with ff objects

2009-11-27 Thread Jens Oehlschlägel
> My script generates a mixture of normal and ff objects. > Very often I would like to save the workspace for > each parameter setting, so that I can get back to it later on. > Is there an easy way to do this, > instead of needing to save individual ff objects separately? With one save() you c

Re: [R] Force a variable substitution in expression

2009-11-27 Thread Saptarshi Guha
> The substitute() function can do it, but the bquote() function is the most > convenient: > >> foo <- function(s) { > +   bquote(expression({ > +    v <- .(s) > +    print(v) > +   })) > + } >> foo(10) > expression({ >   v <- 10 >   print(v) > }) > > Duncan Murdoch > Thanks, after sending this em

[R] changing titlessymbols in plot.spm

2009-11-27 Thread Assa Yeroslaviz
Hi all, I'm trying my best in changing the names of my experiments in the spectral map. My experiments are from six different time points, two as control and two treated experiments. I set have for each time point a different color . Now I want to change the symbols so that I have for the treated

Re: [R] Time Series Data

2009-11-27 Thread chris carleton
I'm not sure what you mean by 'multiplying by 10 and rounding'. I've tried to re-organize the data and then create a ts object that is analogous to the various sunspot datasets provided in the data package with R, but the matrix I generate by binning the data into centuries (which would make th

Re: [R] If condition using accessors

2009-11-27 Thread Vitória Magalhães Piai
Thanx. I'll just post in here the way I solved it in case someone happens to need something similar in the future. It may not be the most elegant way, but it worked! yes = "small" no = "large" data3$group = ifelse(data3$SIE < 50.64593,yes,no) Peter Ehlers wrote: Peter Ehlers wrote: Vitória

Re: [R] Simple Function doesn't work?

2009-11-27 Thread Alain Guillet
Hi, If you execute the following code it works but I wouldn't use grid if I were you as a vector as this name is already used by R (check help(grid)) and it explains why you have to define it in the function. ReturnsGrid = function(x,y,m){ grid <- numeric(m) for (i in 1:m){ grid[i] <- x + (

Re: [R] Force a variable substitution in expression

2009-11-27 Thread Duncan Murdoch
On 27/11/2009 11:04 AM, Saptarshi Guha wrote: Hello, I have a function that creates an expression object with some variables substituted e.g foo <- function(s){ expression({ v <- s print(v) }) } Thus foo returns an expression, however the expression has the symbol 's' contained with

Re: [R] Simple Function doesn't work?

2009-11-27 Thread baptiste auguie
Hi, The error message, Error in grid[i] <- x + (i - 1) * (y - x)/m : object of type 'closure' is not subsettable indicates that "grid" is actually known to R as a function (type grid to see its definition). You can define your own variable with the same name, but that needs to be done before t

Re: [R] Simple Function doesn't work?

2009-11-27 Thread Ista Zahn
Hi, You need to create the grid object before you can assign values to it. Try ReturnsGrid = function(x,y,m){ grid <- numeric() for (i in 1:m){ grid[i] <- x + (i-1)*(y-x)/m } grid } On Fri, Nov 27, 2009 at 11:00 AM, Anastasia wrote: > Hello, > > I am new to R program, therefore, I am sorry i

Re: [R] problem tick marker and text

2009-11-27 Thread Peter Ehlers
Duijvesteijn, Naomi wrote: Hi R-ers, I am struggling with my x-axis in a association plot. What I would like is to place the labels of the x-axis between the tick markers and normally the labels are printed at the place where the tick marker is placed. I don’t want to move the

[R] Force a variable substitution in expression

2009-11-27 Thread Saptarshi Guha
Hello, I have a function that creates an expression object with some variables substituted e.g foo <- function(s){ expression({ v <- s print(v) }) } Thus foo returns an expression, however the expression has the symbol 's' contained within it and thus returns an error when eval'd e.g

Re: [R] matching and extracting data

2009-11-27 Thread David Winsemius
On Nov 27, 2009, at 10:25 AM, ram basnet wrote: Dear all, I have querry on how to extract the data by matching between two data set where one has the same elements multiple times? For example, I have two matrix X and Y. X [,1][,2] [,3] 1 A 5 P 2 B 6

[R] Simple Function doesn't work?

2009-11-27 Thread Anastasia
Hello, I am new to R program, therefore, I am sorry if this is a really stupid question. I wrote a simple function and for some reason it doesn't work ReturnsGrid = function(x,y,m){ for (i in 1:m){ grid[i] <- x + (i-1)*(y-x)/m } grid } xx=ReturnsGrid(0,9,3) Thanks a lot! [[alternati

Re: [R] matching and extracting data

2009-11-27 Thread Chuck Cleland
On 11/27/2009 10:25 AM, ram basnet wrote: > Dear all, > I have querry on how to extract the data by matching between two data set > where one has the same elements multiple times? > > For example, I have two matrix X and Y. >> X >[,1][,2] [,3] > 1 A 5 P > 2 B 6

[R] matching and extracting data

2009-11-27 Thread ram basnet
Dear all, I have querry on how to extract the data by matching between two data set where one has the same elements multiple times?   For example, I have two matrix X and Y. > X    [,1]    [,2]   [,3] 1  A     5 P 2  B  6     P 3  C 7     P 4  D 5     Q 5  E   

Re: [R] problem with "dynformula" from "plm" package [RE-POST]

2009-11-27 Thread David Winsemius
On Nov 27, 2009, at 10:25 AM, Owen Powell wrote: Hello list, I'm following the paper (http://www.jstatsoft.org/v27/i02/paper) on how to use "plm" to run panel regressions, and am having trouble with what I believe should be something very basic. When I run the command (p.9 in the paper): R>

[R] problem with "dynformula" from "plm" package [RE-POST]

2009-11-27 Thread Owen Powell
Hello list, I'm following the paper (http://www.jstatsoft.org/v27/i02/paper) on how to use "plm" to run panel regressions, and am having trouble with what I believe should be something very basic. When I run the command (p.9 in the paper): R> dynformula(emp~wage+capital,log=list(capital=FALSE,TR

Re: [R] Time Series Data

2009-11-27 Thread David Winsemius
On Nov 27, 2009, at 9:55 AM, chris carleton wrote: Hi All, I'm trying to analyze some time series data and I have run into difficulty. I have decadal sun spot data and I want to separate the very regular periodic function from the trend and noise. I looked into using stl(), but the freq

[R] problem with "dynformula" from "plm" package

2009-11-27 Thread Owen Powell
Hello list, I'm following the paper (http://www.jstatsoft.org/v27/i02/paper) on how to use "plm" to run panel regressions, and am having trouble with what I believe should be something very basic. When I run the command (p.9 in the paper): R> dynformula(emp~wage+capital,log=list(capital=FALSE,TR

Re: [R] post

2009-11-27 Thread Peter Dalgaard
Anastasia wrote: > Dear all, > > I have the following problem which I cannot solve: > > data <- numeric (100) > for (i in 1:100){ > p <- runif(1,min=0,max=1) > data[i] <- rnorm(1,mean=2,sd=1) > } Er what is "p" doing in there. As far as I can see, you might as well do data <- rnorm(200, 2,

[R] problem tick marker and text

2009-11-27 Thread Duijvesteijn, Naomi
Hi R-ers, I am struggling with my x-axis in a association plot. What I would like is to place the labels of the x-axis between the tick markers and normally the labels are printed at the place where the tick marker is placed. I don’t want to move the tick marker (it gives the sw

[R] Time Series Data

2009-11-27 Thread chris carleton
Hi All, I'm trying to analyze some time series data and I have run into difficulty. I have decadal sun spot data and I want to separate the very regular periodic function from the trend and noise. I looked into using stl(), but the frequency of the time series data must be greater than 1 for s

[R] post

2009-11-27 Thread Anastasia
Dear all, I have the following problem which I cannot solve: data <- numeric (100) for (i in 1:100){ p <- runif(1,min=0,max=1) data[i] <- rnorm(1,mean=2,sd=1) } ke <- density(data,bw="sj",n=61) How can I now find the value of this density function for example ke(0), ke(2) etc. Maybe this is

Re: [R] If condition using accessors

2009-11-27 Thread David Winsemius
On Nov 27, 2009, at 5:34 AM, Vitória Magalhães Piai wrote: Hi, I'm quite new using R and have got no one to help me get through it. Hopefully someone can help me with one problem I've been struggling with for the last hours!! (Sorry if I'm using the wrong terminology as well!) I have a d

Re: [R] layers in xYplot of Hmisc

2009-11-27 Thread David Winsemius
I took the OP's variables, reversed the "direction" of the second y, bundled them up in a dataframe and then assigned the grp. So starting with the OP's code: x<-seq(1,10,1) y<-seq(1,10,1) ci<-y*.10 ciupper<-y+ci cilower<-y-ci x2<-seq(1,5,.5) y2<-seq(5, 1, -.5) #reverse second variable ci2<-

[R] how to analyze this design using lmer

2009-11-27 Thread Adaikalavan Ramasamy
Dear all, A friend of mine requested me to analyze some data she has generated. I am hoping for some advice on best way of properly analyzing the data as I have never worked with such complicated or nested designs. Here is the setup. She has taken material from 5 animals and each material is

Re: [R] Barplot with confidence intervals

2009-11-27 Thread Marc Schwartz
On Nov 26, 2009, at 3:23 PM, Jonhnny Weslley wrote: I have a file in the following format: Scenario1 Scenario1CIL Scenario1CIU Scenario2 Scenario2CIL Scenario2CIU 60 57 62 45 48 50 110 101 111 51 50 52 120 117 122 64 62 66 192 190 194 79 75 79 where: First column = Scenario1 mean value Se

Re: [R] If condition using accessors

2009-11-27 Thread Peter Ehlers
Peter Ehlers wrote: Vitória Magalhães Piai wrote: Hi, I'm quite new using R and have got no one to help me get through it. Hopefully someone can help me with one problem I've been struggling with for the last hours!! (Sorry if I'm using the wrong terminology as well!) I have a data matrix

Re: [R] If condition using accessors

2009-11-27 Thread Peter Ehlers
Vitória Magalhães Piai wrote: Hi, I'm quite new using R and have got no one to help me get through it. Hopefully someone can help me with one problem I've been struggling with for the last hours!! (Sorry if I'm using the wrong terminology as well!) I have a data matrix in which SIE is one of

[R] plotting two surfaces simultaneously in a single panel

2009-11-27 Thread Umesh Srinivasan
Hi, I have recently begun using the lattice package, and have been using the wireframe command to visualise matrices which are model outputs. I have been trying to plot two surfaces (from two matrices) simultaneously in one panel, to visualise intersections etc., but neither my attempts or trawlin

Re: [R] barchart() {Lattice} help.

2009-11-27 Thread Deepayan Sarkar
On Fri, Nov 27, 2009 at 9:21 AM, Peng Cai wrote: > @ Peter: I got it, thanks a lot for all your help! And yes, as you said the > "title" option in auto.key is redundant. > > > @ All, Hi: I need to add percentage sign to y-axis labels (like 0%, 20%, > ..., 100%). How can I get it. I'm using barchar

Re: [R] Adding text in the panels for Trellis plot ...

2009-11-27 Thread Deepayan Sarkar
On Thu, Nov 26, 2009 at 6:22 AM, ychu066 wrote: > > i was trying to do a for loop for plotting the histograms , but it doesnt > work properly > >> library(lattice) >> columns <- 8:153 >> plots <- vector("list", length(columns)) >> j <- 0 >> for (i in columns) > + { > +   plots[[ j <- j+1 ]] <- his

Re: [R] Concave hull

2009-11-27 Thread baptiste auguie
Thanks for the interesting reference to alphahull. It might be a good starting point for placing e.g. a legend in a plot (I think the usual techniques for this (gregmisc?) are a bit more brute-force.) baptiste 2009/11/27 Kjetil Halvorsen : > There is a package on CRAN implementing such an idea:

Re: [R] order of panels in xyplots

2009-11-27 Thread Deepayan Sarkar
On Thu, Nov 26, 2009 at 12:09 AM, Titus Malsburg wrote: > The documentation of xyplot could be improved here.  It says: > >  "If 'index.cond' is a list, it has to be as long as the number of > conditioning >   variables, and the 'i'-th component has to be a valid indexing vector for > the >   in

Re: [R] lattice --- different properties of lines corresponding to type=c("l", "a") respectively

2009-11-27 Thread Deepayan Sarkar
On Fri, Nov 27, 2009 at 1:13 AM, Kjetil Halvorsen wrote: > I think the subject says it all. I want to make a simple lattice plot, > using xyplot with the > argument     type=c("l","a"). > The problem then is that in the resulting plot it is > difficult/impossible to see which plot corresponds to t

Re: [R] Arguments

2009-11-27 Thread Gabor Grothendieck
See ?commandArgs Also the getopt and optparse CRAN packages. On Fri, Nov 27, 2009 at 4:20 AM, yonosoyelmejor wrote: > > hello,I would like to ask a question,Is there any way to pass arguments to a > script?I have this code: > > Invernadero<-read.table(file.choose(),header=T,sep=",") > > attach(In

[R] Arguments

2009-11-27 Thread yonosoyelmejor
hello,I would like to ask a question,Is there any way to pass arguments to a script?I have this code: Invernadero<-read.table(file.choose(),header=T,sep=",") attach(Invernadero) names(Invernadero) Invernadero<-ts(Invernadero<-argument) //Here introduce the argument plot(Hora,Invernadero,main=ar

Re: [R] R help with princomp and pam clustering

2009-11-27 Thread Tyler82
Hi JorisMeys and thanks! JorisMeys wrote: > > On Thu, Nov 26, 2009 at 1:04 AM, Tyler82 > wrote: >> >> Hi all! >> I am working with R package cluster and I have a little problem: >> let's say I have two datasets...first one ("A") is divided into 4 >> clusters >> by means of Pam algorythm. >> Le

[R] If condition using accessors

2009-11-27 Thread Vitória Magalhães Piai
Hi, I'm quite new using R and have got no one to help me get through it. Hopefully someone can help me with one problem I've been struggling with for the last hours!! (Sorry if I'm using the wrong terminology as well!) I have a data matrix in which SIE is one of my variables. What I need to

Re: [R] Learning R - View datasets

2009-11-27 Thread ojal john owino
Hi Brock, Have you tried View() ? Regards. On Fri, Nov 27, 2009 at 7:46 AM, Brock Tibert wrote: > Hi All, > > I am making a serious effort to try to learn R, but one hurdle I am facing > is that I need to "see" the data as I walk through the examples in the > packages. For instance, many exam

[R] generating a matrix after a for loop..

2009-11-27 Thread John Seppänen
Hi all, I have to ask this and I know that the reason is that I am a newbie with R programming. So apologize if it is too obvious but I didn't find an answer after googling and reading "An introduction to R". So i have return data from 30 instruments and I am fitting a mixture of normal distributi

Re: [R] Learning R - View datasets

2009-11-27 Thread Gabor Grothendieck
Try this: # each of these three show entire data set wines dput(wines) View(wines) # get help ?wines # various info on data set head(wines) tail(wines) summary(wines) str(wines) class(wines) dim(wines) # plotting plot(wines) # for a better plot see the example at the bottom of ?wines On Thu, N

Re: [R] layers in xYplot of Hmisc

2009-11-27 Thread butter
do you mind me asking what code you used to create that data frame and name the groups 1 and 2? David Winsemius wrote: > > > On Nov 26, 2009, at 10:26 PM, Joe King wrote: > >> In the "filled bands" part of xYplot of the Hmisc package, is there >> a way to >> have multiple bands with multiple

Re: [R] Learning R - View datasets

2009-11-27 Thread Liviu Andronic
Hello On Fri, Nov 27, 2009 at 4:46 AM, Brock Tibert wrote: > In short, I know that data() lists all of the available datasets, > data("wines") will load the dataset wines, but how can I look at the raw data? > See this [1]. [1] http://www.mail-archive.com/r-help@r-project.org/msg66111.html Livi

Re: [R] GAM plots

2009-11-27 Thread Simon Wood
Are you using mgcv:::gam? To get plot data suitable for making plots of smooth effects, you probably need to use `predict.gam' to evaluate the smooth curves (and standard errors) at a nice regular set of points for plotting. Also don't forget that the residuals shown on plot.gam are the `partial

[R] VAR forecasts and out-of-sample prediction

2009-11-27 Thread manta
Dear users, I am struggling with this issue. I want to estimate a VAR(1) for three variables, say beta1 beta2 beta3, using monthly observations from January 1984 to September 2009. In-sample period January 1984 to December 2003, out-of-sample January 2004 to September 2009. This is what I have don

Re: [R] How to compute Rolling analysis of Standard Deviation using ZOO package?

2009-11-27 Thread Gabor Grothendieck
We will need a minimal reproducible example (as per last line on every message to r-help) to answer as trying it with made up data seems to work for me. If c is very long try cutting it down to the smallest you can get it to and still produce the error, e.g. cc <- c[1:10] rollapply(cc, ...) and

[R] Question about S4

2009-11-27 Thread Hun S. Tesatte
I'm curious about why no one has answered my question below. I can't imagine it would be because no one knows how to answer, it must be something basic I am ignorant about. But I have never seen such a pattern, it seems strange to me that a class with an empty definition is automatically virtua

[R] multicore: defunct R processes left

2009-11-27 Thread Ramon Diaz-Uriarte
Dear All, At least in three different GNU/Linux systems, the parallel function from the multicore package leaves defunct (zombie) R processes. For instance library(multicore) parallel(1:5) collect() After collect, the child process (with pid as given by collect) is left defunct. (If we run th

[R] Comparing membership of clusters

2009-11-27 Thread Paul Lemmens
Hello, I'm taking several physiological measurements on participants (e.g., skin conductivity, heart rate, etc.). I know that those participants belong to one of three groups (from another measurement), and I'm looking to find the physiological measurement that best describes group membership. The

Re: [R] Best way to preallocate numeric NA array?

2009-11-27 Thread Petr PIKAL
Hi There is one issue which I encountered recently with this type of behaviour, mat<-matrix(NA,5,4) fix(mat) put some number in any cell and close fix mat col1 col2 col3 col4 [1,] NA NA NA NA [2,] NA NA NA NA [3,] NA NA NA NA [4,] NA NA NA NA [5,] NA N

Re: [R] Learning R - View datasets

2009-11-27 Thread Paul Hiemstra
Brock Tibert wrote: Hi All, I am making a serious effort to try to learn R, but one hurdle I am facing is that I need to "see" the data as I walk through the examples in the packages. For instance, many examples on the web start by a command like data("wines"). How can I actually view what

Re: [R] How to compute Rolling analysis of Standard Deviation using ZOO package?

2009-11-27 Thread Karl Ove Hufthammer
On Fri, 27 Nov 2009 15:19:11 +0800 Saji Ren wrote: > I want to get a rolling estimation of the stdev of my data. There is a 'runsd' in the 'caTools' package which does exactly this. -- Karl Ove Hufthammer __ R-help@r-project.org mailing list https://

Re: [R] Learning R - View datasets

2009-11-27 Thread Wincent
There are different ways to inspect the conent of a data frame. For example, >View(CO2) 2009/11/27 Brock Tibert : > Hi All, > > I am making a serious effort to try to learn R, but one hurdle I am facing is > that I need to "see" the data as I walk through the examples in the packages. >  For i

Re: [R] Learning R - View datasets

2009-11-27 Thread guohao.huang
, November 27, 2009 4:00 PM To: Cc: "Brock Tibert" ; Subject: Re: [R] Learning R - View datasets guohao.hu...@gmail.com wrote: > Please check the following pdf file. > http://tw.nextmedia.com/applenews/article/art_id/32119622/IssueID/20091127 > Besides a beautiful blond gir

Re: [R] Overlapping x - axis lables

2009-11-27 Thread Patrick Connolly
On Thu, 26-Nov-2009 at 11:03PM -0800, Julia Cains wrote: |> Dear Sir, |> |> Thanks again for your help. Sir, actually I am working on a data containing 84 months and the corresponding savings (core) balances and non-core balances and I am plotting a stacked bar. |> |> The original code I am us

Re: [R] Learning R - View datasets

2009-11-27 Thread Mario Valle
guohao.hu...@gmail.com wrote: > Please check the following pdf file. > http://tw.nextmedia.com/applenews/article/art_id/32119622/IssueID/20091127 > Besides a beautiful blond girl, I do not recognize any pdf here... Could you provide a direct link? Thanks! mario &