[R] quantmod's addTA plotting functions

2011-05-04 Thread Russ Abbott
Hi, I'm having trouble with quantmod's addTA plotting functions. They seem to work fine when run from the command line. But when run inside a function, only the last one run is visible. Here's an example. test.addTA <- function(from = "2010-06-01") { getSymbols("^GSPC", from = from) GS

Re: [R] best subset regression in R

2011-05-04 Thread Bert Gunter
On Wed, May 4, 2011 at 9:47 AM, FMH wrote: > Dear All, > > Could someone please give some advice the way to do linear modelling via best > subset regression in R?... Yes. Don't do it. -- Bert (Very Brief Explanation: Best subset regression was a questionable approach to parsimonious modeling

Re: [R] bivariate linear interpolation

2011-05-04 Thread Andrew Robinson
The one that gives results that you trust and uses algorithms that you understand! Cheers Andrew On Wed, May 04, 2011 at 12:52:03PM +, Halldór Björnsson wrote: > Hi, > > I have three matrices (X,Y,P) with the same dimension. The X,Y grid is > regular and I want to > perform linear interpola

Re: [R] SAPPLY function XXXX

2011-05-04 Thread Petr PIKAL
Hi r-help-boun...@r-project.org napsal dne 04.05.2011 22:26:59: > Erik Iverson > Odeslal: r-help-boun...@r-project.org > > 04.05.2011 22:26 > > Komu > > Dan Abner > > > > Ultimately, I would like for this to be 1 conponent in a larger function > > that will produce PROC CONTENTS style ou

Re: [R] Uniform Gaussian Kernel

2011-05-04 Thread Andrew Robinson
I can't honestly say that I grasp what you're trying to do, but that said, I wonder if the curve() function will help you? Cheers Andrew On Wed, May 04, 2011 at 07:30:20AM -0700, blutack wrote: > I have a vector with lots of different numbers. I need to make a graph > showing the Uniform Distri

Re: [R] split character vector by multiple keywords simultaneously

2011-05-04 Thread Andrew Robinson
A hack would be to use gsub() to prepend e.g. XXX to the keywords that you want, perform a strsplit() to break the lines into component strings, and then substr() to extract the pieces that you want from those strings. Cheers Andrew On Wed, May 04, 2011 at 04:08:40PM -0700, sunny wrote: > Hi. I

Re: [R] two-way group mean prediction in survreg with three factors

2011-05-04 Thread Andrew Robinson
I hope not! Facetiousness aside, the model that you have fit contains C, and, indeed, an interaction between A and C. So, the effect of A upon the response variable depends on the level of C. The summary you want must marginalize C somehow, probably by a weighted or unweighted average across its

Re: [R] fGarch

2011-05-04 Thread Andrew Robinson
Hi Paul, I suggest that you should send us commented, minimal, self-contained, reproducible code. That means, in essence, developing the simplest possible representation of your problem. In the process of developing the simplest possible representation, you may learn more about the problem. May

Re: [R] Obtaining the name of an object XXXX

