Re: [R] strsplit question

2011-10-12 Thread Joshua Wiley
unlist(strsplit(Block[1:5], -.+$)) if you are going to want the other pieces later, the most efficient way depends on the assumptions you can make about your data. If there are always two elements from the split: matrix(unlist(strsplit(Block[1:5], -)), ncol = 2, byrow = TRUE) ## or

Re: [R] strsplit question

2011-10-12 Thread Jeff Newmiller
sapply(strsplit(Block[1:5],-), function (x) {x[1]}) comes to mind... --- Jeff Newmiller The . . Go Live... DCN:jdnew...@dcn.davis.ca.us Basics: ##.#. ##.#. Live Go... Live: OO#.. Dead: OO#.. Playing Research Engineer

[R] Using a for loop

2011-10-12 Thread JohnnyJames
I am using the normtemp data set in UsingR package. I want to store all normtemp[,1] from 1 to 65 as x and normtemp[,1] from 66to 130 as y. How do I write code to do that from here: data(normtemp) -- View this message in context:

[R] PMML for Cox Regression

2011-10-12 Thread Rick Pechter
Hi, We have created a cox regression model (coxph) and exported from R as PMML using the Rattle package. The resulting PMML is a standard Regression model in the format y=b1*x1+b2*x2+…….bn*xn: risk= (0.027968680616809*dropped_Calls)+(0.00777220409115466*helpdesk_Calls) The PMML is included at

[R] ploting pie charts centered at given points

2011-10-12 Thread JeffND
Hi to all, Suppose we have a group of points on the plane, at each point, we want to draw a pie chart centered at that point. I have found that pie.labels() function might be useful. This function allows us to achieve my goal, but before using the function, we have to use the function plot(),

Re: [R] strsplit question

