Re: [R] Connect R and Lyx in UBUNTU

2012-10-10 Thread Rainer M Krug
-BEGIN PGP SIGNED MESSAGE- Hash: SHA1 I have pointed out already: This is a LyX question - please ask on their mailing list (http://www.lyx.org/MailingLists#toc2 and http://dir.gmane.org/gmane.editors.lyx.general) There are many users who use LyX / sweave or knitr / R under Ubuntu! Rai

Re: [R] Exporting summary plm results to latex

2012-10-10 Thread Duncan Mackay
Hi Sebastian I think I found the package by accident when I did a search of the Cran package page for"latex" but did not use it as it could not do a very particular problem. If there was no other alternative use the add.to.row argument of xtable A while ago I needed to add some info from the

Re: [R] Connect R and Lyx in UBUNTU

2012-10-10 Thread Yihui Xie
It is actually much easier to do it under Ubuntu; see a video here: http://yihui.name/knitr/demo/lyx/ If you want to use Sweave instead of knitr, there is also a module for it. The official documentation is here: - https://github.com/downloads/yihui/lyx/sweave.pdf - https://github.com/downloads/y

Re: [R] Connect R and Lyx in UBUNTU

2012-10-10 Thread ATANU
By "Connect" I meant to say that I was able to write code chunks in LYX and compile them within LYX( using R) to produce results along with other stuffs. There are many tutorials available for doing this under Windows but I could not solve the problem for linux (UBUNTU). -Atanu -- View thi

Re: [R] replacing ugly for loops

2012-10-10 Thread Bert Gunter
I am not sure you have expressed what you wanjt to do correctly. See inline: On Wed, Oct 10, 2012 at 9:10 PM, andrewH wrote: > I have a couple of hundred American Community Survey Summary Files files > containing rectangular arrays of data, mainly though not exclusively > numeric. Each file is r

[R] replacing ugly for loops

2012-10-10 Thread andrewH
I have a couple of hundred American Community Survey Summary Files files containing rectangular arrays of data, mainly though not exclusively numeric. Each file is referred to as a sequence (henceforth "seq"). From these files I am trying to extract particular subsets (tables) consisting of a set

Re: [R] svyplot

2012-10-10 Thread Muhuri, Pradip (SAMHSA/CBHSQ)
Dear Anthony, You have been so so helpful! The par() example code has worked well. Thanks, Pradip From: Anthony Damico [ajdam...@gmail.com] Sent: Wednesday, October 10, 2012 5:25 PM To: Muhuri, Pradip (SAMH/CBHSQ) Cc: R help Subject: Re: [R] svyplot ht

Re: [R] multiple t-tests across similar variable names

2012-10-10 Thread arun
HI, If you have a lot of variables and in no order, then it would be better to order the data by column names. For e.g. set.seed(432) dat2<-data.frame(apple_pre=sample(10:20,5,replace=TRUE),orange_post=sample(18:28,5,replace=TRUE),banana_pre=sample(25:35,5,replace=TRUE),apple_post=sample(20:30,5,

Re: [R] Converting factors to bounded random numerical data

2012-10-10 Thread David L Carlson
I forgot that you were looking for whole numbers: > Zfact <- data.frame(factor=LETTERS[1:3], low=c(.5, 4.5, 9.5), high=c(4.499, 9.499, 12.499)) > Zfact factor low high 1 A 0.5 4.499 2 B 4.5 9.499 3 C 9.5 12.499 > set.seed(42) > a$Zval <- round(runif(nrow(a), Zfact$low[as

Re: [R] Converting factors to bounded random numerical data

2012-10-10 Thread David L Carlson
This should work: > a <- data.frame(X=rnorm(25, 20, 4), Y=rnorm(25, 15, 3), Z=sample(c(LETTERS[1:3]), 25, replace=TRUE)) > head(a) X Y Z 1 17.77449 14.425221 B 2 19.95400 13.408439 A 3 13.40162 12.219984 A 4 15.89822 19.214026 B 5 18.55717 14.568691 B 6 19.86619 11.60

Re: [R] multiple t-tests across similar variable names

2012-10-10 Thread arun
HI, A typo in my solution: row.names(res2)<-unlist(unique(lapply(strsplit(colnames(dat2),"_"),`[`,1)))         A.K. - Original Message - From: "Nundy, Shantanu" To: "r-help@r-project.org" Cc: Se

Re: [R] multiple t-tests across similar variable names

2012-10-10 Thread arun
HI, May be this helps you. set.seed(1) dat2<-data.frame(apple_pre=sample(10:20,5,replace=TRUE),banana_pre=sample(25:35,5,replace=TRUE),apple_post=sample(20:30,5,replace=TRUE),banana_post=sample(40:50,5,replace=TRUE)) list2<-list(dat2[regmatches(colnames(dat2),regexpr("apple.*",colnames(dat2)))],dat

Re: [R] glmmPQL and spatial correlation

2012-10-10 Thread Ben Bolker
Luis Huckstadt ucsc.edu> writes: > I'm running into some computer issues when trying to run a binomial model > for spatially correlated data using glmmPQL and was wondering if anyone > could help me out. > My whole dataset consists of about 300,000 points for which I have a suite > of environment

Re: [R] Contacting Delphi ??

2012-10-10 Thread Ben Bolker
Rui Barradas sapo.pt> writes: > It's in the source code for `?`, file src/library/utils/R/question.R, > lines 32 to 35. > See: https://github.com/wch/r-source/commit/34b3998c928fbf50e24ab0e33c7d72ab8c944330 __ R-help@r-project.org mailing list htt

Re: [R] practical to loop over 2million rows?

2012-10-10 Thread David Winsemius
On Oct 10, 2012, at 6:45 PM, jim holtman wrote: > This is a classic example from my tag line: > > Tell me what you want to do, not how you want to do it. > > For example you provided no information as to what the objects were. > I hope that 'stratID' is at least of length one greater than 'x'

Re: [R] Exporting summary plm results to latex

2012-10-10 Thread Sebastian Barfort
Hi David, this looks very promising, but I am afraid I can't see clearly how this would work. If possible, would you mind a short explanation? perhaps using the attached exampled? Thank you for your advice, Sebastian On Oct 10, 2012, at 9:16 PM, David Winsemius wrote: > > On Oct 10, 2012

Re: [R] "nlmnib" Package + Hessian Output

2012-10-10 Thread Ben Bolker
nserdar hotmail.com> writes: > > > I need a "Hessian" matrix in "nlmnib" package to discuss whether parameters > are significant or not. > > Please let me know how to obtain hessian matrix and how to evaluate the > significancy of parameters. > You can get a finite-difference approximatio

Re: [R] practical to loop over 2million rows?

2012-10-10 Thread jim holtman
This is a classic example from my tag line: Tell me what you want to do, not how you want to do it. For example you provided no information as to what the objects were. I hope that 'stratID' is at least of length one greater than 'x' based on your loops. Also on the last iteration you are tryin

Re: [R] multiple t-tests across similar variable names

2012-10-10 Thread Rui Barradas
Hello, Could you post a data example? Using, with data.frame named 'dat' dput( head(dat, 30) ) # paste the output of this in a post I have written code that creates pairs pre/post columns but it can't really be tested. Hope this helps, Rui Barradas Em 11-10-2012 00:09, Nundy, Shantanu escr

Re: [R] Converting factors to bounded random numerical data

2012-10-10 Thread jim holtman
You need to include how many random numbers you want to create; you are asking for only 1, that is why they are are the same. You probably want something like: a$Z<-ifelse(a$Z=="L" ,sample(1:4, nrow(a), TRUE) ,ifelse(a$Z=="M" ,sample(5:9, nrow(a), TRUE) ,ifelse(a$Z=="U"

Re: [R] Converting factors to bounded random numerical data

2012-10-10 Thread Jeff Newmiller
Please don't double post. --- Jeff NewmillerThe . . Go Live... DCN:Basics: ##.#. ##.#. Live Go... Live: OO#.. Dead: OO#.. Playin

Re: [R] Exporting summary plm results to latex

2012-10-10 Thread David Winsemius
On Oct 10, 2012, at 4:45 PM, Sebastian Barfort wrote: > I am also interested in the standard errors, but beneath not next to the > point estimates which is standard in the xtable package. Last year Mark Difford offered code to do that. From: Mark Difford Subject:Re:

Re: [R] Converting factors to bounded random numerical data

2012-10-10 Thread William Dunlap
sample(Range,1) generates 1 random number. Change that (in all 3 places) to sample(Range, length(a$Z), replace=TRUE) to get length(a$Z) random numbers in the same range. There are other ways to do this with less typing. Bill Dunlap Spotfire, TIBCO Software wdunlap tibco.com > -Original Me

[R] Converting factors to bounded random numerical data

2012-10-10 Thread KoopaTrooper
I have a data set (a) with three columns (X,Y,Z). The first 2 columns are numeric. The third (Z) is a factor with three levels A,B,C. I want to turn each A into a different random number between 1 and 4, each B into a different random number between 5 and 8, etc. I tried this: a$Z<-ifelse(a$Z=="L

[R] Converting factors to bounded random numerical data

2012-10-10 Thread KoopaTrooper
I have a data set (a) with 3 columns (X,Y,Z). The first 2 columns are numerical. The third column (Z) is a factor with three levels ("A","B","C"). What I want to do is turn each of the "A's" into different random numbers between 1 and 4, "B's" into a random number between 5 and 8, etc. I tried thi

Re: [R] Exporting summary plm results to latex

2012-10-10 Thread Sebastian Barfort
I am also interested in the standard errors, but beneath not next to the point estimates which is standard in the xtable package. If you by any chance remember the name of the package or how to do it that would be much appreciated! Cheers, Sebastian On Oct 10, 2012, at 7:10 PM, Duncan Mackay

Re: [R] Exporting summary plm results to latex

2012-10-10 Thread Sebastian Barfort
Hi A.K, thanks for the suggestion. Unfortunately, it's not really what I'm looking for. I'm looking for a package or some function that can create the point estimate in the first row with the standard error in parenthesis in the second which is not what you normally get from the xtable package.

Re: [R] multiple t-tests across similar variable names

2012-10-10 Thread Nundy, Shantanu
Hi everyone- I have a dataset with multiple "pre" and "post" variables I want to compare. The variables are named "apple_pre" or "pre_banana" with the corresponding post variables named "apple_post" or "post_banana". The variables are in no particular order. apple_pre orange_pre orange_post pr

[R] model selection with spg and AIC (or, convert list to fitted model object)

2012-10-10 Thread Adam Zeilinger
Dear R Help, I have two nested negative log-likelihood functions that I am optimizing with the spg function [BB package]. I would like to perform model selection on these two objective functions using AIC (and possibly anova() too). However, the spg() function returns a list and I need a fi

Re: [R] smoothScatter plot

2012-10-10 Thread JiangZhengyu
Hi John,I installed & but somehow it did not work on my computer, while I tried another computer - it works. Thanks for all your communications.Best,Zhengyu Date: Tue, 9 Oct 2012 05:16:59 -0800 From: jrkrid...@inbox.com Subject: RE: [R] smoothScatter plot To: zhyjiang2...@hotmail.com Glad i

Re: [R] Exporting summary plm results to latex

2012-10-10 Thread Duncan Mackay
Hi If you just want the coefficients. xtable(summary(fe)$coef) % latex table generated in R 2.15.1 by xtable 1.7-0 package % Thu Oct 11 09:04:59 2012 \begin{table}[ht] \begin{center} \begin{tabular}{r} \hline & Estimate & Std. Error & t-value & Pr($>$$|$t$|$) \\ \hline x & 0.12 & 0.07

Re: [R] Send Email from R

2012-10-10 Thread Kjetil Halvorsen
see inline below. On Tue, Oct 9, 2012 at 10:35 AM, Rantony wrote: > Hi, > > i wanted to send a mail from R (using Eclips). > Currently i installed the following packages > base64_1.1 > sendmailR_1.1-1 > mail_1.0.tar.gz > Rmail_1.1.tar.gz > > But while installing package some error was occuring.

Re: [R] get: problem with environments

2012-10-10 Thread Thomas Lumley
On Thu, Oct 11, 2012 at 11:18 AM, R. Michael Weylandt wrote: > Thanks Prof Lumley, > > I'm still not sure how this gets to the call stack of 5 or 6 the OP > reported or the difference between GUI & Terminal. Any thoughts there? I don't see how the terminal version is getting 5 and 6 rather than

Re: [R] How to replicate SAS by group processing in R

2012-10-10 Thread Ista Zahn
Hi, here is one way using ddply (from the plyr package): dat <- read.table(text="tdate stock_symbolexpiration strike 9/11/2012 C 9/16/201211 9/11/2012 C 9/16/201212 9/11/2012 C 9/16/201

[R] GAM without intercept

2012-10-10 Thread SAEC
Hi everybody, I am trying to fit a GAM model without intercept using library mgcv. However, the result has nothing to do with the observed data. In fact the predicted points are far from the predicted points obtained from the model with intercept. For example: #First I generate some simulated

Re: [R] How to replicate SAS by group processing in R

2012-10-10 Thread arun
Hi, You can also try with aggregate() or ddply() dat2<-aggregate(dat,list(dat$stock_symbol,dat$tdate),FUN=function(x) head(x,1))  dat2[,3:6] #  tdate stock_symbol expiration strike #1 9/11/2012    C  9/16/2012 11 #2 9/12/2012    C  9/16/2012 14 library(plyr)  ddply(da

Re: [R] lm on matrix data

2012-10-10 Thread R. Michael Weylandt
On Wed, Oct 10, 2012 at 3:35 PM, Baoqiang Cao wrote: > Hi, > > I have a question about using lm on matrix, have to admit it is very > trivial but I just couldn't find the answer after searched the mailing > list and other online tutorial. It would be great if you could help. > > I have a matrix "t

Re: [R] warning in summary(aov())

2012-10-10 Thread R. Michael Weylandt
On Mon, Oct 8, 2012 at 8:29 PM, Jhope wrote: > Hi Michael, > > I am very willing to subscribe to where you think is appropriate but I > thought I was subscribed. > > I have a password that I can log into Nabble with and I receive update > emails from the R-help forum constantly in to my inbox [R]

Re: [R] Error in matrix (unlist(value, recursive = FALSE, use.names = FALSE), nrow = nr, : attempt to set an attribute on NULL

2012-10-10 Thread R. Michael Weylandt
On Tue, Oct 9, 2012 at 10:54 PM, marsyxp wrote: > I am using Donlp2 package to solve a non-linear problem, but there's an error > I always meet: > > Error in matrix(unlist(value, recursive = FALSE, use.names = FALSE), nrow = > nr, : > attempt to set an attribute on NULL > > I have been sufferin

Re: [R] get: problem with environments

2012-10-10 Thread R. Michael Weylandt
On Wed, Oct 10, 2012 at 10:04 PM, Thomas Lumley wrote: > On Thu, Oct 11, 2012 at 9:05 AM, R. Michael Weylandt > wrote: >> Can someone more capable than I help Martin out with this? I'm feeling >> out of my league (that or I've missed something obvious) >> >> Shot in the dark: you aren't running t

Re: [R] histogram & corelation plots

2012-10-10 Thread David L Carlson
As for interactive plots, there is also aplpack (sliders for changing the width of histograms, and density plots; iPlots; and iWebPlots. -- David L Carlson Associate Professor of Anthropology Texas A&M University College Station, TX 77843-4352 > -

Re: [R] How to replicate SAS by group processing in R

2012-10-10 Thread David Winsemius
On Oct 10, 2012, at 11:09 AM, ramoss wrote: > Hello, > > I am trying to re-code all my programs from SAS into R. > > In SAS I use the following code: > > proc sort data=upper; > by tdate stock_symbol expire strike; > run; > data upper1; > set upper; > by tdate stock_symbol expire strike;

Re: [R] svyplot

2012-10-10 Thread Anthony Damico
https://stat.ethz.ch/pipermail/r-help/2012-October/324944.html On Wed, Oct 10, 2012 at 5:05 PM, Muhuri, Pradip (SAMHSA/CBHSQ) < pradip.muh...@samhsa.hhs.gov> wrote: > Hello, > > Using the svyplot () function, I have plotted four graphs that are saved > in four different .png files. > > I am look

Re: [R] practical to loop over 2million rows?

2012-10-10 Thread David Winsemius
On Oct 10, 2012, at 1:31 PM, Jay Rice wrote: > New to R and having issues with loops. I am aware that I should use > vectorization whenever possible and use the apply functions, however, > sometimes a loop seems necessary. > > I have a data set of 2 million rows and have tried run a couple of lo

Re: [R] practical to loop over 2million rows?

2012-10-10 Thread Joshua Wiley
Hi Jay, A few comments. 1) As you know, vectorize when possible. Even if you must have a loop, perhaps you can avoid nested loops or at least speed each iteration. 2) Write your loop in a function and then byte compile it using the cmpfun() function from the compiler package. This can help dram

[R] svyplot

2012-10-10 Thread Muhuri, Pradip (SAMHSA/CBHSQ)
Hello, Using the svyplot () function, I have plotted four graphs that are saved in four different .png files. I am looking for examples how to redraw the same four graphs within grid viewports so that they stay together on a page. The goal is to create one .png file that will include all four

Re: [R] get: problem with environments

2012-10-10 Thread Thomas Lumley
On Thu, Oct 11, 2012 at 9:05 AM, R. Michael Weylandt wrote: > Can someone more capable than I help Martin out with this? I'm feeling > out of my league (that or I've missed something obvious) > > Shot in the dark: you aren't running this in some sort of debug mode, are you? > > RMW > > On Sun, Oct

Re: [R] quantreg Wald-Test

2012-10-10 Thread Massimo Schiavo
Hi Stephan, Now I have your same problem on Granger causality test, have you got some news about it? The anova.rq is the right one to do it? Regards, Massimo -- View this message in context: http://r.789695.n4.nabble.com/quantreg-Wald-Test-tp4638198p4645773.html Sent from the R help mailing lis

[R] "optim" and "nlminb"

2012-10-10 Thread nserdar
#optim package estimate<-optim(init.par,Linn,hessian=TRUE, method=c("L-BFGS-B"),control = list(trace=1,abstol=0.001),lower=c(0,0,0,0,-Inf,-Inf,-Inf,-Inf,-Inf,-Inf,-Inf,-Inf,-Inf),upper=c(1,1,1,1,Inf,Inf,Inf,Inf,Inf,Inf,Inf,Inf,Inf)) #nlminb package estimate<-nlminb(init.par,Linn,gr=NULL,hessian=T

[R] "nlmnib" Package + Hessian Output

2012-10-10 Thread nserdar
I need a "Hessian" matrix in "nlmnib" package to discuss whether parameters are significant or not. Please let me know how to obtain hessian matrix and how to evaluate the significancy of parameters. Regards Serdar -- View this message in context: http://r.789695.n4.nabble.com/nlmnib-Packag

[R] practical to loop over 2million rows?

2012-10-10 Thread Jay Rice
New to R and having issues with loops. I am aware that I should use vectorization whenever possible and use the apply functions, however, sometimes a loop seems necessary. I have a data set of 2 million rows and have tried run a couple of loops of varying complexity to test efficiency. If I do a v

[R] How to replicate SAS by group processing in R

2012-10-10 Thread ramoss
Hello, I am trying to re-code all my programs from SAS into R. In SAS I use the following code: proc sort data=upper; by tdate stock_symbol expire strike; run; data upper1; set upper; by tdate stock_symbol expire strike; if first.expire then output; rename strike=astrike; run; on the

[R] combine unadjusted and adjusted forest plots

2012-10-10 Thread Romita Mukerjee, M.D.
Hello, I am learning to use the metafor package to conduct meta-regression analyses for a systematic review on multidisciplinary care interventions in chronic kidney disease. For the forest plots, I can't figure out how to plot unadjusted and adjusted models on the same plot. From top to bott

Re: [R] get: problem with environments

2012-10-10 Thread R. Michael Weylandt
Can someone more capable than I help Martin out with this? I'm feeling out of my league (that or I've missed something obvious) Shot in the dark: you aren't running this in some sort of debug mode, are you? RMW On Sun, Oct 7, 2012 at 5:10 PM, Martin Ivanov wrote: > Thank You very much for Your

Re: [R] histogram & corelation plots

2012-10-10 Thread R. Michael Weylandt
On Wed, Oct 10, 2012 at 4:47 PM, sagarnikam123 wrote: > i want to show histogram in innovative way (good if interactive) in html > report..is threre any tutorials/hint > also there is any package to show correalation plot The histogram is well defined, so I'm not really sure how much innovation y

Re: [R] own function: computing time

2012-10-10 Thread William Dunlap
Your original method would be the following function f <- function (x, y) { xy <- cbind(x, y) outside <- function(z) { !any(x > z[1] & y > z[2]) } j <- apply(xy, 1, outside) which(j) } and the following one quickly computes the same thing as the above as long as there

Re: [R] package request

2012-10-10 Thread Greg Snow
There are packages for big data analysis, which is best depends on what you want to do. The High Performance Computing task view on CRAN has a section on packages that deal with big data which gives some more detail and may help you choose which package(s) to use. On Wed, Oct 10, 2012 at 12:36 AM

Re: [R] own function: computing time

2012-10-10 Thread Jan van der Laan
Did not see a simple way to make it faster. However, this is a piece of code which can be made to run much faster in C. See below. I don't know if you are familiar with running c-code from R. If not, the official documentation is in the R Extensions manual. However, this is not the most easy

Re: [R] Numeric "Label" of Factor value?

2012-10-10 Thread Rui Barradas
Sorry, not one of my days. Forgot to Cc the list. Rui barradas Em 10-10-2012 20:28, Sarah Goslee escreveu: Sent just to me? On Wed, Oct 10, 2012 at 3:26 PM, Rui Barradas wrote: You're right, apologies to the op and the list. I was thinking of the more complicated as.numeric(levels(abc)[abc[1

Re: [R] Contacting Delphi ??

2012-10-10 Thread R. Michael Weylandt
I think David was saying the humor was unsuccessful (i.e., Christian didn't get the joke), not that the code was. Cheers, RMW On Wed, Oct 10, 2012 at 7:16 PM, Rui Barradas wrote: > Unsuccessfull? Why unsuccessfull? Have you noticed the call to rpois? It > beats runif by potential infinity. > > R

Re: [R] Exporting summary plm results to latex

2012-10-10 Thread arun
HI, May be you can use library(texreg): library(plm) #generating some data x <- rnorm(270) y <- rnorm(270) t <- rep(1:3,30) i <- rep(1:90, each=3) data <- data.frame(i,t,x,y) fe <- plm(y~x,data=data,model="within") summary(fe) library(texreg) fe1<-extract.plm(fe) #extract the plm object libra

Re: [R] Numeric "Label" of Factor value?

2012-10-10 Thread Sarah Goslee
Rui, that doesn't answer the question as I understood it: Your suggestion returns the numeric value of the second value of the levels: > as.numeric(levels(abc)[2]) [1] 3 But I read the question as wanting the numeric value of the second element of abc: > as.numeric(as.character(abc[2])) [1] 2 O

Re: [R] Numeric "Label" of Factor value?

2012-10-10 Thread Rui Barradas
Hello, Try instead ?levels abc <- factor(c(2,2,3,4,7,7)) as.numeric(levels(abc)[1]) Hope this helps, Rui Barradas Em 10-10-2012 19:39, Brigid Mooney escreveu: Sorry, I'm sure I'm not using the appropriate vocab here, which is undoubtedly why I can't seem to find a fix to this (hopefully very

Re: [R] Numeric "Label" of Factor value?

2012-10-10 Thread Sarah Goslee
Hi Brigid, as.numeric() extracts the index of the factor level, which is the way R handles the likelihood that a factor is not actually numeric. Try: > as.numeric(as.character(abc[1])) [1] 2 and see also ?factor particularly the section on the interpretation of a factor. Sarah On Wed, Oct 10,

[R] Numeric "Label" of Factor value?

2012-10-10 Thread Brigid Mooney
Sorry, I'm sure I'm not using the appropriate vocab here, which is undoubtedly why I can't seem to find a fix to this (hopefully very easy) problem. Suppose you have a factor abc <- factor(c(2,2,3,4,7,7)) And you want to know what the number in the nth spot in that would be abc[1] [1] 2 Levels:

Re: [R] reading in a (very simple) list from a file

2012-10-10 Thread David L Carlson
It might help if you would bore us with at least one or two of the things you have tried. It seems logical to read the file into a data frame using read.table(). Then you can change it into any format you want: > Listname <- read.table(text="key value Key1 1 Key2 2 Key3 3", header=TRUE) > Listname

Re: [R] Contacting Delphi ??

2012-10-10 Thread Rui Barradas
Unsuccessfull? Why unsuccessfull? Have you noticed the call to rpois? It beats runif by potential infinity. Rui Barradas Em 10-10-2012 18:30, David Winsemius escreveu: On Oct 10, 2012, at 9:59 AM, Christian Hoffmann wrote: What does the sudden appearance of "Contacting Delphi ..the oracle

[R] interaction testing in GAM

2012-10-10 Thread anna freni sterrantino
Hello! I'd like to know if it is correct to test with anova two models specified like this: m1=y~x1+s(x2,by=x3),family="poisson" m0=y~x1+s(x2),family="poisson" anova(m1,m0) Cheers Anna   Anna Freni Sterrantino Department of Statistics University of Bologna, Italy via Belle Arti 41, 40124 BO.

Re: [R] RMySQL install on windows

2012-10-10 Thread Greg Snow
I finally was able to compile/load it under windows 7. I had similar problems to what you show below. I set the MYSQL_HOME environmental variable through windows (start button > control panel > System and Security > system > Advanced System Settings > Environmental variables). I had to set it to

Re: [R] Contacting Delphi ??

2012-10-10 Thread Jeff Newmiller
Its suddenness is only as sudden as your input that triggered it. Its oblique nature is a reflection of the nature of your input. http://lmgtfy.com/?q=R+contacting+Delphi+oracle+is+unavailable --- Jeff Newmiller

Re: [R] reading in a (very simple) list from a file

2012-10-10 Thread arun
HI, By modifying the earlier solution using sapply() set.seed(1)  dat1<-data.frame(keys=paste0("key",1:5),value=sample(1:15,5,replace=TRUE)) list2<-sapply(split(dat1,dat1$keys),`[`,2)  names(list2)<-dat1[,1] list2$key2 #[1] 6 A.K. - Original Message - From: VA Smith To: r-help@r-project

Re: [R] append for .Rdata?

2012-10-10 Thread J Toll
On Tue, Oct 9, 2012 at 10:35 AM, Jessica Streicher wrote: > Can i somehow append objects to an .Rdata file? > > I didn't see an option for it in the save() method. > > dump() won't work since i have s4 objects in there. I'm not sure I completely understand the issues you're trying to work around

Re: [R] reading in a (very simple) list from a file

2012-10-10 Thread arun
Hi, May be this: set.seed(1)  dat1<-data.frame(keys=paste0("key",1:5),value=sample(1:15,5,replace=TRUE)) list1<-lapply(split(dat1,dat1$keys),`[`,2) list1$key2 #  value #2 6 A.K. - Original Message - From: VA Smith To: r-help@r-project.org Cc: Sent: Wednesday, October 10, 2012 1:29

Re: [R] Contacting Delphi ??

2012-10-10 Thread Rui Barradas
Hello, It's in the source code for `?`, file src/library/utils/R/question.R, lines 32 to 35. 32cat("Contacting Delphi...") 33 flush.console() 34Sys.sleep(2+rpois(1,2)) 35cat("the oracle is unavailable.\nWe apologize for any inconvenience.\n") 36return(in

[R] reading in a (very simple) list from a file

2012-10-10 Thread VA Smith
Apologies - I feel this is a very simple thing to do yet I am failing massively. I keep finding information about how to do much more complicated things (usually on this mailing list!), which then fail when I try to apply it to my simple task. Anyway, all I want to do is read in a series of key-va

Re: [R] Contacting Delphi ??

2012-10-10 Thread David Winsemius
On Oct 10, 2012, at 9:59 AM, Christian Hoffmann wrote: > What does the sudden appearance of "Contacting Delphi ..the oracle is > unavailable. > We apologize for any inconvenience." mean? A bug? It appears at plotting. In this instance an unsuccessful attempt at humor: http://markmail.org/s

Re: [R] contour plot help

2012-10-10 Thread David Winsemius
On Oct 10, 2012, at 7:58 AM, namrata mohapatra wrote: > Hello > > I am interested to plot a contour plot using the colour : rainbow , however I > want to reverse the order of the colour ( such that red represents max value > and blue min value) and also remove the lines in white in the plot .

Re: [R] Summary using by() returns character arrays in a list

2012-10-10 Thread arun
HI, May be this helps you: Using the dataset iris: by.list<-by(iris, iris$Species, summary) dat1<-do.call(rbind,lapply(by.list,function(x) gsub(".*\\:","",x))) row.names(dat1)<-paste(rep(unlist(dimnames(by.list),use.names=F),each=6),unlist(lapply(lapply(by.list,`[`,1:6),function(x) gsub("\\:.*"

[R] Exporting summary plm results to latex

2012-10-10 Thread Sebastian Barfort
Dear all, I am trying to export my fixed effect results to Latex. I am using the plm package with the summary function. However, it does not look like apsrtable, stargazer, or any other package can accompany using the plm package. I am interested in a classic table with the coefficient in one

[R] history() does not work? 2

2012-10-10 Thread Christian Hoffmann
Hi, history() gives Error in savehistory(file) : no history available to save although I can scroll throu history with C^uparrow an C^downarrow. How can I make history() work and/or show the current history in a *file*, so that I can choose from previous commands? The notion of having re

Re: [R] history() does not work?

2012-10-10 Thread Steve Lianoglou
On Wed, Oct 10, 2012 at 12:52 PM, Christian Hoffmann wrote: > > Am 10.10.12 18:17, schrieb Steve Lianoglou: > >> Hi, >> >> On Wed, Oct 10, 2012 at 12:03 PM, Christian Hoffmann >> wrote: >>> >>> Hi, >>> history() >>> >>> gives Error in savehistory(file) : no history available to save >>> >>>

[R] glmmPQL and spatial correlation

2012-10-10 Thread Luis Huckstadt
Hi all, I'm running into some computer issues when trying to run a binomial model for spatially correlated data using glmmPQL and was wondering if anyone could help me out. My whole dataset consists of about 300,000 points for which I have a suite of environmental variables (I'm trying to come up

[R] Contacting Delphi ??

2012-10-10 Thread Christian Hoffmann
What does the sudden appearance of "Contacting Delphi ..the oracle is unavailable. We apologize for any inconvenience." mean? A bug? It appears at plotting. Thanks Christian -- Christian W. Hoffmann, CH - 8915 Hausen am Albis, Switzerland Rigiblickstrasse 15 b, Tel.+41-44-7640853 c-w.hoffm

Re: [R] own function: computing time

2012-10-10 Thread William Dunlap
No, the desired points are not a subset of the convex hull. E.g., x=c(0,1:5), y=c(0,1/(1:5)). Bill Dunlap Spotfire, TIBCO Software wdunlap tibco.com > -Original Message- > From: William Dunlap > Sent: Wednesday, October 10, 2012 9:46 AM > To: 'tonja.krue...@web.de'; r-help@r-project.org

Re: [R] own function: computing time

2012-10-10 Thread William Dunlap
Are the points you are looking for (those data points with no other data points above or to the right of them) a subset of the convex hull of the data points? If so, chull(x,y) can quickly give you the points on the convex hull (typically a fairly small number) and you can look through them for th

Re: [R] a merge() problem

2012-10-10 Thread Sam Steingold
> * Prof Brian Ripley [2012-10-08 06:37:07 +0100]: > > On 08/10/2012 02:57, Peter Ehlers wrote: >> On 2012-10-07 14:44, Sam Steingold wrote: * Peter Ehlers [2012-10-07 10:03:42 -0700]: On 2012-10-07 08:34, Sam Steingold wrote: > I know it does not look very good - using the sam

Re: [R] Connect R and Lyx in UBUNTU

2012-10-10 Thread Yihui Xie
What do you mean by "connect"? And how did you succeed under Windows? i.e. what is your expectation? Regards, Yihui -- Yihui Xie Phone: 515-294-2465 Web: http://yihui.name Department of Statistics, Iowa State University 2215 Snedecor Hall, Ames, IA On Wed, Oct 10, 2012 at 4:57 AM, ATANU wrote:

Re: [R] synthetic distribution built upon set of discrete values

2012-10-10 Thread Greg Snow
In addition to Bert's answer. If the 0 and/or 100 are hard boundaries (you know that values cannot be outside those values) and you have data points near one or both of the bounds, then the functions in the logspline package may be of use. On Tue, Oct 9, 2012 at 2:38 PM, Bert Gunter wrote: > ??

Re: [R] history() does not work?

2012-10-10 Thread Steve Lianoglou
Hi, On Wed, Oct 10, 2012 at 12:03 PM, Christian Hoffmann wrote: > Hi, > >> history() > > gives Error in savehistory(file) : no history available to save > > although I can scroll throu history with C^uparrow an C^downarrow. > > How can I make history() work and/or show the current history in a fi

[R] history() does not work?

2012-10-10 Thread Christian Hoffmann
Hi, > history() gives Error in savehistory(file) : no history available to save although I can scroll throu history with C^uparrow an C^downarrow. How can I make history() work and/or show the current history in a file, so that I can choose from previous commands? The web did not throw up a

Re: [R] plot separate groups with plotmeans()

2012-10-10 Thread Bart Toonen
Hi, This answer may be a bit overdue, but I had a similar problem and it took me quite some time to find a good solution. When searching with google, your post kept showing up, so I post this to help others. The command that worked for me was interaction.plot. Example: x <- data.frame(Score=rno

[R] histogram & corelation plots

2012-10-10 Thread sagarnikam123
i want to show histogram in innovative way (good if interactive) in html report..is threre any tutorials/hint also there is any package to show correalation plot -- View this message in context: http://r.789695.n4.nabble.com/histogram-corelation-plots-tp4645726.html Sent from the R help mailing

Re: [R] se's and CI's for fitted lines in multivariate regression analysis

2012-10-10 Thread Rui Barradas
Hello, Your model is equivalent to the model below. As for standard errors, try predict.lm with the appropriate argument. ?predict.lm model <- lm(decrease ~ rowpos + colpos*treatment, data = OrchardSprays) predict(model, se.fit = TRUE, interval = "confidence") Hope this helps, Rui Barradas E

Re: [R] Summary using by() returns character arrays in a list

2012-10-10 Thread Rui Barradas
Hello, If 'by' is giving you trouble, why not 'aggregate'? agg.df <- aggregate(iris, list(iris$Species), FUN = summary) str(agg.df) Hope this helps, Rui Barradas Em 10-10-2012 15:02, Alex van der Spek escreveu: Thank you Petr, Try this str(by(iris, iris$Species, summary)) and you will see

Re: [R] write.csv with append = TRUE

2012-10-10 Thread Duncan Murdoch
On 10/10/2012 10:56 AM, R. Michael Weylandt wrote: write.table() should append just fine. write.csv() doesn't append I'm told because there's no guarantee it matches the column headings of the original file. That said, write.table(..., append = TRUE, sep = ",") is a pretty good approximation. Y

Re: [R] se's and CI's for fitted lines in multivariate regression analysis

2012-10-10 Thread Sigrid
Hi Bert, I just looked at An Introduction to R - and I do apologize if my questions are trivial. I see that they use predict as a function in lm, but I'm not sure how to incorporate it into a command. Thank you, S -- View this message in context: http://r.789695.n4.nabble.com/se-s-and-CI-s-f

Re: [R] Generating random geographical coordinates

2012-10-10 Thread R. Michael Weylandt
On Wed, Oct 10, 2012 at 4:09 PM, Poizot Emmanuel wrote: > Le 10/10/2012 17:02, R. Michael Weylandt a écrit : >> >> Change 1 to some other number to get more points from runif() >> >> More generally, take a look at "An Introduction to R" and read most >> everything you can find on the topic of vect

Re: [R] Generating random geographical coordinates

2012-10-10 Thread Berend Hasselman
On 10-10-2012, at 16:15, Poizot Emmanuel wrote: > Dear all, > > I have two coordinates vectors, say X and Y of length n. > I want to generate for each couple of coordinates X1,Y1 X2,Y2 X3,Y3Xn,Yn > a random coordinate which is located in a square define as X +/- dx and Y +/- > dy. > I saw

Re: [R] Generating random geographical coordinates

2012-10-10 Thread R. Michael Weylandt
Change 1 to some other number to get more points from runif() More generally, take a look at "An Introduction to R" and read most everything you can find on the topic of vectorization. If you don't know how to get "An Introduction to R", try typing help.start() at your prompt and it should happen

Re: [R] write.csv with append = TRUE

2012-10-10 Thread R. Michael Weylandt
write.table() should append just fine. write.csv() doesn't append I'm told because there's no guarantee it matches the column headings of the original file. That said, write.table(..., append = TRUE, sep = ",") is a pretty good approximation. Cheers, Michael On Wed, Oct 10, 2012 at 3:44 PM, Alex

  1   2   >