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 1. First

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 using

Re: [R] Learning R - View datasets

2009-11-27 Thread guohao.huang
/20091127 Besides a beautiful blond girl, I do not recognize any pdf here... Could you provide a direct link? Thanks! mario 1. First install.packages(Flury) 2. library(Flury) 3. data(wines) 'wines’ is a data frame with 26 observations, one factor denoting the country of origin and 15

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 btibe...@yahoo.com: 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

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 saji@gmail.com 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

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 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

[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.

[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

[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

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] 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

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

Re: [R] Learning R - View datasets

2009-11-27 Thread Liviu Andronic
Hello On Fri, Nov 27, 2009 at 4:46 AM, Brock Tibert btibe...@yahoo.com 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]

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 lines? or

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,

[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

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 btibe...@yahoo.com 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

[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] 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 procaccianti.clau...@gmail.com 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

[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

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 yonosoyelme...@hotmail.com wrote: hello,I would like to ask a question,Is there any way to pass arguments to a script?I have this code:

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 kjetilbrinchmannhalvor...@gmail.com 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

Re: [R] order of panels in xyplots

2009-11-27 Thread Deepayan Sarkar
On Thu, Nov 26, 2009 at 12:09 AM, Titus Malsburg malsb...@gmail.com 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

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 kjetilbrinchmannhalvor...@gmail.com: There is a package

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 ychu...@aucklanduni.ac.nz 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 ]] -

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

2009-11-27 Thread Deepayan Sarkar
On Fri, Nov 27, 2009 at 9:21 AM, Peng Cai pengcaimaill...@gmail.com 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.

[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

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

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] 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

[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

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-y2*.10

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

[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 an

[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

[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

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, 1) ke -

[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

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

[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

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] 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] 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 P

[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!

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] 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] 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 nast...@gmail.com wrote: Hello, I am new to R program, therefore,

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

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 the

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

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] 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] 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

[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

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 email, i used

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 can

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 each

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] 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] 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 to me.

[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

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 =

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:

[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

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 lags,

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? ~Owen

[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

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 romildo.mart...@gmail.com: Hello, how to replace the NA by number zero? matrizt                 [,1]           [,2]             [,3]        [,4] [1,] 1.000            NA             NA         NA

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.) Is

[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

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.6717685

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)

[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

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

[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

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,

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 dwinsem...@comcast.net wrote: The situation that I see (after looking at the documentation for adehabit) is that for some as yet unarticulated reason, you have decided that the methods used in prior publications in your domain are not the best and

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

[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'

[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] 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] 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:

[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

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

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] How to combine matrix with sapply function?

2009-11-27 Thread Ben Bolker
Peng Yu pengyu.ut at 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] 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: not a generic function 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

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: not a generic function sessionInfo() R version 2.10.0 Patched (2009-10-29 r50258) x86_64-apple-darwin9.8.0 locale: [1]

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 ggrothendi...@gmail.com 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

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 (solo

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 jasonkrup...@yahoo.com To: R-help@r-project.org Sent:

[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