[R] Multi-output regression

2011-09-10 Thread zeec
Hi, I have a question regarding the modeling methodology of the following problem: * I have two data sets {X_i,y_i} {X_i,z_i}, i=1..N, where y_i = f(X_i) + i.i.d. Gaussian noise and z_i = g(X_i) + i.i.d. Gaussian noise * I apply bayesian linear regression to each of them and obtain p(y|X

[R] closed testing procedure

2011-09-10 Thread array chip
Hi, are the methods for multiple testing p value adjustment (Shaffer, Westfall, free) implemented in the function adjusted() in multcomp package so called closed testing procedure? what about those methods (holm, hochberg, hommel, BH, BY) implemented in the p.adjust() in the stats package?   Tha

Re: [R] how to get character from keyboard when run a file

2011-09-10 Thread Greg Snow
Are you running a script using 'source' or in batch mode? And do you want the script to get input from a user instead of from the script? If that is the case then you need to work a little harder because R assumes that if you are running in a non-interactive mode (source and batch mode) that i

Re: [R] DBS to R

2011-09-10 Thread Ben qant
The below is how scan() printed it in R...no obscured data here. I'm guessing that the 'XX's are place fillers because I read in only one ticker. The true 'source' is an SQL query. The other program isn't anything anyone would know about here. Thanks for any suggestions, ben On Sat, Sep 10, 2011

Re: [R] DBS to R

2011-09-10 Thread Sarah Goslee
I'm assuming you don't have access to the original program and source, but do you know what it is, and what version? (If you do have access to the original, then just export it in a different format.) The identity of the program would definitely help in reverse-engineering the format or finding a w

Re: [R] DBS to R

2011-09-10 Thread Ben qant
Recap: trying to get a 'dbs' extension file into R. I tried most of the functions in the 'forgein' package. Nothing there worked. Scan seemed to do the best. See below... I've tired using readBin() for hours and I haven't extracted the value I want, price. This example is a tiny sample dbs file t

Re: [R] ordering rows within CrossTable

2011-09-10 Thread Marc Schwartz
On Sep 10, 2011, at 4:03 PM, Luca Meyer wrote: > Hi, > > I am running the following -masked- code: > > set.seed(23) > city <- sample(c("C1","C2"),size=100,replace=T) > reason <- sample(c("R1","R2","R3","R4"),size=100,replace=T) > df <- data.frame(city,reason) > library(gmodels) > CrossTable(df$r

Re: [R] Code for The R Book

2011-09-10 Thread John C Frain
Peter, Paul Why not complain to Wiley. Uwe is 100% correct. John On Saturday, 10 September 2011, Uwe Ligges wrote: > > > On 08.09.2011 10:56, Peter Raundal wrote: >> >> Hi Paul, >> >> I'm struggling to find the R code as well on the books website. Did you >> managed to get it and are you willi

[R] ordering rows within CrossTable

2011-09-10 Thread Luca Meyer
Hi, I am running the following -masked- code: set.seed(23) city <- sample(c("C1","C2"),size=100,replace=T) reason <- sample(c("R1","R2","R3","R4"),size=100,replace=T) df <- data.frame(city,reason) library(gmodels) CrossTable(df$reason,df$city,prop.r=F,prop.c=F,prop.t=F,prop.chisq=F) And I get th

Re: [R] Accessing graphs via url

2011-09-10 Thread Dennis Murphy
Hi: Start here: http://www.quantmod.com/examples/ HTH, Dennis On Fri, Sep 9, 2011 at 10:44 PM, veepsirtt wrote: > Hi > I want to bring the graph from the site www.avasaram.com > for the symbol GOOG > here is the link > > http://www.avasaram.com/historyAction!addPatterns.action?model.ticker=GOOG

[R] npreg: plotting out of sample, extremely large bandwidths

2011-09-10 Thread Bert Lloyd
Hello r-help, I am using the excellent np package to conduct a nonparametric kernel regression and am having some trouble plotting the results. I have 2 covariates, x1 and x2, and a continuous outcome variable y. I am conducting a nonparametric regression of y on x1 and x2. The one somewhat unusu

Re: [R] Error : subscript out of bounds

2011-09-10 Thread David Winsemius
On Sep 10, 2011, at 4:33 PM, avinash barnwal wrote: Hi , I appreciate ur help, i do my homework carefully but here the problem is i m declaring a m as null matrix Yes. That IS the problem. to be filled by columns not by rows and afterwards i m storing the values directly into the matrix

Re: [R] Error : subscript out of bounds

2011-09-10 Thread avinash barnwal
Hi , I appreciate ur help, i do my homework carefully but here the problem is i m declaring a m as null matrix to be filled by columns not by rows and afterwards i m storing the values directly into the matrix and i had gone through same above mentioned error in previous r-help queries. I am not

Re: [R] Error : subscript out of bounds

2011-09-10 Thread jim holtman
?debug ?browser Read up on debugging your script. On Saturday, September 10, 2011, avinash barnwal < avinashbarnwal...@gmail.com> wrote: > m<-matrix(byrow=FALSE) > t<-as.list(na.exclude(x))_ > j<-0 > o<-0 > for(i in 1:998) > { > d<- 5*(i-1)+3 > if(t[[d]][[1]]>80) > { > j<-j+1 > e[j]<-d > l<-lengt

Re: [R] Error : subscript out of bounds

2011-09-10 Thread R. Michael Weylandt
Take a look at "m". Your error should be quite clear More generally: we don't do homework help on this list. Michael Weylandt On Sat, Sep 10, 2011 at 12:59 PM, avinash barnwal < avinashbarnwal...@gmail.com> wrote: > m<-matrix(byrow=FALSE) > t<-as.list(na.exclude(x)) > j<-0 > o<-0 > for(i in 1:9

[R] Error : subscript out of bounds

2011-09-10 Thread avinash barnwal
m<-matrix(byrow=FALSE) t<-as.list(na.exclude(x)) j<-0 o<-0 for(i in 1:998) { d<- 5*(i-1)+3 if(t[[d]][[1]]>80) { j<-j+1 e[j]<-d l<-length(t[[d]]) u<-t[[d]] price_rand<-t[[d-1]] n<-0 for(k in 1:l) { if((u[k]>49)&&(u[k]<51)) { n<-n+1 m[n,j]<-price_rand[k] } } } } I am getting error in assigning

Re: [R] DBS to R

2011-09-10 Thread R. Michael Weylandt
No experience with dbs files, but the foreign package might be able to help you out. Michael On Sat, Sep 10, 2011 at 2:44 PM, Ben qant wrote: > Hello, > > I have a bunch of data files all with "dbs" file extensions. They are > generated via a SQL query from another program and source. Does anyo

[R] DBS to R

2011-09-10 Thread Ben qant
Hello, I have a bunch of data files all with "dbs" file extensions. They are generated via a SQL query from another program and source. Does anyone know (or have ideas) how to get the data from a dbs file type into R (or into some other format that can imported to R)? I've searched online for 4 ho

Re: [R] simple t test for multiple colimns

2011-09-10 Thread Daniel Malter
It will never hurt your chances to get an answer if you tell what the code in the different language actually does. Is it doing t-est for all of the variables and the group indicators are in a separate column? If, imagine you have the following data: id<-rep(c(1,2),each=100) x<-c(rnorm(100,0,1),rn

Re: [R] General help - online statistics courses?

2011-09-10 Thread Kjetil Halvorsen
On Thu, Sep 8, 2011 at 5:58 PM, Rolf Turner wrote: > On 09/09/11 06:46, Spencer Graves wrote: > > >> >> I'd rather ask a stupid question than make a stupid mistake. > > > > Fortune? > Its already there, more or less: library(fortunes) fortune("spencer") kjetil >    cheers, > >        Rolf Tur

Re: [R] prcomp: results with reversed sign in output?

2011-09-10 Thread Vokey, John
In addition, many PCA packages follow the convention that if the majority of weights are negative for that component, reverse the sign. On 2011-09-10, at 4:00 AM, r-help-requ...@r-project.org wrote: > The point is that a principal component vector is a solution, > say V, of a matrix equation A%*

Re: [R] Very slow using S4 classes

2011-09-10 Thread Martin Maechler
> Martin Morgan > on Sat, 10 Sep 2011 10:18:11 -0700 writes: > On 09/10/2011 08:08 AM, André Rossi wrote: >> Hi everybody! >> >> I'm creating an object of a S4 class that has two slots: ListExamples, >> which is a list, and idx, which is an integer (as the code be

Re: [R] how to get character from keyboard when run a file

2011-09-10 Thread R. Michael Weylandt
scan(, what=character()) Michael Weylandt On Sat, Sep 10, 2011 at 10:16 AM, bebece wrote: > when run a R file, i can get numeric by scan(), but scan() can't read > character or i make a wrong. > so how can i get character from keyboard when i run a R file? > > -- > View this message in context

Re: [R] Variable name with a '-'

2011-09-10 Thread peter dalgaard
On Sep 10, 2011, at 16:33 , Adam Moore wrote: > My statistics teacher wrote a textbook and included .rdata files for the > problem set data. I downloaded these and opened one of them up. Then I used > the ls() command, and it gave me this: > >> ls() > [1] "ex02-12" > > However, when i try to us

Re: [R] dgamma in jags within r

2011-09-10 Thread Martin Maechler
> Marc Girondot > on Sat, 10 Sep 2011 11:43:20 +0200 writes: > Marc Girondot > on Sat, 10 Sep 2011 11:43:20 +0200 writes: > I define priors in jags within r using a gamma distribution. I would > like to control the shape but I have problems. Any help will be usef

[R] how to get character from keyboard when run a file

2011-09-10 Thread bebece
when run a R file, i can get numeric by scan(), but scan() can't read character or i make a wrong. so how can i get character from keyboard when i run a R file? -- View this message in context: http://r.789695.n4.nabble.com/how-to-get-character-from-keyboard-when-run-a-file-tp3803944p3803944.html

Re: [R] General help - online statistics courses?

2011-09-10 Thread kydaviddoyle
Hi, Another option is http://www.khanacademy.org/ Very interesting site. The guy started out helping out some of his relatives with school work and then started to do them as videos. His relatives liked them even better being they could pause/reply them if they didn't understand. They told fri

[R] Variable name with a '-'

2011-09-10 Thread Adam Moore
My statistics teacher wrote a textbook and included .rdata files for the problem set data. I downloaded these and opened one of them up. Then I used the ls() command, and it gave me this: > ls() [1] "ex02-12" However, when i try to use the object, it interprets it as ex02 - 12, and I get the foll

Re: [R] Very slow using S4 classes

2011-09-10 Thread Martin Morgan
On 09/10/2011 08:08 AM, André Rossi wrote: Hi everybody! I'm creating an object of a S4 class that has two slots: ListExamples, which is a list, and idx, which is an integer (as the code below). Then, I read a data.frame file with 1 (ten thousands) of lines and 10 columns, do some pre-proce

Re: [R] Very slow using S4 classes

2011-09-10 Thread Uwe Ligges
On 10.09.2011 17:08, André Rossi wrote: Hi everybody! I'm creating an object of a S4 class that has two slots: ListExamples, which is a list, and idx, which is an integer (as the code below). Then, I read a data.frame file with 1 (ten thousands) of lines and 10 columns, do some pre-proces

Re: [R] Save Workspace Image? problem at start up

2011-09-10 Thread Uwe Ligges
So you have added a q() or quit() in one of the startup scripts? Uwe Ligges On 07.09.2011 11:50, Chandrasekhar Rudrappa wrote: Dear Sir, I am using R 2.13.1. Since yesterday when I start the Rgui the message "Save Workspace Image?" shows up immediately on opening. Even if I cancel, it does

Re: [R] error in knn: too many ties in knn

2011-09-10 Thread Uwe Ligges
On 08.09.2011 14:19, Максим Иванов wrote: Hello. I found the behavior of knn( http://stat.ethz.ch/R-manual/R-devel/library/class/html/knn.html) function looking very strange. Consider the toy example. library(class) train<- matrix(nrow=5000,ncol=2,data=rnorm(1,0,1)) test<- matrix(nrow=10,nc

Re: [R] Code for The R Book

2011-09-10 Thread Uwe Ligges
On 08.09.2011 10:56, Peter Raundal wrote: Hi Paul, I'm struggling to find the R code as well on the books website. Did you managed to get it and are you willing to share it with me? So time to ask the author who probably does not necessarily read this list. Or just buy another book that is

Re: [R] error building package: packaging into .tar.gz failed

2011-09-10 Thread Uwe Ligges
On 07.09.2011 19:02, Juliet Hannah wrote: To follow up (because I received a few emails off-list), things work now. I'm not sure what I did differently. In case it is helpful: I reinstalled R tools, and edited by path so that C:\Rtools\bin;C:\Rtools\perl\bin;C:\Rtools\MinGW\bin;C:\Program fi

[R] Very slow using S4 classes

2011-09-10 Thread André Rossi
Hi everybody! I'm creating an object of a S4 class that has two slots: ListExamples, which is a list, and idx, which is an integer (as the code below). Then, I read a data.frame file with 1 (ten thousands) of lines and 10 columns, do some pre-processing and, basically, I store each line as an

Re: [R] Editing the variables attributes section in the netCDF header of netCDF files created using the package ncdf.

2011-09-10 Thread David William Pierce
On Fri, Sep 9, 2011 at 3:18 PM, RMaidment wrote: > Another question which you might be able to answer is, is it possible to set > the order of the variables attributes in which they appear? [...] Hi Ross, The ncdf package makes the units and longname attributes automatically, so those will alwa

Re: [R] reshape data from long to wide format

2011-09-10 Thread Gabor Grothendieck
On Fri, Sep 9, 2011 at 7:28 AM, maxbre wrote: > This is my reproducible example: > > example<-structure(list(SENSOR = structure(1:6, .Label = c("A", "B", "C", > "D", "E", "F"), class = "factor"), VALUE = c(270, 292.5, 0, 45, > 247.5, 315), DATE = structure(1:6, .Label = c(" 01/01/2010 1", > " 01/0

Re: [R] reshape data from long to wide format

2011-09-10 Thread maxbre
thank you to both you, this is exactly the solution I was looking for max -- View this message in context: http://r.789695.n4.nabble.com/reshape-data-from-long-to-wide-format-tp3801381p3803759.html Sent from the R help mailing list archive at Nabble.com.

Re: [R] R-help Digest, Vol 103, Issue 9

2011-09-10 Thread mihalicza . peter
Szeptember 12-től 26-ig irodán kívül vagyok, és az emailjeimet nem érem el. Sürgős esetben kérem forduljon Kárpáti Edithez (karpati.e...@gyemszi.hu). Üdvözlettel, Mihalicza Péter I will be out of the office from 12 till 26 September with no access to my emails. In urgent cases please contact

Re: [R] Get previous R digest volumes/issues

2011-09-10 Thread Uwe Ligges
On 10.09.2011 12:57, (Ted Harding) wrote: On 10-Sep-11 10:19:54, soeren.vo...@uzh.ch wrote: Hello, (how) can I download/re-retrieve/order previous R-** digest volumes/issues to my mailbox for local browsing? Thank you, Sören Taking as example the R-help digest: https://stat.ethz.ch/piper

Re: [R] Rgl and plotmath symbols (via sprites): a trial

2011-09-10 Thread Duncan Murdoch
On 11-09-09 2:03 PM, Duncan Murdoch wrote: On 09/09/2011 12:19 PM, Marius Hofert wrote: Dear Duncan, attached my next trial. It seems to work :-) Unfortunately, not the printing to pdf. Not getting a high-quality pdf is certainly one of the major drawbacks. If anybody knows how to fix that, p

Re: [R] Get previous R digest volumes/issues

2011-09-10 Thread Ted Harding
On 10-Sep-11 10:19:54, soeren.vo...@uzh.ch wrote: > Hello, (how) can I download/re-retrieve/order previous R-** digest > volumes/issues to my mailbox for local browsing? Thank you, Sören Taking as example the R-help digest: https://stat.ethz.ch/pipermail/r-help/ for each month listed, you will

[R] Get previous R digest volumes/issues

2011-09-10 Thread soeren . vogel
Hello, (how) can I download/re-retrieve/order previous R-** digest volumes/issues to my mailbox for local browsing? Thank you, Sören __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http:

[R] dgamma in jags within r

2011-09-10 Thread Marc Girondot
I define priors in jags within r using a gamma distribution. I would like to control the shape but I have problems. Any help will be usefull. From help of dgamma ___ The Gamma distribution with parameters shape = a and scale = s has density f(x)= 1/(s^a Gamma(a)) x^(a-1) e^-(x/s)

[R] Mutiple vertical reference lines with a dotplot

2011-09-10 Thread markm0705
Dear R helpers I would like to plot vertical reference lines for yearly averages the variable being plotted. In this case there are three years 2002, 2009 and 2010 for a variable Resc_Gt. I've attempted to get the 2002 mean plotting (2.27) using the 'abline' command but no joy. Ideally I woul

Re: [R] Different results with arima in R 2.12.2 and R 2.11.1

2011-09-10 Thread Uwe Ligges
On 10.09.2011 01:54, Stephen Sefick wrote: Those are not "very" different. I would suspect they are within computer rounding error. Only numbers of log2 can be represented exactly, could that be it? Did you compiled the two versions differently, or could something have changed between 2.11.1