2011-10-12 Thread Remko Duursma
Hi Erin, this is one way: Block - c(5600-5699,6100-6199,9700-9799,9400-9499,8300-8399) splBlock - strsplit(Block,-) sapply(splBlock, [, 1) greetings, Remko -- View this message in context: http://r.789695.n4.nabble.com/strsplit-question-tp3896847p3896850.html Sent from the R help

Re: [R] Using a for loop

2011-10-12 Thread Joshua Wiley
You should read the posting guide again (at all?). Among other things, it suggests you should show what you have tried doing that did not work (worst case it shows us you tried to figure it out on your own before asking other people to freely give their time to solve your problems, at best case,

Re: [R] strsplit question

2011-10-12 Thread David Winsemius
On Oct 12, 2011, at 1:20 AM, Erin Hodgess wrote: Dear R People: I have the following set of data Block[1:5] [1] 5600-5699 6100-6199 9700-9799 9400-9499 8300-8399 and I want to split at the - strsplit(Block[1:5],-) [[1]] [1] 5600 5699 [[2]] [1] 6100 6199 [[3]] [1] 9700 9799 [[4]] [1]

Re: [R] ploting pie charts centered at given points

2011-10-12 Thread Jim Lemon
On 10/12/2011 03:15 PM, JeffND wrote: Hi to all, Suppose we have a group of points on the plane, at each point, we want to draw a pie chart centered at that point. ... Hi Jeff, I think what you are looking for is the floating.pie function (plotrix). Jim

Re: [R] Creating the mean using algebra matrix

2011-10-12 Thread Timothy Bates
On Oct 12, 2011, at 12:04 AM, Rolf Turner wrote: On 12/10/11 08:31, Timothy Bates wrote: To do matrix multiplication: m x n, the Rows and columns of m must be equal to the columns and rows of n, respectively. No. The number of columns of m must equal the number of rows of n, that's all.

Re: [R] p adjustment on 4thcorner results

2011-10-12 Thread Christoph Molnar
Hi, ?p.adjust Christoph 2011/10/12 Cristina Ramalho cristina.rama...@grs.uwa.edu.au Hi all, This is probably a very simple question but I cannot figure out how to do it. I run the fourthcorner method with my data and would like to adjust the p values for multiple comparisons using Holm

[R] goofy class question

2011-10-12 Thread Erin Hodgess
Dear R People: Here is a really goofy question. I have some objects which have 2 classes: data.frame and ucr. Also, the classes will always be in that order. I have tried all sorts of things, but to no avail. listucrModels - function(envir=.GlobalEnv, ...) { objects - ls(envir=envir, ...)

Re: [R] goofy class question

2011-10-12 Thread Joshua Wiley
names(lapply(.GlobalEnv, function(x) inherits(x, ucr))) HTH, Josh On Wed, Oct 12, 2011 at 12:46 AM, Erin Hodgess erinm.hodg...@gmail.com wrote: Dear R People: Here is a really goofy question. I have some objects which have 2 classes: data.frame and ucr. Also, the classes will always be

Re: [R] Parallel processing for R loop

2011-10-12 Thread Eik Vettorazzi
Hi Sandeep, still missing an answer? Perhaps you cross check your post with the rules of the posting guide and find what is missing at all here. Anyway, depending on your OS, package multicore, snow/snowfall may fit your needs - but you have to re-formulate your loop using adequate multicore

Re: [R] goofy class question

2011-10-12 Thread Barry Rowlingson
On Wed, Oct 12, 2011 at 8:46 AM, Erin Hodgess erinm.hodg...@gmail.com wrote: Dear R People: Here is a really goofy question. I have some objects which have 2 classes: data.frame and ucr. Also, the classes will always be in that order. I have tried all sorts of things, but to no avail.

[R] Usng MCMCpack,error is \\\initial value in vmmin is not finite\\\

2011-10-12 Thread yiy83102
__ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.

Re: [R] Problem executing function

2011-10-12 Thread Divyam
Hi Michael, Thanks for the reply, but still the problem exists. When I list the objects and return the result, they get printed on the console but the objects do not get created. I am really baffled and clueless as to what the problem is. Divya -- View this message in context:

[R] plot probability density function (pdf)

2011-10-12 Thread pigpigmeow
I have 2 series of variables, I want to plot the probability density function of these 2 variabels (i.e. two curves in one graph), I just want to compare these two variable distribution. what should I do? can I use ggplot2 package? -- View this message in context:

Re: [R] Parallel processing for R loop

2011-10-12 Thread Paul Hiemstra
On 10/11/2011 12:13 PM, Sandeep Patil wrote: I have an R script that consists of a for loop that repeats a process for many different files. I want to process this parallely on machine with multiple cores, is there any package for it ? Thanks ...I mostly use the foreach package...

Re: [R] p adjustment on 4thcorner results

2011-10-12 Thread Christoph Molnar
Hi, the function fourthcorners() returns a list. You can access the pvalues shown in the summary by four1[[tabGProb]]. This is a matrix containing the pvalues in the summary. I did some copy and paste from the summary.4thcorners() function (which is called, by using summary on your four1

Re: [R] Problem executing function

2011-10-12 Thread Jim Holtman
show the work you did. what is happening, probably, is that you are returning the values, but not assigning them. Sent from my iPad On Oct 12, 2011, at 3:47, Divyam divyamural...@gmail.com wrote: Hi Michael, Thanks for the reply, but still the problem exists. When I list the objects and

Re: [R] Problem executing function

2011-10-12 Thread Petr PIKAL
Hi Michael, Thanks for the reply, but still the problem exists. When I list the objects and return the result, they get printed on the console but the objects do not get created. I am really baffled and clueless as to what the problem is. I bet you do not understand basic operations

[R] Append data to vector form a column of a dataframe

2011-10-12 Thread behave
Dear R-community, When doing this: test-data.frame(a=c(1,2,3)) rbind(test$a, 3) I expect something like: 1 2 3 2 but get: [,1] [,2] [,3] [1,]123 [2,]222 the same for: rbind(test[[a]], 2) or rbind(as.vector(test[[a]]), 2) or rbind(t(as.vector(test[[a]])),

Re: [R] Append data to vector form a column of a dataframe

2011-10-12 Thread Ivan Calandra
Hi Dom, This is because 3 is recycled. It is necessary because the number of columns have to be the same for every row. Try this instead: c(test$a, 2) ## you wrote 3 but meant 2 I guess HTH, Ivan Le 10/12/2011 10:47, behave a écrit : Dear R-community, When doing this:

Re: [R] extra digits added to data

2011-10-12 Thread Jim Holtman
what are you going to do with the data? If just for presentation, then keep as character. If you are going to compute on the data, then keep as numeric. Since you are using floating point, FAQ 7.31 reminds you that the data is kept as inputted to the best that can be done with 54 bits of

Re: [R] p adjustment on 4thcorner results

2011-10-12 Thread Cristina Ramalho
Hi Christoph, Thank you soo much, this is great! All the very best, Cristina On Wed, Oct 12, 2011 at 4:52 PM, Christoph Molnar christoph.mol...@googlemail.com wrote: Hi, the function fourthcorners() returns a list. You can access the pvalues shown in the summary by four1[[tabGProb]].

Re: [R] strsplit question

2011-10-12 Thread Gabor Grothendieck
On Wed, Oct 12, 2011 at 1:20 AM, Erin Hodgess erinm.hodg...@gmail.com wrote: Dear R People: I have the following set of data Block[1:5] [1] 5600-5699 6100-6199 9700-9799 9400-9499 8300-8399 and I want to split at the - strsplit(Block[1:5],-) [[1]] [1] 5600 5699 [[2]] [1] 6100 6199

Re: [R] Chi-Square test and survey results

2011-10-12 Thread Jean V Adams
gheine wrote on 10/11/2011 02:31:46 PM: An organization has asked me to comment on the validity of their recent all-employee survey. Survey responses, by geographic region, compared with the total number of employees in each region, were as follows: ByRegion All.Employees

Re: [R] R and Forex

2011-10-12 Thread Wolfgang Wu
The quantmod package might be a good start.  http://cran.r-project.org/web/packages/quantmod/index.html Regards,   Wolfgang Wu - Ursprüngliche Message - Von: Yves S. Garret yoursurrogate...@gmail.com An: r-help@r-project.org Cc: Gesendet: 2:29 Mittwoch, 12.Oktober 2011 Betreff:

[R] Problem with pmax and matrix to calculate row maxima

2011-10-12 Thread Wolfgang Wu
I am having the following problem. I want to calculate the maximum of each row in a matrix. If I pass in the matrix split up by each column then this is no problem and works great. However I don't know how many columns I have in advance. In the example below I have 3 columns, but the number of

Re: [R] Append data to vector form a column of a dataframe

2011-10-12 Thread Wolfgang Wu
In addition to Ivan, test$a is not a data.frame anymore but a numerical vector. class(test$a) [1] numeric test$a [1] 1 2 3 So adding a row to your data.frame would be  rbind(test, 2)   a 1 1 2 2 3 3 4 2   Wolfgang Wu - Ursprüngliche Message - Von: Ivan Calandra

Re: [R] Problem executing function

2011-10-12 Thread Divyam
ok. I tested it in two ways. I want to externalise my odbcConnection details dsn, uid, and pwd. Hence I created a csv file to have these information. Like I showed in the sample function initially, the order of the steps were 1) loading of the packages, 2) fetching the csv file, 3)assigning the

Re: [R] Problem executing function

2011-10-12 Thread Divyam
sorry correction in the first and result hence received code: ok. I tested it in two ways. I want to externalise my odbcConnection details dsn, uid, and pwd. Hence I created a csv file to have these information. Like I showed in the sample function initially, the order of the steps were 1)

Re: [R] plot probability density function (pdf)

2011-10-12 Thread jdospina
x=rnorm(100,1,0.8) # A series. y=rnorm(100,0,0.5) # Another series with different mean and variance. plot(density(x),ylim=c(0,1)) lines(density(y),col=red) Remember that density() is a nonparametric estimator. You should properly choose the bandwith. -- View this message in context:

[R] monotonic factors

2011-10-12 Thread Jeffrey Pollock
Hello all, I have an ordered factor that I would like to include in the linear predictor of a binomial glm, where the estimated coefficients are constrained to be monotonic. Does anyone know how to do this? I've tried using an ordered factor but this does not have the desired effect, an

[R] Cross posting (was Restricted Cubic Splines within survfit.cph)

2011-10-12 Thread Terry Therneau
You wrote: It may be best to either write to the package maintainer (me, as you did) or post to the group but not both. This is just a note that I disagree wrt my own packages: I go on vacation or trips, or have other projects so won't always answer Other folks on the list often have

Re: [R] Nonlinear regression aborting due to error

2011-10-12 Thread Jean V Adams
Dennis Fisher wrote on 10/11/2011 07:20:35 PM: Colleagues, I am fitting an Emax model using nls. The code is: START - list(EMAX=INITEMAX, EFFECT=INITEFFECT, C50=INITC50) CONTROL - list(maxiter=1000, warnOnly=T) #FORMULA - as.formula(YVAR ~ EMAX - EFFECT *

[R] labels in a boxplot

2011-10-12 Thread Francesco Sarracino
Dear R-listers, I have a little problem with a boxplot and I hope you can help me figuring it out. I'll try to make up some data to illustrate the issue. Sorry, if my procedures look naive, but these are my first steps in R. Any comments and/or suggestions are very welcome. let's create a vector

Re: [R] Cross posting (was Restricted Cubic Splines within survfit.cph)

2011-10-12 Thread Frank Harrell
Terry I would just add that if someone contacts the maintainer and does not follow the advice of the maintainer, still has problems, and posts a message to the list, then some time is wasted. Frank On 10/12/2011 07:38 AM, Terry Therneau wrote: You wrote: It may be best to either write to

Re: [R] Problem with pmax and matrix to calculate row maxima

2011-10-12 Thread Enrico Schumann
Hi Wolfgang, how about a loop? matRandom - matrix(runif(n=60), ncol=6) ## variant 1 system.time(test1 - pmax(matRandom[,1], matRandom[,2], matRandom[,3], matRandom[,4], matRandom[,5], matRandom[,6])) User System verstrichen 0.010.000.01 ##

Re: [R] stop()

2011-10-12 Thread Doran, Harold
Thank you, Greg. This indeed works well for this purpose. -Original Message- From: Greg Snow [mailto:greg.s...@imail.org] Sent: Tuesday, October 11, 2011 4:27 PM To: Doran, Harold; r-help@r-project.org Subject: RE: stop() Replace stop() with break to see if that does what you

[R] Applying function to only numeric variable (plyr package?)

2011-10-12 Thread Michael . Laviolette
My data frame consists of character variables, factors, and proportions, something like c1 - c(A, B, C, C) c2 - factor(c(1, 1, 2, 2), labels = c(Y,N)) x - c(0.5234, 0.6919, 0.2307, 0.1160) y - c(0.9251, 0.7616, 0.3624, 0.4462) df - data.frame(c1, c2, x, y) pct - function(x) round(100*x, 1) I

Re: [R] Problem executing function

2011-10-12 Thread Duncan Murdoch
Comments inline below. On 12/10/2011 7:47 AM, Divyam wrote: sorry correction in the first and result hence received code: ok. I tested it in two ways. I want to externalise my odbcConnection details dsn, uid, and pwd. Hence I created a csv file to have these information. Like I showed in the

Re: [R] Chi-Square test and survey results

2011-10-12 Thread Jan van der Laan
George, Perhaps the site of the RISQ project (Representativity indicators for Survey Quality) might be of use: http://www.risq-project.eu/ . They also provide R-code to calculate their indicators. HTH, Jan Quoting ghe...@mathnmaps.com: An organization has asked me to comment on the

Re: [R] Applying function to only numeric variable (plyr package?)

2011-10-12 Thread Christoph Molnar
Hi, if the rows in your data.frame are numeric, this solution will work. (numeric.index - unlist(lapply(df, is.numeric))) df[, numeric.index] - apply(df[,numeric.index], 2, pct) This does not work for the example you gave, unless you coerce the columns with the your numerics to numeric: c1 -

Re: [R] Mean or mode imputation fro missing values

2011-10-12 Thread Petr PIKAL
Hi Yes thank you Gu… I am just trying to do this as a rough step and will try other imputation methods which are more appropriate later. I am just learning R, and was trying to do the for loop and f-statement by hand but something is going wrong… This is what I have until now:

Re: [R] Applying function to only numeric variable (plyr package?)

2011-10-12 Thread Jan van der Laan
plyr isn't necessary in this case. You can use the following: cols - sapply(df, is.numeric) df[, cols] - pct(df[,cols]) round (and therefore pct) accepts a data.frame and returns a data.frame with the same dimensions. If that hadn't been the case colwise might have been of help:

[R] Odp: labels in a boxplot

2011-10-12 Thread Petr PIKAL
Dear R-listers, I have a little problem with a boxplot and I hope you can help me figuring it out. I'll try to make up some data to illustrate the issue. Sorry, if my procedures look naive, but these are my first steps in R. Any comments and/or suggestions are very welcome. let's

Re: [R] labels in a boxplot

2011-10-12 Thread Francesco Sarracino
Thanks a lot Andrés. It was easier than I expected. f. 2011/10/12 Andrés Aragón armand...@gmail.com Francesco, Try cex.axis=0.6 Regards, Andrés AM 2011/10/12, Francesco Sarracino f.sarrac...@gmail.com: Dear R-listers, I have a little problem with a boxplot and I hope you can help

Re: [R] labels in a boxplot

2011-10-12 Thread Francesco Sarracino
Hi Petr, thanks a lot for your reply. Unfortunately, your suggestion does not work for me. I even tried larger boxes such as 15,15 , but the result does not change. Is there some setting that I am missing? However, once more thanks a lot for your help. f. On 12 October 2011 15:58, Petr PIKAL

Re: [R] Labels in ICLUST

2011-10-12 Thread William Revelle
On Oct 11, 2011, at 2:28 AM, Steve Powell wrote: Dear all, I can't get the labels slot in ICLUST to accept a character vector. library(psych) test.data - Harman74.cor$cov ic.out - ICLUST(test.data,nclusters =4,labels=letters[1:ncol(test.data)]) ## Error in !labels : invalid argument type

Re: [R] plot probability density function (pdf)

2011-10-12 Thread pigpigmeow
however, if i have an excel file, but there have 6 variables, a,b,c,d,e,f. how to plot the probability density function of a and d in one graph, b and e in another graph? -- View this message in context:

[R] Generelized Negative Binomial model in R

2011-10-12 Thread Akram Khaleghei Ghosheh balagh
Hello; Does anybody knows that R have a function for Generelized Negative Binomial model, something like gnbreg in STATA where dispersion parameter itself is a function of covaraites ? Thanks; [[alternative HTML version deleted]] __

Re: [R] how to calculate the statistics of a yearly window with a rolling step as 1 day?

2011-10-12 Thread ecoc
This doesn't work becaues the rollappy is non-overlapping. My rolling step is 1-day and rolling window is 1-year, so there is 364 days overlapping. -- View this message in context:

Re: [R] R and Forex

2011-10-12 Thread Yves S. Garret
No, that's not what I meant. I was curious if anyone has ever done this before and how well it worked. Any tips for a novice? On Wed, Oct 12, 2011 at 12:19 AM, Liviu Andronic landronim...@gmail.comwrote: On Wed, Oct 12, 2011 at 3:29 AM, Yves S. Garret yoursurrogate...@gmail.com wrote: Hi

[R] R loop within a loop

2011-10-12 Thread saliluna
Hi, I'm working on a loop function for a large dataset which contains 1000 different groups. I would like to reconstruct the order of events within each group by using a loop function in R. (Currently the order of events are based on the ascending order of prev_event within the group) A demo

[R] Error in Rcpp/inline (Windows XP)

2011-10-12 Thread allenhahaha
Hi, everyone, I am just trying to use Rcpp in my computer, and I would like to try a simple example from website, but R keeps reporting me error. I am using Windows XP, and has installed Rtools and GSI. Here is the response: src = ' + Rcpp::NumericVector xa(a); + Rcpp::NumericVector xb(b); +

Re: [R] Problem executing function

2011-10-12 Thread Divyam
Thanks a lot Duncan! It works fine now. Divya -- View this message in context: http://r.789695.n4.nabble.com/Problem-executing-function-tp3894359p3898054.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org mailing

Re: [R] R and Forex

2011-10-12 Thread Yves S. Garret
When can I send stuff to the mailing list without having moderator approval? Is that possible? On Wed, Oct 12, 2011 at 12:19 AM, Liviu Andronic landronim...@gmail.comwrote: On Wed, Oct 12, 2011 at 3:29 AM, Yves S. Garret yoursurrogate...@gmail.com wrote: Hi all, I recently started

Re: [R] labels in a boxplot

2011-10-12 Thread Andrés Aragón
Francesco, Try cex.axis=0.6 Regards, Andrés AM 2011/10/12, Francesco Sarracino f.sarrac...@gmail.com: Dear R-listers, I have a little problem with a boxplot and I hope you can help me figuring it out. I'll try to make up some data to illustrate the issue. Sorry, if my procedures look

Re: [R] R loop within a loop

2011-10-12 Thread saliluna
I've also tried to make the function work for one particular group, then apply the same function to the whole data frame with all groups using by() or lapply() as follow. But I'm still receiving error messages. Could someone please explain what is happening here? dfdfdf - function(localdata){

[R] CVbinary - Help

2011-10-12 Thread anamiguita
Hey, I need some help. I want to obtain a cross validation for a regression model (binary response) but I got an error with CVbinary. Well I did this: fit - lm(resp ~ PC1 + PC2 + PC3 + PC4 + PC5 + PC6 + PC7 + PC8 +

Re: [R] plot probability density function (pdf)

2011-10-12 Thread pigpigmeow
i want to plot probability density function,predictvalue has missing value and observevalue has not missing value, I tried .. attach(test) #test is the name of the data file names(test) plot(density(predictvalue,na.rm=TRUE)) lines(density(observevalue)) is it correct? -- View this message in

Re: [R] plot probability density function (pdf)

2011-10-12 Thread R. Michael Weylandt
Assuming you mean you want them on the same device: layout(1:2) plot(density(a)) lines(density(d),col=2) plot(density(b)) lines(density(e),col=2) Getting your data into R is more of a challenge, but if you want my unsolicited advice, you can do far worse than saving as CSV and using read.csv()

[R] loop function within a loop

2011-10-12 Thread Sally Zhen
Hi all, I'm working on a loop function for a large dataset which contains 1000 different groups. I would like to reconstruct the order of events within each group by using a loop function in R. (Currently the order of events are based on the ascending order of prev_event within the group) A

Re: [R] plot probability density function (pdf)

2011-10-12 Thread jdospina
x11() plot(density(a)) lines(density(d)) x11() plot(density(b)) lines(density(e)) 2011/10/12 pigpigmeow [via R] ml-node+s789695n3898183...@n4.nabble.com: however, if i have an excel file, but there have 6 variables, a,b,c,d,e,f. how to plot the probability density function of a and d in one

Re: [R] Error in Rcpp/inline (Windows XP)

2011-10-12 Thread Prof Brian Ripley
Did you read the posting guide (it says R-help is not for questions about compiled code)? Or the rw-FAQ (it warned you that lots of contributed packages did not take account of spaces in path names, so suggested you install R in a path without one)? This is a bug in one of the packages you

Re: [R] Minimum cutsets

2011-10-12 Thread malhomidi
Hi Gabor, I'm looking for minimum cutsets in the igraph manual but I didn't find the functions you mentioned above. Also, how can I see their source code. Thanks, Mohammed -- View this message in context: http://r.789695.n4.nabble.com/Minimum-cutsets-tp885346p3898347.html Sent from the

Re: [R] plot probability density function (pdf)

2011-10-12 Thread pigpigmeow
x11() what does it mean? if my data has missing value, can I plot the graph? -- View this message in context: http://r.789695.n4.nabble.com/plot-probability-density-function-pdf-tp3897055p3898362.html Sent from the R help mailing list archive at Nabble.com.

[R] exclude columns with at least three consecutive zeros

2011-10-12 Thread Samir Benzerfa
Hi everyone, I have a large data set with about 3'000 columns and I would like to exclude all columns which include three or more consecutive zeros (see below example). A further issue is that it should just jump NA values if any. How can I do this? In the below example R should exclude

[R] Read in cluster files (cdt, atr, gtr)

2011-10-12 Thread Mark Ebbert
Does anyone know if there is a method available to read in cluster files (cdt, atr and gtr)? I found one method (http://bioinformatics.holstegelab.nl/manuals/R/library/integromicsMethods/html/iMethods.read.tv.html), but it doesn't seem to create an object that can be used by heatmap. The

Re: [R] R and Forex

2011-10-12 Thread R. Michael Weylandt
This being what exactly? Traded in FX using R? Yes, its done everyday, even as I type Michael On Wed, Oct 12, 2011 at 8:10 AM, Yves S. Garret yoursurrogate...@gmail.com wrote: No, that's not what I meant.  I was curious if anyone has ever done this before and how well it worked.  Any tips

[R] ARMA and prediction

2011-10-12 Thread Ritz
Hello, I am running an ARMA model to run forecast for changes in SP 500 prices. My ARMA calculations look as follows armacal - arma( spdata, order = c(0,4), lag = list(ma = c(1,2,4)) ) Output: Call: arma(x = spdata, order = c(0, 4), lag = list(ma = c(1, 2, 4)) ) Coefficient(s): ma1

[R] treat NA's as zero when summed up with numbers

2011-10-12 Thread Samir Benzerfa
Hello everybody, is there any way to treat NA's as zero when they are summed up with numbers, but to treat them as NA's when summed up only with NA's. Specifically want that: 5+NA=5, but NA+NA=NA (and not zero). Any ideas? Best, S.B. [[alternative HTML version deleted]]

Re: [R] plot probability density function (pdf)

2011-10-12 Thread Bert Gunter
Please -- this is not your personal help advisor. Use R's Help facilities before posting. -- Bert On Wed, Oct 12, 2011 at 8:20 AM, pigpigmeow gloryk...@hotmail.com wrote: x11() what does it mean? ?x11 if my data has missing value, can I plot the graph? Try it and see. -- View this

[R] SARIMAX as extended AR(I)MAX?

2011-10-12 Thread testrider
I am currently trying to make a forecast based on past observations of the dependent variable AND external variables at the same time. I know that ARIMAX allows you to do this, however when I use this function it fits the model using the last k lags. What i actually want is to decide on the best

Re: [R] CVbinary - Help

2011-10-12 Thread Weidong Gu
Maybe it was caused by your modeling binary variable using lm rather than glm. Weidong Gu On Wed, Oct 12, 2011 at 9:59 AM, anamiguita anamigu...@gmail.com wrote: Hey, I need some help. I want to obtain a cross validation for a regression model (binary response) but I got an error with

Re: [R] Generelized Negative Binomial model in R

2011-10-12 Thread Steve Lianoglou
Hi, On Wed, Oct 12, 2011 at 11:23 AM, Akram Khaleghei Ghosheh balagh a.khaleg...@gmail.com wrote: Hello; Does anybody knows that R have a function for Generelized Negative Binomial model, something like gnbreg in STATA where dispersion parameter itself is a function of covaraites ? Take a

Re: [R] Nonlinear regression aborting due to error

2011-10-12 Thread Jean V Adams
Dennis Fisher fis...@plessthan.com wrote on 10/12/2011 08:06:12 AM: jean initial values: INITEMAX - -25 INITEFFECT - 25 INITC50 - 14 GAMMA - INITGAMMA - 500 see below for other issues. dennis Dennis Fisher MD P (The P Less

Re: [R] Minimum cutsets

2011-10-12 Thread Gábor Csárdi
Hi Mohammed, http://igraph.sourceforge.net/doc/R/graph.maxflow.html For directed graphs, and s-t cuts you need the development version, from igraph.sf.net. The source code is either here: http://cran.r-project.org/web/packages/igraph/index.html or here:

Re: [R] loop function within a loop

2011-10-12 Thread Weidong Gu
It's better to avoid loop in this situation. If you want to reorder subsets of the data based on event, the follow works df-read.table('clipboard',header=TRUE) sp.or-lapply(split(df,df$group),function(ldf) ldf[order(ldf$event),]) new.df-do.call('rbind',sp.or) Weidong Gu On Wed, Oct 12, 2011 at

Re: [R] exclude columns with at least three consecutive zeros

2011-10-12 Thread William Dunlap
First define a function that returns TRUE if a column should be dropped. E.g., has3Zeros.1 - function(x) { x - x[!is.na(x)] == 0 # drop NA's, convert 0's to TRUE, others to FALSE if (length(x) 3) { FALSE # you may want to further test short vectors } else {

Re: [R] treat NA's as zero when summed up with numbers

2011-10-12 Thread David Winsemius
On Oct 12, 2011, at 12:45 PM, Samir Benzerfa wrote: Hello everybody, is there any way to treat NA's as zero when they are summed up with numbers, but to treat them as NA's when summed up only with NA's. Specifically want that: 5+NA=5, but NA+NA=NA (and not zero). sum(x , na.rm=TRUE)

Re: [R] high and lowest with names

2011-10-12 Thread Ben qant
Hello, This is my solution. This is pretty fast (tested with a larger data set)! If you have a more elegant way to do it (of similar speed), please reply. Thanks for the help! ## get highest and lowest values and names of a matrix # create sample data x - swiss$Education[1:25]

Re: [R] loop function within a loop

2011-10-12 Thread David Winsemius
On Oct 12, 2011, at 10:55 AM, Sally Zhen wrote: Hi all, I'm working on a loop function for a large dataset which contains 1000 different groups. I would like to reconstruct the order of events within each group by using a loop function in R. Not generally a good first strategy in R.

Re: [R] treat NA's as zero when summed up with numbers

2011-10-12 Thread David Winsemius
On Oct 12, 2011, at 1:33 PM, David Winsemius wrote: On Oct 12, 2011, at 12:45 PM, Samir Benzerfa wrote: Hello everybody, is there any way to treat NA's as zero when they are summed up with numbers, but to treat them as NA's when summed up only with NA's. Specifically want that:

[R] Minimization/Optimization under functional constraints

2011-10-12 Thread forget_f1
Hi, I hope someone can help me with the following issue. I need find the minimum beta that satisfies the following: inf{beta0 | f(x+beta*f(x))*f(x)=0} where f() is a function and x is a sample statistic. Functions such as nlminb and constrOptim minimize a function and output the parameter

Re: [R] Minimum cutsets

2011-10-12 Thread malhomidi
Hi Gabor, Thanks for the reply. I'm actually working on directed graphs and using Windows. Please, send me the windows version of the source code. Regards, Mohammed -- View this message in context: http://r.789695.n4.nabble.com/Minimum-cutsets-tp885346p3899012.html Sent from the R

[R] Simple question regarding to specifying model

2011-10-12 Thread BvZ
I have hunted around but cannot find the command which allows me to specify parameters of a model. For example, model.m1 - nls(y ~ alpha * x1/(beta + x1), data = data, start = list(beta = 20, alpha = 120), trace = TRUE) This will estimate the parameters, which allows to investigate the

Re: [R] How to run Rcmdr with OS 10.4?

2011-10-12 Thread Robert Baer
As the posting guide tells you, the first thing to try is to install the latest version of R. -Original Message- From: roche...@free.fr Sent: Tuesday, October 11, 2011 5:44 AM To: r-help@r-project.org Subject: [R] How to run Rcmdr with OS 10.4? I've installed Rcmdr package and it

Re: [R] Simple question regarding to specifying model

2011-10-12 Thread David Winsemius
On Oct 12, 2011, at 2:25 PM, BvZ wrote: I have hunted around but cannot find the command which allows me to specify parameters of a model. For example, model.m1 - nls(y ~ alpha * x1/(beta + x1), data = data, start = list(beta = 20, alpha = 120), trace = TRUE) This will estimate the

[R] Tinn-R change editor background color

2011-10-12 Thread Ben qant
Hello, Does anyone know how to change the Tinn-R editor background color? White is rough on the eyes... Thanks, Ben [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE

[R] strange behavior with R.dll

2011-10-12 Thread Erin Hodgess
Dear R People: I am using R-2.13.2 on a Windows 7 machine. I compiled from source on 32 bit a couple of weeks ago and the R.dll got removed by my anti-virus software. Same thing on 64 bit today. Is anyone having this problem, please? I'm using Norton AV. Thanks, Erin -- Erin Hodgess

Re: [R] Tinn-R change editor background color

2011-10-12 Thread Ben qant
Never mind: option color preference Sorry...overlooked that 10 times I guess. regards On Wed, Oct 12, 2011 at 12:54 PM, Ben qant ccqu...@gmail.com wrote: Hello, Does anyone know how to change the Tinn-R editor background color? White is rough on the eyes... Thanks, Ben

Re: [R] R and Forex

2011-10-12 Thread R. Michael Weylandt
As was pointed out to you before, this is really more of an R-SIG-Finance question, but I wouldn't expect too much explanation there either, just people pointing you to the standard R finance tools (quantmod, zoo/xts, TTR, RBloomberg, and the Rmetrics suite; there's also some fantastic tools in

Re: [R] strange behavior with R.dll

2011-10-12 Thread Duncan Murdoch
On 12/10/2011 3:03 PM, Erin Hodgess wrote: Dear R People: I am using R-2.13.2 on a Windows 7 machine. I compiled from source on 32 bit a couple of weeks ago and the R.dll got removed by my anti-virus software. Same thing on 64 bit today. Is anyone having this problem, please? I'm using

Re: [R] Simple question regarding to specifying model

2011-10-12 Thread BvZ
Tried this, it does not seem to work. It is really simple what I am trying to do. I have a pre-specified best-fit line, and wish to run some diagnostic tests for goodness of fit. I will play around with the predict function, thanks a lot David! -- View this message in context:

[R] large numbers

2011-10-12 Thread Veerappa Chetty
Hi, When I import an excel CSV file, large numbers such as 43988014.3 is imported as 43988014, leaving out the decimal .4. How to import keeping the fraction? Thanks. Chetty -- Professor of Family Medicine Boston University Tel: 617-414-6221, Fax:617-414-3345 emails:

Re: [R] large numbers

2011-10-12 Thread R. Michael Weylandt
Are you sure its being imported into R without the decimal and that it's not just a print option? I can't off the cuff think of a reason why that would happen... Try print(valueFromImport - 43988014) and see what you get Michael On Wed, Oct 12, 2011 at 3:16 PM, Veerappa Chetty

[R] reading large numbers

2011-10-12 Thread Veerappa Chetty
Hi, This happens when I read in large numbers; as.numeric(4398801.3) [1] 4398801 as.numeric(439880.3) [1] 439880.3 Please help to read in numbers with more than 8 characters! Thanks. Chetty -- Professor of Family

  1   2   >