2011-05-04 Thread Rolf Turner
On 05/05/11 11:16, Dan Abner wrote: Hello everyone, How does one write a function to return the name of an input object (that is assumed to be a data frame) as a character string? I tired using the get(), but this does not work as I had hoped. For example: myfn<-function(x){ output<-data.frame

Re: [R] print elements' placings from vector

2011-05-04 Thread derek
Yes, thank you. -- View this message in context: http://r.789695.n4.nabble.com/print-elements-placings-from-vector-tp3496645p3497080.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org mailing list https://stat.ethz.ch

[R] split character vector by multiple keywords simultaneously

2011-05-04 Thread sunny
Hi. I have a character vector that looks like this: > temp <- c("Company name: The first company General Manager: John Doe I > Managers: John Doe II, John Doe III","Company name: The second company > General Manager: Jane Doe I","Company name: The third company Managers: > Jane Doe II, Jane Doe II

[R] Obtaining the name of an object XXXX

2011-05-04 Thread Dan Abner
Hello everyone, How does one write a function to return the name of an input object (that is assumed to be a data frame) as a character string? I tired using the get(), but this does not work as I had hoped. For example: myfn<-function(x){ output<-data.frame(Attribute="Data Set Name",Value=as.ch

Re: [R] Panels order in lattice graphs

2011-05-04 Thread Andrew Robinson
Hi Cristina, you can probably hack your own solution using the index.cond argument. Cheers Andrew On Wed, May 04, 2011 at 04:50:53PM +0100, Cristina Silva wrote: > Hi all, > > In lattice graphs, panels are drawn from left to right and bottom to > top. The flag "as.table=TRUE" changes to left

Re: [R] Problems saving ff objects

2011-05-04 Thread Andrew Robinson
I wonder if this question should be directed to the package maintainer? Best wishes, Andrew On Wed, May 04, 2011 at 02:31:51PM +0100, Jannis wrote: > Just did some more testing.May the problem be due to the fact that I am > using a windows machine? I just ran the same code on a Linux machin

Re: [R] help with the maxBHHH routine

2011-05-04 Thread Andrew Robinson
Hi Rohit, actually, the request for simple reproducible code means that you have to find the simplest possible representation of the problem. What happens if you simplify the observation level gradient and the likelihood function? Eg to trivial examples? If you still get the error, then simpl

Re: [R] combine lattice plot and standard R plot

2011-05-04 Thread Dennis Murphy
Hi: If it's doable, you'll probably need the gridBase package. Fortunately, it has a nice vignette to get you started, which tells you at the end that there are limitations in compatibility between base and grid graphics (lattice is built on the latter). HTH, Dennis 2011/5/4 Lucia Cañas : > Dear

Re: [R] how to replace all variable values?

2011-05-04 Thread Dennis Murphy
Hi: > x1 <- 1:5 > x2 <- 6:10 > x3 <- c(1, 2, 5) > x4 <- 3 > x5 <- 4 > x2[x1 %in% x3] [1] 6 7 10 > x2[x1 %in% x4] [1] 8 > x2[x1 %in% x5] [1] 9 HTH, Dennis On Wed, May 4, 2011 at 2:20 PM, blue100 wrote: > dear all, > > Im a complete R newby with the following question. > > I have a dataset whe

Re: [R] print elements' placings from vector

2011-05-04 Thread Dennis Murphy
Hi: Is this what you're after? > z=c(1,3,0.5,6,8,10,2,2,3,4,7,3) > which(z > 2) [1] 2 4 5 6 9 10 11 12 HTH, Dennis On Wed, May 4, 2011 at 1:02 PM, derek wrote: > Dear R, > > Here is a code: > z=c(1,3,0.5,6,8,10,2,2,3,4,7,3) > z[z>2] > > I dont want print the elements, but theirs placings

Re: [R] best subset regression in R

2011-05-04 Thread Frank Harrell
Beware - this approach is a statistical train wreck. Been there, done that. If all you want is "an" answer it will save you a lot of time thinking, however. Frank FMH-4 wrote: > > Dear All, >   > Could someone please give some advice the way to do linear modelling via > best subset regression

Re: [R] Box-Cox transformation in R

2011-05-04 Thread Dennis Murphy
Hi: Start here: library(sos)# Install first if necessary findFn('Box-Cox') This search finds 131 matches; the basic Box-Cox transformations for regression are found in the MASS and car packages. For other situations, consult the packages and functions identified from the sos search. HTH, De

Re: [R] combine lattice plot and standard R plot

2011-05-04 Thread Abhijit Dasgupta
Hi, All of the components in the grid.arrange statement are from ggplot2 and lattice, which are both in turn based on the grid package. What grid.arrange is able to do is use the grid framework to arrange the individual plots on a page. The base graphics are not based on grid, and so won't wor

[R] how to replace all variable values?

2011-05-04 Thread blue100
dear all, Im a complete R newby with the following question. I have a dataset where my variable values are incorrectly numbered. it has to be something like this Where x is variable name, y=actual value which must become corresponding z-value x3 x4 x5

[R] fGarch

2011-05-04 Thread Paul Ossenbruggen
Hi, I am attempting to fit a ARMA/GARCH regression model without success. ### ARIMA-GARCH model with regressor ### ### Time series data: A multivariate data set. cov.ts.dq = cov.ts[1:4,"dq1"][!is.na(cov.ts[,"dq1"])] cov.ts.day = ts.intersect(dq = diff(q.ts), day = lag(q.t

[R] print elements' placings from vector

2011-05-04 Thread derek
Dear R, Here is a code: z=c(1,3,0.5,6,8,10,2,2,3,4,7,3) z[z>2] I dont want print the elements, but theirs placings in vector. -- View this message in context: http://r.789695.n4.nabble.com/print-elements-placings-from-vector-tp3496645p3496645.html Sent from the R help mailing list archive at

Re: [R] ANOVA 1 too few degrees of freedom

2011-05-04 Thread Rovinpiper
And I responded as follows: Hi, Thanks for your advice. I tried using table() to check for missing data. Here are the results: > table(Combined.Plot) Combined.Plot 60m A1 B1 B3 B4 C5 C9 D2 D9 F60m F8 Q7 34 34 3434 343434 3434 34 34 34 >

Re: [R] merging multiple columns from two dataframes

2011-05-04 Thread Scott Chamberlain
merge(dat, pheno, by.x = "sid", by.y = "sen") Scott On Wednesday, May 4, 2011 at 12:53 PM, 1Rnwb wrote: > Hello, > > I have data in a dataframe with 139104 rows which is multiple of 96x1449. i > have a phenotype file which contains the phenotype information for the 96 > samples. the snp name i

Re: [R] ANOVA 1 too few degrees of freedom

2011-05-04 Thread Rovinpiper
This response went to my email: Without your data it's hard to say, but one possibility is that your plots are nested within treatments instead of crossed, or that you have something rather more cunning going on involving the Days. For example if you had 8 days for six of your plots and another 8

Re: [R] Bigining with a Program of SVR

2011-05-04 Thread ypriverol
How can I apply feature selection with caret and support vector regression. -- View this message in context: http://r.789695.n4.nabble.com/Bigining-with-a-Program-of-SVR-tp3484476p3496685.html Sent from the R help mailing list archive at Nabble.com. ___

Re: [R] SAPPLY function XXXX

2011-05-04 Thread Dan Abner
Perfect Erik! Thank you! On Wed, May 4, 2011 at 4:22 PM, Erik Iverson wrote: > Dan, > > > > I am attempting to write a function to count the number of non-missing >> values of each column in a data frame using the sapply function. I have >> the >> following code which is receiving the error me

Re: [R] help with the maxBHHH routine

2011-05-04 Thread Rohit Pandey
Hi Andrew, Ravi and Arne, Thank you so much for your prompt replies. I see that all of you mention the need for simple, reproducible code. I had thought of doing this, but the functions I was using for the observation level gradient and likelihood function were very long. I will paste them below h

Re: [R] what happens when I store linear models in an array?

2011-05-04 Thread Andrew D. Steen
Thanks all, this is very helpful. --Andrew Steen > -Original Message- > From: David Winsemius [mailto:dwinsem...@comcast.net] > Sent: Wednesday, May 04, 2011 4:35 PM > To: Andrew D. Steen > Cc: r-help@r-project.org > Subject: Re: [R] what happens when I store linear models in an array? >

Re: [R] Regexp question

2011-05-04 Thread johannes rara
Thank you all! 2011/5/4 David Wolfskill : > On Wed, May 04, 2011 at 10:41:36PM +0300, johannes rara wrote: >> I have a string like this >> >> st <- "SELECT COUNT(empid), COUNT(mgrid), COUNT(empname), >> COUNT(salary), FROM Employees" >> >> How can I remove the last comma before the FROM statement?

Re: [R] Box-Cox transformation in R

2011-05-04 Thread Patrick Breheny
On 05/04/2011 12:53 PM, FMH wrote: Hi, Could any one please help how I can transform data based on Box-Cox Transformations in R. Any helps will be much appreciated. thanks, Kagba [[alternative HTML version deleted]] See the boxcox function in the MASS package.

Re: [R] Regexp question

2011-05-04 Thread David Wolfskill
On Wed, May 04, 2011 at 10:41:36PM +0300, johannes rara wrote: > I have a string like this > > st <- "SELECT COUNT(empid), COUNT(mgrid), COUNT(empname), > COUNT(salary), FROM Employees" > > How can I remove the last comma before the FROM statement? This doesn't use a regex, per se, but: > st <-

Re: [R] Str info. Thanks for helping

2011-05-04 Thread Ivan Calandra
It looks from str(SA) that Response IPS1 is a data.frame of class "anova", which probably cannot be coerced to vector. Maybe you can use unlist() instead of as.vector() Or something like SA[["Response IPS1"]]["as.factor(WSD)",] ## to select the first row only, even maybe with unlist() Without

Re: [R] best subset regression in R

2011-05-04 Thread Jeremy Miles
On 4 May 2011 09:47, FMH wrote: > Dear All, > > Could someone please give some advice the way to do linear modelling via best > subset regression in R? I'd really appreciate for your kindness. > Google is your friend here: http://www.google.com/search?q=best+subsets+regression+R , and sends me

Re: [R] SAPPLY function XXXX

2011-05-04 Thread Erik Iverson
Ultimately, I would like for this to be 1 conponent in a larger function that will produce PROC CONTENTS style output. Something like... data1.contents<-data.frame(Variable=names(data1), Class=sapply(data1,class), n.valid=sapply(data1,sum(!is.na)), n.miss=sapply(data1,sum(is.na))) data1.cont

Re: [R] Regexp question

2011-05-04 Thread Jerome Asselin
On Wed, 2011-05-04 at 22:41 +0300, johannes rara wrote: > I have a string like this > > st <- "SELECT COUNT(empid), COUNT(mgrid), COUNT(empname), > COUNT(salary), FROM Employees" > > How can I remove the last comma before the FROM statement? gsub(",[^,]*FROM ", " FROM ", st) HTH, Jerome __

Re: [R] SAPPLY function XXXX

2011-05-04 Thread Erik Iverson
Dan, I am attempting to write a function to count the number of non-missing values of each column in a data frame using the sapply function. I have the following code which is receiving the error message below. n.valid<-sapply(data1,sum(!is.na)) Error in !is.na : invalid argument type Tha

Re: [R] join tables in R

2011-05-04 Thread Steve Lianoglou
Hi, On Wed, May 4, 2011 at 10:32 AM, Roccato Alfredo (UniCredit) wrote: > I'd to match-merge 2 tables in such a manner that I keep all the rows in > table 1, but not the rows that are in both table 1 and 2. > Thank you for your help, > Alfredo > >> master <- data.frame(ID=2001:2011) >> train   <

Re: [R] Regexp question

2011-05-04 Thread Henrique Dallazuanna
Try this: gsub(",\\s+FROM", " FROM", st) On Wed, May 4, 2011 at 4:41 PM, johannes rara wrote: > I have a string like this > > st <- "SELECT COUNT(empid), COUNT(mgrid), COUNT(empname), > COUNT(salary), FROM Employees" > > How can I remove the last comma before the FROM statement? > > -J > > _

Re: [R] join tables in R

2011-05-04 Thread Petr Savicky
On Wed, May 04, 2011 at 04:32:14PM +0200, Roccato Alfredo (UniCredit) wrote: > I'd to match-merge 2 tables in such a manner that I keep all the rows in > table 1, but not the rows that are in both table 1 and 2. > Thank you for your help, > Alfredo > > > master <- data.frame(ID=2001:2011) > > tra

Re: [R] Saving Values in a Vector from a For Loop

2011-05-04 Thread Steve Lianoglou
Hi, On Wed, May 4, 2011 at 10:19 AM, blutack wrote: > Hi, > I have a created a function, but now I need to call it about a hundred times > and store the results as a vector. > I think doing a for loop would work, but I cant work out how to save the > values generated from the function as a vector

Re: [R] join tables in R

2011-05-04 Thread Steve_Friedman
Look at the merge command ?merge Steve Friedman Ph. D. Ecologist / Spatial Statistical Analyst Everglades and Dry Tortugas National Park 950 N Krome Ave (3rd Floor) Homestead, Florida 33034 steve_fried...@nps.gov Office (305) 224 - 4282 Fax (305) 224 - 4147

[R] SAPPLY function XXXX

2011-05-04 Thread Dan Abner
Hello everyone, I am attempting to write a function to count the number of non-missing values of each column in a data frame using the sapply function. I have the following code which is receiving the error message below. > n.valid<-sapply(data1,sum(!is.na)) Error in !is.na : invalid argument ty

[R] Box-Cox transformation in R

2011-05-04 Thread FMH
Hi,   Could any one please help how I can transform data based on Box-Cox Transformations in R.   Any helps will be much appreciated.   thanks, Kagba [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.c

[R] best subset regression in R

2011-05-04 Thread FMH
Dear All,   Could someone please give some advice the way to do linear modelling via best subset regression in R? I'd really appreciate for your kindness.   Thanks, Kagba [[alternative HTML version deleted]] __ R-help@r-project.org mailing list

[R] Outlier removal by Principal Component Analysis : error message

2011-05-04 Thread boule
Hi, I am currently analysis Raman spectroscopic data with the hyperSpec package. I consulted the documentation on this package and I found an example work-flow dedicated to Raman spectroscopy (see the address : http://hyperspec.r-forge.r-project.org/chondro.pdf) I am currently trying to remove ou

[R] two-way group mean prediction in survreg with three factors

2011-05-04 Thread Pang Du
I'm fitting a regression model for censored data with three categorical predictors, say A, B, C. My final model based on the survreg function is Surv(..) ~ A*(B+C). I know the three-way group mean estimates can be computed using the predict function. But is there any way to obtain two-way group

[R] Fwd: simple question

2011-05-04 Thread Nilaya Sharma
Sorry I had typo in previous email, this typo corrected version: Dear R experts I have simple question, please execuse me: #example data, the real data consists of 2 pairs of variables K1 <- c(1,2,1, 1, 1,1); K2 <- c(1, 1,2,2, 1,2); K3 <- c(3, 1, 3, 3, 1, 3) M1a <- rep( K1, 100); M1b <- re

Re: [R] Storing data from a test as a vector or matrix

2011-05-04 Thread wwreith
I figured out that attributes is the command that I was trying to find. It allowed me to find out that I was needing to use "stats" not "Df" or "Pillai" etc. Following command worked. > S1<-as.vector(S$stats[1,]) However when I try the same thing with summary.aov it is not working. >SA<-summary.

Re: [R] Storing data from a test as a vector or matrix

2011-05-04 Thread wwreith
SA gives the output: Response IPS1 : Df Sum Sq Mean Sq F value Pr(>F) as.factor(WSD) 3 3.3136 1.10455 23.047 5.19e-12 *** Residuals 129 6.1823 0.04793 --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 . . . There are 11 more just

Re: [R] Superscript number before letter

2011-05-04 Thread Janhal
Yes thats it :-) Thank you very much! Janine -- View this message in context: http://r.789695.n4.nabble.com/Superscript-number-before-letter-tp3495577p3495812.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org mailing

[R] Uniform Gaussian Kernel

2011-05-04 Thread blutack
I have a vector with lots of different numbers. I need to make a graph showing the Uniform Distribution of the figures. I have created a graph showing all the different values, but now want individual Gaussian Kernel round each point. This is what I have but each time it comes up with an error as I

[R] Distance measure in heatmaps

2011-05-04 Thread Muhidini
Hi, First- I am relatively new to R and would therefore greatly appreciate any kind of help from you ! I am currently trying to use R to make a heatmap using the Pheatmap package. My question: I want to include a different distance measure than the ones given by the implemented hclust function

[R] Instrumental variable quantile estimation of spatial autoregressive models

2011-05-04 Thread Marie-Line Glaesener
Dear all, I would like to implement a spatial quantile regression using instrumental variable estimation (according to Su and Yang (2007), Instrumental variable quantile estimation of spatial autoregressive models, SMU economics & statistis working paper series, 2007, 05-2007, p.35 ). I am app

[R] Saving Values in a Vector from a For Loop

2011-05-04 Thread blutack
Hi, I have a created a function, but now I need to call it about a hundred times and store the results as a vector. I think doing a for loop would work, but I cant work out how to save the values generated from the function as a vector. Any ideas? Thanks. -- View this message in context: http:/

Re: [R] nls problem with R

2011-05-04 Thread Mike Marchywka
> Date: Wed, 4 May 2011 07:07:44 -0700 > From: sterles...@hotmail.com > To: r-help@r-project.org > Subject: Re: [R] nls problem with R > > Thanks Andrew. > I am sorry for some typos that I omit some numbers of T2. > Based on your suggestion,I think the problem is in the initial values. > And I w

[R] join tables in R

2011-05-04 Thread Roccato Alfredo (UniCredit)
I'd to match-merge 2 tables in such a manner that I keep all the rows in table 1, but not the rows that are in both table 1 and 2. Thank you for your help, Alfredo > master <- data.frame(ID=2001:2011) > train <- data.frame(ID=2004:2006) > valid <- ??? in this example table valid should have th

[R] merging multiple columns from two dataframes

2011-05-04 Thread 1Rnwb
Hello, I have data in a dataframe with 139104 rows which is multiple of 96x1449. i have a phenotype file which contains the phenotype information for the 96 samples. the snp name is repeated 1449X96 samples. I haveto merge the two dataframes based on sid and sen. this is how my two dataframes loo

[R] Regexp question

2011-05-04 Thread johannes rara
I have a string like this st <- "SELECT COUNT(empid), COUNT(mgrid), COUNT(empname), COUNT(salary), FROM Employees" How can I remove the last comma before the FROM statement? -J __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo

Re: [R] combine lattice plot and standard R plot

2011-05-04 Thread Scott Chamberlain
What about the example in gridExtra package: require(ggplot2); require(lattice); require(gridExtra) grid.arrange(qplot(1:10), xyplot(1:10~1:10), tableGrob(head(iris)), nrow=2, as.table=TRUE, main="test main", sub=textGrob("test sub", gp=gpar(font=2))) On Wednesday, May 4, 2011 at 1:44 PM, Jonatha

Re: [R] combine lattice plot and standard R plot

2011-05-04 Thread Jonathan Daily
If you read the help documentation, lattice is not really compatible with standard graphics. library("lattice") ?lattice 2011/5/4 Lucia Cañas : > Dear R users, > > I would like to combine lattice plot (xyplot) and standard R plot (plot and > plotCI) in an unique figure. > > I use the function "p

[R] combine lattice plot and standard R plot

2011-05-04 Thread Lucia Cañas
Dear R users, I would like to combine lattice plot (xyplot) and standard R plot (plot and plotCI) in an unique figure. I use the function "par()" to combine plot and plotCI and I use the function "print()" to combine xyplot. I tried to use these functions to combine xyplot and plotCI and plots

Re: [R] tryCatch?

2011-05-04 Thread Mikkel Grum
Beautiful Prof. This worked: for (i in seq(along = tbl$key)) { if (inherits(try(sqlSave(pg, tbl[i, ], "tbl", append = TRUE, rownames = FALSE), silent = TRUE), "try-error", TRUE)) sqlUpdate(pg, tbl[i, ], index = "key1") ) } --- On Wed, 5/4/11, Prof Brian Ripley wrote:

Re: [R] natural cubic splines

2011-05-04 Thread David Winsemius
On May 4, 2011, at 8:34 AM, Joshua Wiley wrote: Hi Ondrej, What documentation have you looked at? Does this help at all? require(splines) ?ns ## one example summary(lm(y ~ ns(x, df = 3), data = data.frame(y = runif(100), x = rbinom(100, 9, .25)^2))) ## built in examples example(ns) Also,

Re: [R] Superscript number before letter

2011-05-04 Thread David Winsemius
On May 4, 2011, at 8:43 AM, Richard M. Heiberger wrote: > David, > > That is not clear from the original request. The request was for > {}^18*O I see your point. The use of the phantom would also allow a pre- superscript with no preceding "delta". Thanks for bearing with my obtuseness. --

Re: [R] problem with package "adapt" for R in Mac

2011-05-04 Thread Uwe Ligges
On 04.05.2011 19:22, Matías Ramírez Salgado wrote: The package is called: "adapt" (adapt_1.0-4.tgz is the version of package for mac, search it in google) is not in all CRAN servers, because is an old package. Yes, there is a reason why the package was archived (which I said already!): It does

[R] problem with package "adapt" for R in Mac

2011-05-04 Thread Matías Ramírez Salgado
The package is called: "adapt" (adapt_1.0-4.tgz is the version of package for mac, search it in google) is not in all CRAN servers, because is an old package. I tried to install it from package manager as CRAN (binaries), Other repositorys, and as a local source packages, and allways returns this

Re: [R] tryCatch?

2011-05-04 Thread Prof Brian Ripley
Start with try(): you may find it easier to understand. if(inherits(try(), "try-error")) so in your case if(inherits(try(sqlSave(pg, tbl[i, ], "tbl", append = TRUE, rownames = FALSE sqlUpdate(pg, tbl[i, ], index = "key") or some such. On Wed, 4 May 2011, Mikke

Re: [R] Simple loop

2011-05-04 Thread Petr Savicky
On Wed, May 04, 2011 at 08:52:07AM -0700, William Dunlap wrote: > > -Original Message- > > From: r-help-boun...@r-project.org > > [mailto:r-help-boun...@r-project.org] On Behalf Of Petr Savicky > > Sent: Wednesday, May 04, 2011 12:51 AM > > To: r-help@r-project.org > > Subject: Re: [R] Sim

[R] tryCatch?

2011-05-04 Thread Mikkel Grum
I would like to do inserts into a database table, but do updates in the fairly rare cases in which the inserts fail. I thought tryCatch might be the way to do it, but I honestly do not understand the help file for tryCatch at all. I thought something like this might work: for (i in seq(along = t

Re: [R] Error Rscript: No such file or directory

2011-05-04 Thread xavier abulker
Thanks a lot Duncan, I have removed the system variable incorrectly setup and it works now. Have a nice day Xavier From: Duncan Murdoch Cc: r-help@r-project.org Sent: Wed, 4 May, 2011 17:37:55 Subject: Re: [R] Error Rscript: No such file or directory On 04/05

Re: [R] Error Rscript: No such file or directory

2011-05-04 Thread Duncan Murdoch
On 04/05/2011 11:49 AM, xavier abulker wrote: Hello, I'm trying to build a simple cpp file using the R CMD SHLIB command and I always receive the same error message: cygwin warning: MS-DOS style path detected: C:/PROGRA~1/R/R-212~1.1/etc/i386/Makeconf Preferred POSIX equivalent is: /cygdri

[R] xtable without a loop alongside a ggplot

2011-05-04 Thread Justin Haynes
I would like to create a table of my points and identify which 'quadrant' of a plot they are in with the 'origin' at the means. the kicker is i would like to display it right next to or below a ggplot of the data. Maybe xtable isnt the right thing to use, but its the only thing i can think of. A

[R] Panels order in lattice graphs

2011-05-04 Thread Cristina Silva
Hi all, In lattice graphs, panels are drawn from left to right and bottom to top. The flag "as.table=TRUE" changes to left to right and top to bottom. Is there any way to change to first top to bottom and then left to right? didn´t find anything neither in Help pages nor Lattice book. Cristi

Re: [R] Simple loop

2011-05-04 Thread William Dunlap
> -Original Message- > From: r-help-boun...@r-project.org > [mailto:r-help-boun...@r-project.org] On Behalf Of Petr Savicky > Sent: Wednesday, May 04, 2011 12:51 AM > To: r-help@r-project.org > Subject: Re: [R] Simple loop > > On Tue, May 03, 2011 at 12:04:47PM -0700, William Dunlap wrote

Re: [R] Superscript number before letter

2011-05-04 Thread Richard M. Heiberger
David, That is not clear from the original request. The request was for {}^18*O It wasn't for delta^18 Therefore I put the space there to be sure that the 18 was seen as pre-superscript of O, not as a post-superscript of delta. I probably should also have used ~ as plot(1:10, xlab=expression(d

[R] Error Rscript: No such file or directory

2011-05-04 Thread xavier abulker
Hello, I'm trying to build a simple cpp file using the R CMD SHLIB command and I always receive the same error message: cygwin warning: MS-DOS style path detected: C:/PROGRA~1/R/R-212~1.1/etc/i386/Makeconf Preferred POSIX equivalent is: /cygdrive/c/PROGRA~1/R/R-212~1.1/etc/i386/Makeconf

Re: [R] natural cubic splines

2011-05-04 Thread Joshua Wiley
Hi Ondrej, What documentation have you looked at? Does this help at all? require(splines) ?ns ## one example summary(lm(y ~ ns(x, df = 3), data = data.frame(y = runif(100), x = rbinom(100, 9, .25)^2))) ## built in examples example(ns) Also, I am very fond of the book, Modern Applied Statisti

[R] natural cubic splines

2011-05-04 Thread Ondřej Mikula
Dear R-helpers, I need to fit natural cubic spline with specified number of knots. I expected 'splines' package will be helpful, but I am confused by its help. Is more detailed documentation available for it or could you recommend another R function? Best regards Ondrej Mikula

[R] hurdle, simulated power

2011-05-04 Thread David Atkins
Hi all-- We are planning an intervention study for adolescent alcohol use, and I am planning to use simulations based on a hurdle model (using the hurdle() function in package pscl) for sample size estimation. The simulation code and power code are below -- note that at the moment the "powe

Re: [R] create a folder with mode '0777'

2011-05-04 Thread Prof Brian Ripley
Please read the comments in the help about umask (and in the posting guide about the 'at a minimum' information required in postings, for the details do depend on the version of R and it seems yours is not current). In R 2.13.0: ‘dir.create’ creates the last element of the path, unless

Re: [R] Is this confict of different versions of R or something else?

2011-05-04 Thread David Winsemius
On May 4, 2011, at 3:44 AM, TheSavageSam wrote: Hello! I have had some problems lately with use of R at home and school. At my home laptop (Ubuntu linux 64bit & R-2.12.0) R works just fine. But when I take my codes to school(Windows XP 32bit & R-2.10.1 I think) and run those there those

Re: [R] Superscript number before letter

2011-05-04 Thread David Winsemius
On May 4, 2011, at 7:28 AM, Richard M. Heiberger wrote: Dos this do what you want? plot(1:10, xlab=expression(delta*{}^18*"O" * " VSMOW [‰]")) The specific is to put an empty item there to hold the superscript. I do not think that is necessary: plot(1:10, xlab=expression(delta^18*O~VSMOW["

Re: [R] Storing data from a test as a vector or matrix

2011-05-04 Thread Ivan Calandra
Hi, I would suggest you to check the structure of your summary object with str(), like this: S <- summary(M, test="Pillai") str(S) You will then see how to access each element of it. If you cannot manage to do it yourself, then provide an example, or at least the output from str(s). By the wa

Re: [R] Format ddmmYYYY in date

2011-05-04 Thread David Winsemius
On May 4, 2011, at 5:57 AM, Ben Bolker wrote: smoff boku.ac.at> writes: My problem is that I have a table containing dates in the first column of 10 years. These dates have the format ddmm at least in the csv- file. After importing the file using read.table() R deletes the first charac

Re: [R] what happens when I store linear models in an array?

2011-05-04 Thread David Winsemius
On May 4, 2011, at 4:49 AM, Andrew D. Steen wrote: I've got a bunch of similar datasets, all of which I've fit to linear models. I'd like to easily create arrays of a specific parameter from each linear model (e.g., all of the intercepts in one array). I figured I'd put the model objects

Re: [R] Format ddmmYYYY in date

2011-05-04 Thread David Winsemius
On May 4, 2011, at 5:57 AM, Ben Bolker wrote: smoff boku.ac.at> writes: My problem is that I have a table containing dates in the first column of 10 years. These dates have the format ddmm at least in the csv- file. After importing the file using read.table() R deletes the first charac

Re: [R] create a folder with mode '0777'

2011-05-04 Thread Patrick Breheny
Linux systems have a user mask that limits the file mode creation possibilities of any processes launched from that shell. If you check your /etc/profile file, you will see the line umask 022 This prevents you by default from creating files with write access for everyone except the user. In

Re: [R] Simple General Statistics and R question (with 3 line example) - get z value from pairwise.wilcox.test

2011-05-04 Thread peter dalgaard
On May 4, 2011, at 15:11 , JP wrote: > Peter thanks for the fantastically simple and understandable explanation... > > To sum it up... to find the z values of a number of pairwise wilcox > tests do the following: > > # pairwise tests with bonferroni correction > x <- pairwise.wilcox.test(a, b,

[R] Odp: Is this confict of different versions of R or something else?

2011-05-04 Thread Petr PIKAL
Hi r-help-boun...@r-project.org napsal dne 04.05.2011 12:44:04: > TheSavageSam > Odeslal: r-help-boun...@r-project.org > > > Hello! > > I have had some problems lately with use of R at home and school. At my home > laptop (Ubuntu linux 64bit & R-2.12.0) R works just fine. But when I take m

Re: [R] Superscript number before letter

2011-05-04 Thread Richard M. Heiberger
Dos this do what you want? plot(1:10, xlab=expression(delta*{}^18*"O" * " VSMOW [‰]")) The specific is to put an empty item there to hold the superscript. On Wed, May 4, 2011 at 9:37 AM, Janhal wrote: > Salut, > I have been struggling to superscript the 18 before the O without ^ visible > an

Re: [R] nls problem with R

2011-05-04 Thread Ravi Varadhan
In addition to the suggestion about finding a good initial value, you should also scale your response V2 (and, of course, V0). Divide V2 by 10^4, for example. Now your V0 should also be scaled by this factor. This would likely help with convergence. Ravi. ---

Re: [R] Watts Strogatz game

2011-05-04 Thread Gábor Csárdi
On Wed, May 4, 2011 at 9:28 AM, Robert Baer wrote: >>> I have a erdos-renyi game with 6000 nodes and probability 0.003. >>> >>> g1 = erdos.renyi.game(6000, 0.003) >>> >>> How to create a Watts Strogatz game with the same probability. >>> >>> g1 = watts.strogatz.game(1, 6000, ?, ?) >>> What should

Re: [R] nls problem with R

2011-05-04 Thread sterlesser
Thanks Andrew. I am sorry for some typos that I omit some numbers of T2. Based on your suggestion,I think the problem is in the initial values. And I will read more theory about the non-linear regression. -- View this message in context: http://r.789695.n4.nabble.com/nls-problem-with-R-tp3494454p

Re: [R] ID parameter in model

2011-05-04 Thread Göran Broström
On Wed, May 4, 2011 at 2:18 PM, Mike Harwood wrote: > Thank you, Goran.  Please see the package details below: Thanks, I have uploaded a corrected version of eha to CRAN. Should be available soon. Göran [...] __ R-help@r-project.org mailing list http

Re: [R] nls problem with R

2011-05-04 Thread sterlesser
Thanks Ruben. Your suggestion about more deeper analysis about the model itself is really helpful. I am trying out some new initial values based on the analysis of the special T2 in the model. -- View this message in context: http://r.789695.n4.nabble.com/nls-problem-with-R-tp3494454p3495663.html

[R] create a folder with mode '0777'

2011-05-04 Thread Xian Zhang
Dear list, I am trying to create a folder structure, say 'test/sub', and set the folder and sub folder to be writable to everyone. By default dir.create('test/sub', recursive=TRUE, mode='0777') creates folders with mode: drwxr-xr-x After Sys.chmod('test/sub',mode='0777') The folder 'test' is

Re: [R] Unexp. behavior from boot with multiple statistics

2011-05-04 Thread algorimancer
Thanks, that clears things up quite a bit. Now I'm left wondering why there is so much bias, but that's a separate issue. -- View this message in context: http://r.789695.n4.nabble.com/Unexp-behavior-from-boot-with-multiple-statistics-tp3493300p3495590.html Sent from the R help mailing list arch

  1   2   >