Re: [R] How to re-combine values based on an index?

2012-12-02 Thread signal
Thank you I believe that is just what I was looking for. Brian On Dec 2, 2012, at 12:37 AM, Berend Hasselman [via R] ml-node+s789695n4651599...@n4.nabble.com wrote: On 02-12-2012, at 06:06, Brian Feeny wrote: Thank you for your response, here is a better example of what I am trying

Re: [R] Line numbers with errors and warnings?

2012-12-02 Thread Worik R
What I mean is how do I get the R compilation or execution process to spit out a line number with errors and warnings? I should not have mentioned ESS, that is a distraction. option(error=browser) is a help. But it still does not say what piece of code caused the error. This is costing me a

[R] How can i change the number of iterations in lme function

2012-12-02 Thread Haris Rhrlp
lme2-lme(y~ age + time + timeDay,data= DF, random=~time |id/timeDay) Error in lme.formula(y ~ age + time + timeDay, data = DF, random = ~time |  :    nlminb problem, convergence error code = 1   message = iteration limit reached without convergence (10) [[alternative HTML version

Re: [R] Nightingale’s Rose chart-any suggestion?

2012-12-02 Thread Jim Lemon
On 11/30/2012 11:24 PM, Henry Smith wrote: Hello, Everyone. Does anyone know how to create a Nightingale’s Rose chart by using R? Hopefully, the graph could be displayed like this: http://mbostock.github.com/protovis/ex/crimea-rose.html Hi Henry, I cited Florence in a paper with Anupam Tyagi

Re: [R] x axis labels

2012-12-02 Thread Jim Lemon
On 12/02/2012 09:49 AM, apprentice wrote: Hi, I have a problem with plotting my data as a histogram using barplot2. I am plotting the x axis using axis(). I want to label only every fifth bar and not all of them as they get too busy but I don't know how to do it. I was trying to set the labels

Re: [R] instrumental variables regression using ivreg (AER) or tsls (sem)

2012-12-02 Thread Achim Zeileis
On Thu, 29 Nov 2012, Ranjan Maitra wrote: Dear friends, I am trying to understand and implement instrumental variables regression using R. I found a small (simple) example here which purportedly illustrates the mechanics (using 2-stage least-squares):

Re: [R] [mgcv][gam] Manually defining my own knots?

2012-12-02 Thread Simon Wood
Hi Andrew, mgcv matches the knots to the smooth arguments by name. If an element of 'knots' has no name it will be ignored. The following will do what you want... dumb.example = gam(y~s(x,k=3),knots=list(x=dumb.knots)) best, Simon On 29/11/12 23:44, Andrew Crane-Droesch wrote: Dear List,

Re: [R] How can i change the number of iterations in lme function

2012-12-02 Thread Berend Hasselman
On 02-12-2012, at 10:15, Haris Rhrlp wrote: lme2-lme(y~ age + time + timeDay,data= DF, random=~time |id/timeDay) Error in lme.formula(y ~ age + time + timeDay, data = DF, random = ~time | : nlminb problem, convergence error code = 1 message = iteration limit reached without convergence

Re: [R] cubic spline

2012-12-02 Thread Hans W Borchers
but definitely *no* need to use a function from an extra CRAN package .. as someone else ``erronously'' suggested. Except that Matlab's interp1() 'cubic' method does not use cubic spline interpolation, but Moler's 'pchip' approach, a piecewise cubic Hermite interpolation. Thus the results are

Re: [R] Line numbers with errors and warnings?

2012-12-02 Thread Steve Lianoglou
Hi, On Sun, Dec 2, 2012 at 12:31 AM, Worik R wor...@gmail.com wrote: What I mean is how do I get the R compilation or execution process to spit out a line number with errors and warnings? As Duncan mentioned already, you can't *always* get a line number. You can, however, usually get enough

Re: [R] Line numbers with errors and warnings?

2012-12-02 Thread Milan Bouchet-Valat
Le dimanche 02 décembre 2012 à 06:02 -0500, Steve Lianoglou a écrit : Hi, On Sun, Dec 2, 2012 at 12:31 AM, Worik R wor...@gmail.com wrote: What I mean is how do I get the R compilation or execution process to spit out a line number with errors and warnings? Indeed, I often suffer from the

Re: [R] Line numbers with errors and warnings?

2012-12-02 Thread Duncan Murdoch
On 12-12-02 12:31 AM, Worik R wrote: What I mean is how do I get the R compilation or execution process to spit out a line number with errors and warnings? I should not have mentioned ESS, that is a distraction. I think ESS is not a distraction here, it is the cause of the problem. If you

Re: [R] Line numbers with errors and warnings?

2012-12-02 Thread Duncan Murdoch
On 12-12-02 8:33 AM, Milan Bouchet-Valat wrote: Le dimanche 02 décembre 2012 à 06:02 -0500, Steve Lianoglou a écrit : Hi, On Sun, Dec 2, 2012 at 12:31 AM, Worik R wor...@gmail.com wrote: What I mean is how do I get the R compilation or execution process to spit out a line number with errors

Re: [R] Getting all possible contingency tables

2012-12-02 Thread Christofer Bogaso
Thanks Ted (and other) for your suggestion. Here I have implemented following: Tab - matrix(c(8, 10, 12, 6), nr = 2) Simu_Number - 5 Tab_Simulate - vector(list, length = Simu_Number) for (i in 1:Simu_Number) { Tab_Simulate[[i]] - matrix(rmultinom(1, sum(Tab), rep(0.25, 4)), nrow =

Re: [R] reading json tables

2012-12-02 Thread Michael Friendly
On 12/1/2012 4:08 PM, Duncan Temple Lang wrote: Hi Michael The problem is that the content of the .js file is not JSON, but actual JavaScript code. You could use something like the following tt = readLines(http://mbostock.github.com/protovis/ex/wheat.js;) txt = c([, gsub(;, ,, gsub(var

Re: [R] reading json tables

2012-12-02 Thread Duncan Temple Lang
Hi Michael The actual result I want is two data frames, wheat and monarch, whereas fromJSON returns a list of lists. I'll try to figure that part out. do.call(rbind, data[[1]]) will do the job, but there are elements in each of data[[1]] and data[[2]] that are incomplete and which need

Re: [R] Line numbers with errors and warnings?

2012-12-02 Thread Milan Bouchet-Valat
Le dimanche 02 décembre 2012 à 09:02 -0500, Duncan Murdoch a écrit : On 12-12-02 8:33 AM, Milan Bouchet-Valat wrote: Le dimanche 02 décembre 2012 à 06:02 -0500, Steve Lianoglou a écrit : Hi, On Sun, Dec 2, 2012 at 12:31 AM, Worik R wor...@gmail.com wrote: What I mean is how do I get the

Re: [R] Getting all possible contingency tables

2012-12-02 Thread Ted Harding
On 02-Dec-2012 14:17:20 Christofer Bogaso wrote: Thanks Ted (and other) for your suggestion. Here I have implemented following: Tab - matrix(c(8, 10, 12, 6), nr = 2) Simu_Number - 5 Tab_Simulate - vector(list, length = Simu_Number) for (i in 1:Simu_Number) {

Re: [R] How to use doSMP(revoIPC) with R 2.15.x version

2012-12-02 Thread Uwe Ligges
On 29.11.2012 17:56, 박상규 wrote: Hello, I'd like to use package 'doSMP'. But I can only found source codes in the CRAN. I tried to build source code using 'r CMD build doSMP' in the source directory, the following error is shown: ERROR: dependencies 'foreach', 'iterators', 'revoIPC' are

[R] overlapping graphs in logarithmic y-axis

2012-12-02 Thread eliza botto
dear useRs, i want to overlap graphs of two matrices in such a way that the y-axis of graph should be logarithmic against normal x-axis. i am, unsuccessfully, trying the followings matplot(mata, log=mata,type = l, col=red)lines(mata, log=matb,type = l, col=yellow) could you please help me out

Re: [R] Getting all possible contingency tables

2012-12-02 Thread Christofer Bogaso
Thanks Ted for your correction. I was depressed thinking that I did not understand the theory. However now it comes as just a plotting mistake! Thanks, On 02 December 2012 20:47:48, (Ted Harding) wrote: On 02-Dec-2012 14:17:20 Christofer Bogaso wrote: Thanks Ted (and other) for your

Re: [R] overlapping graphs in logarithmic y-axis

2012-12-02 Thread Uwe Ligges
On 02.12.2012 16:19, eliza botto wrote: dear useRs, i want to overlap graphs of two matrices in such a way that the y-axis of graph should be logarithmic against normal x-axis. i am, unsuccessfully, trying the followings matplot(mata, log=mata,type = l, col=red)lines(mata, log=matb,type =

Re: [R] Stripchart colors don't vary after I sort a data frame

2012-12-02 Thread Uwe Ligges
On 28.11.2012 18:27, Thomas Levine wrote: # Hi, # This plot has two colors. overflow - read.csv('http://chainsaw.thomaslevine.com/overflow.csv', stringsAsFactors = F) png('original.png') stripchart(overflow$precipi ~ overflow$after.9.am, method='stack', pch = 22, bg = overflow$overflow + 1,

Re: [R] overlapping graphs in logarithmic y-axis

2012-12-02 Thread eliza botto
Dear ligges,thankyou very much. Actually majority of the questions i asked on this forum were without sense, as quoted by majority of experts. i actually wanted to over lap mata and matb. by replacing your code with matb, it worked!!! thankseliza Date: Sun, 2 Dec 2012 16:53:24 +0100 From:

Re: [R] overlapping graphs in logarithmic y-axis

2012-12-02 Thread Uwe Ligges
On 02.12.2012 17:03, eliza botto wrote: Dear ligges,thankyou very much. Actually majority of the questions i asked on this forum were without sense IOnm this case the reason is that you do not ask the question precisely. Your exampole had matplot(mata, .) lines(mata, ) Please

Re: [R] Plot from a jpeg

2012-12-02 Thread Uwe Ligges
On 27.11.2012 22:38, bgnumis wrote: Hi all, I know that I can for instance draw to plots in one using nf-layout(matrix(c(1,2),1,2,byrow=FALSE)) Imagine I have 3 files: plot1.jpeg plot2.jpeg plot3.jpeg Anyone knows if I can read them and put on one colum and three rows reading directly

[R] simple subset question

2012-12-02 Thread Felipe Carrillo
 Hi, Consider the small dataset below, I want to subset by two variables in one line but it wont work...it works though if I subset separately. I have to be missing something obvious that I did not realize before while using subset.. fish - structure(list(IDWeek = c(27L, 28L, 29L, 30L, 31L, 32L,

Re: [R] simple subset question

2012-12-02 Thread Anthony Damico
shouldn't you just change b to x and winter to fish? :) On Sun, Dec 2, 2012 at 12:21 PM, Felipe Carrillo mazatlanmex...@yahoo.comwrote: Hi, Consider the small dataset below, I want to subset by two variables in one line but it wont work...it works though if I subset separately. I have to

Re: [R] simple subset question

2012-12-02 Thread David Winsemius
On Dec 2, 2012, at 9:21 AM, Felipe Carrillo wrote: Hi, Consider the small dataset below, I want to subset by two variables in one line but it wont work...it works though if I subset separately. I have to be missing something obvious that I did not realize before while using subset..

Re: [R] simple subset question

2012-12-02 Thread R. Michael Weylandt
On Sun, Dec 2, 2012 at 5:21 PM, Felipe Carrillo mazatlanmex...@yahoo.com wrote: Hi, Consider the small dataset below, I want to subset by two variables in one line but it wont work...it works though if I subset separately. I have to be missing something obvious that I did not realize before

Re: [R] simple subset question

2012-12-02 Thread Gerrit Eichner
Hi, Felipe, two typos? See below! On Sun, 2 Dec 2012, Felipe Carrillo wrote:  Hi, Consider the small dataset below, I want to subset by two variables in one line but it wont work...it works though if I subset separately. I have to be missing something obvious that I did not realize before

[R] How to calculate mean of every nth time series data with zoo or xts ?

2012-12-02 Thread 박상규
Hello, I have 1-minute time series stock data and I'd like to calculate mean of every n-th candle data of m-days. result = c(mean of 1th data, mean of 2nd data, ...) mean of 1th data = (1th data of 2012-1-1 + 1th data of 2012-1-2 + 1th data of 2012-1-3) / 3 mean of 2nd data = (2nd data of

Re: [R] missed values

2012-12-02 Thread Joshua Ulrich
It's simpler, easier, and more correct to use zoo::na.approx. R library(zoo) R dat1-read.table(text= +date value + 80 2006-10 0.1577647 + 81 2006-11 NaN + 82 2006-12 NaN + 83 2007-01 NaN + 84 2007-02 NaN + 85 2007-03 0.2956429, +

Re: [R] simple subset question

2012-12-02 Thread Felipe Carrillo
Sorry, I was trying it to subset from a bigger dataset called 'winter' and forgot to change the variable names when I asked the question. David W suggestion works but the strange part is that I am still getting an error message with :   x - subset(fish,Year==2012 Total==max(Total));x I get: [1]

Re: [R] How to calculate mean of every nth time series data with zoo or xts ?

2012-12-02 Thread Rui Barradas
Hello, You should include a data example, using ?dput. Anyway, making up some data, t1 - as.POSIXct(2012-1-1 0:0:0) t2 - as.POSIXct(2012-1-4 0:0:0) d - seq(t1, t2, by = 1 min) x - rnorm(length(d)) z - zoo(x, d) aggregate(z, list(format(index(z), %H:%M)), FUN = mean) Hope this helps, Rui

Re: [R] simple subset question

2012-12-02 Thread Felipe Carrillo
Works with the small dataset (2 years) but I get the error message with the whole dataset (12 years of data). I am going to have to check what's wrong with it...Thanks Felipe D. Carrillo Supervisory Fishery Biologist Department of the Interior US Fish Wildlife Service California, USA

Re: [R] simple subset question

2012-12-02 Thread R. Michael Weylandt
On Sun, Dec 2, 2012 at 6:46 PM, Felipe Carrillo mazatlanmex...@yahoo.com wrote: Works with the small dataset (2 years) but I get the error message with the whole dataset (12 years of data). I am going to have to check what's wrong with it...Thanks Off the cuff guess: there's a NA in Total so

Re: [R] simple subset question

2012-12-02 Thread William Dunlap
I am still getting an error message with : x - subset(fish,Year==2012 Total==max(Total));x I get: [1] IDWeek Total FrySmolt FryEq Year 0 rows (or 0-length row.names) The above is not an error message. It says that there are no rows satisfying your criteria. Note that

Re: [R] Line numbers with errors and warnings?

2012-12-02 Thread Duncan Murdoch
On 12-12-02 9:52 AM, Milan Bouchet-Valat wrote: Le dimanche 02 décembre 2012 à 09:02 -0500, Duncan Murdoch a écrit : On 12-12-02 8:33 AM, Milan Bouchet-Valat wrote: Le dimanche 02 décembre 2012 à 06:02 -0500, Steve Lianoglou a écrit : Hi, On Sun, Dec 2, 2012 at 12:31 AM, Worik R

Re: [R] simple subset question

2012-12-02 Thread Felipe Carrillo
  Using my whole dataset I get: library(plyr) ddply(winter,Year,summarise,maxTotal=max(Total))    fish - structure(list(Year = 2002:2012, maxTotal = c(1464311L, 1071051L, 714837L, 2115018L, 850491L, 207537L, 321195L, 935599L, 194429L, 157260L, 303259L)), .Names = c(Year, maxTotal), row.names =

Re: [R] e1071 SVM: Cross-validation error confusion matrix

2012-12-02 Thread Frank Harrell
What do you mean by accuracy? Proportion classified correctly is not a good index of accuracy if that's the problem. Frank rahul143 wrote Hi, I ran two svm models in R e1071 package: the first without cross-validation and the second with 10-fold cross-validation. I used the following

Re: [R] simple subset question

2012-12-02 Thread David Winsemius
The reason I suggested the alternative that I did was because your code could fail when the max(Total) value was not in the subset where Year==2012. -- David On Dec 2, 2012, at 11:34 AM, Felipe Carrillo wrote: Using my whole dataset I get: library(plyr)

Re: [R] simple subset question

2012-12-02 Thread Felipe Carrillo
I also thought that I could use it but ggplot didn't like it by itself,,,it wants to whole rowI guess it is hard to understand because I didn't post all the ggplot code of what I am trying to accomplish. The strangest thing is that ddply can extract all my Totals by year but when I try to

Re: [R] Line numbers with errors and warnings?

2012-12-02 Thread Milan Bouchet-Valat
Le dimanche 02 décembre 2012 à 14:21 -0500, Duncan Murdoch a écrit : On 12-12-02 9:52 AM, Milan Bouchet-Valat wrote: Le dimanche 02 décembre 2012 à 09:02 -0500, Duncan Murdoch a écrit : On 12-12-02 8:33 AM, Milan Bouchet-Valat wrote: Le dimanche 02 décembre 2012 à 06:02 -0500, Steve

Re: [R] cubic spline

2012-12-02 Thread Steve Stephenson
Dear all, many thanks for your answer and support!! Steve -- View this message in context: http://r.789695.n4.nabble.com/cubic-spline-tp4651537p4651698.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org mailing

[R] What is ts command for?

2012-12-02 Thread Hard Core
Hello, i'm here to have an explenation about time series I have to analize a time series of a stock prices and i don't understand the meaning of ts .. i load, for instance mps - get.hist.quote(instrument = bmps.mi , quote=Close, start=2006-01-01, end=2012-09-30, compression=d) and i get a

[R] function stays in loop

2012-12-02 Thread Mkkl
function - function(p){ i - 2 r - 0 while(i p) { q - (p/i) while(q=1) { q=(q-1) } if (q==0) { r=1 } else i=i+1 } return(r) } This function stays in some kind of loop

[R] Fwd: How to calculate different groups of varialbes importance level?

2012-12-02 Thread Solmaz Filiz KARABAĞ
Dear R user! I have a small question! I have calculated the relative importance of the variables. However I would like to compare the relative importance of two different groups of variables (i.e Strategy and industry) For example let me say that strategy has 2 sub varialbes and industry has

Re: [R] Daily Time Series, patterns.

2012-12-02 Thread arun
HI, I guess you can check library(forecast). You could also use ?stl(), though not sure whether this will help you in this case. For e.g. set.seed(5) quantity-sample(c(120:220,NA),699,replace=TRUE) Date=seq(as.Date(2011-01-01),len=699,by=1 day) dat3-data.frame(Date=Date,quantity=quantity)

[R] help setting up crossed data

2012-12-02 Thread derrick
Hello, and thanks for your time reading this. I'm trying to test interactions of my dataset, in which the all of the factors are within the same column. TypeVol 1 CMass -4.598 2 BBack -4.605 3 BMass -4.602 4 CMass -4.601 5 CBack -4.605 6

Re: [R] CreateThread failure since R 2.15.2 (32-bit)

2012-12-02 Thread Martin Studer
I tried with Java 1.6.0_20 and 1.7.0_09. It actually doesn't matter what command I'm running with system. Any command will fail with the same error message. But for the example I posted, yes, I was using ls from Rtools. I found that I can only reproduce this with R 2.15.2+ 32-bit on Windows

Re: [R] simple subset question

2012-12-02 Thread arun
Hi, I am getting this: x-subset(fish,Year==2012 Total==max(Total))  x #   IDWeek  Total   Fry  Smolt  FryEq Year #21 47 303259 34008 269248 491733 2012 A.K. - Original Message - From: Felipe Carrillo mazatlanmex...@yahoo.com To: R. Michael Weylandt michael.weyla...@gmail.com Cc:

Re: [R] How to calculate mean of every nth time series data with zoo or xts ?

2012-12-02 Thread arun
Hi, You could use: set.seed(5) dat1-data.frame(Date=seq(from=as.POSIXct(2012-11-01 00:00:00),to=as.POSIXct(2012-11-03 23:59:00),  by=1 min),col2=rnorm(4320,0,1)) library(xts) dat2-xts(dat1[,-1],order.by=dat1[,1]) res-ave(dat2[,1],format(index(dat2),%H:%M),FUN=mean) #or

[R] Warning message: In scan(file, what, nmax...)

2012-12-02 Thread F86
Dear R-users, When i try - Data1-read.table(/Users/kama/Analysis/GDP10.csv,header=TRUE,sep=;) i am getting this error: Warning message: In scan(file, what, nmax, sep, dec, quote, skip, nlines, na.strings, : number of items read is not a multiple of the number of columns I wonder what Iam

[R] Repeated-measures anova with a within-subject covariate (or varying slopes random-effects?)

2012-12-02 Thread Giuseppe Pagnoni
Dear all, I am having quite a hard time in trying to figure out how to correctly spell out a model in R (a repeated-measures anova with a within-subject covariate, I guess). Even though I have read in the posting guide that statistical advice may or may not get an answer on this list, I decided

[R] Change case of factor in data frame

2012-12-02 Thread Audrey
I am trying to write a function to change the case of all of the text in a data frame to lower case. I do not have foreknowledge of the data frame names or the data types of each column. It seems that if one references the data frame by index, then it returns class data.frame but if it is

[R] marginal structural modeling

2012-12-02 Thread John Sorkin
Does anyone know an R package that implements marginal structural modeling as described by Jamie Robins? A similar question was asked of the list approximately two-years ago and it does not appear that anyone responded. Thanks, John Sorkin John David Sorkin M.D., Ph.D. Chief, Biostatistics

Re: [R] simple subset question

2012-12-02 Thread arun
Hi, From the ddply() output, you could get the whole row by:  fish1 - structure(list(Year = 2002:2012, maxTotal = c(1464311L, 1071051L, 714837L, 2115018L, 850491L, 207537L, 321195L, 935599L, 194429L, 157260L, 303259L)), .Names = c(Year, maxTotal), row.names = c(NA, -11L), class = data.frame)

Re: [R] Line numbers with errors and warnings?

2012-12-02 Thread John Sorkin
Gentleman, This thread has been of great interest. Perhaps I missed part of it, but do far I have not seen an example of code that has line numbers that demonstrates how one can (in some instances) recover the line number of an error. Can I impose upon the people who contributed to this thread to

Re: [R] Line numbers with errors and warnings?

2012-12-02 Thread Duncan Murdoch
On 12-12-02 5:02 PM, John Sorkin wrote: Gentleman, This thread has been of great interest. Perhaps I missed part of it, but do far I have not seen an example of code that has line numbers that demonstrates how one can (in some instances) recover the line number of an error. Can I impose upon the

Re: [R] simple subset question

2012-12-02 Thread Felipe Carrillo
Actually, I dput(winter) and is not that big.so here is the whole thing I am trying to accomplish.  winter - structure(list(IDWeek = c(27L, 28L, 29L, 30L, 31L, 32L, 33L, 34L, 35L, 36L, 37L, 38L, 39L, 40L, 41L, 42L, 43L, 44L, 45L, 46L, 47L, 48L, 49L, 50L, 51L, 52L, 27L, 28L, 29L, 30L, 31L,

[R] Problem with figures

2012-12-02 Thread Shige Song
I am having problem making ggplot2, tikzDevice, and knitr working together. I used a very simple example: ---example.Rnw- \documentclass[preview]{standalone} \begin{document} \begin{figure} fig1,eval=TRUE,echo=FALSE,dev='tikz'= library(ggplot2)

Re: [R] Line numbers with errors and warnings?

2012-12-02 Thread Steve Lianoglou
Similar to Duncan's example, if you have a script test.R which looks like so: start script = a1 - 1:10 a2 - 101:122 plot(a1, a1) plot(a1, a2) end script == You can source it one way: R source('test.R', keep.source=TRUE) Error in xy.coords(x, y, xlabel, ylabel, log) : 'x' and

Re: [R] simple subset question

2012-12-02 Thread David L Carlson
As David W. guessed. The maximum is in year 2005 not 2012 so no row from 2012 matches the maximum. subset(winter,Year==2012 Total==max(Total)) [1] IDWeek Total FrySmolt FryEq Year 0 rows (or 0-length row.names) winter[which(winter$Total==max(winter$Total)),] IDWeek Total Fry

[R] Object Browser

2012-12-02 Thread rahul143
Dear all, I have tried all the popular R IDE or editors like Eclipse, Komodo, JGR, Revolution... They all have fancy fucntions like auto completion, syntax highlight BUT, I JUST WANT A OBJECT BROWSER! The easiest way to view objects in R console is fix(), but you have no global view of

[R] Can you turn a string into a (working) symbol?

2012-12-02 Thread rahul143
Dear folks-- Suppose I have an expression that evaluates to a string, and that that string, were it not a character vector, would be a symbol. I would like a function, call it doppel(), that will take that expression as an argument and produce something that functions exactly like the symbol

[R] Conjoint Analysis in R??

2012-12-02 Thread rahul143
Pls advise how I can use R in conjoint analysis?? regds Faisal Afzal Siddiqui Karachi, Pakistan Looking for last minute shopping deals? - TO GET MORE DETAILS CLICK HERE -- View this

[R] postForm() in RCurl and library RHTMLForms

2012-12-02 Thread rahul143
Hi RUsers, Suppose I want to see the data on the website url - http://www.nseindia.com/content/indices/ind_histvalues.htm; for the index SP CNX NIFTY for dates FromDate=01-11-2010,ToDate=02-11-2010 then read the html table from the page using readHTMLtable() I am using this code webpage

[R] Stepwise analysis with fixed variables

2012-12-02 Thread rahul143
Hello, How can I run a backward stepwise regression with part of the variables fixed, while the others participate in the backward stepwise analysis? Thank you - TO GET MORE DETAILS CLICK HERE -- View this message in context:

[R] Reading PDF files

2012-12-02 Thread rahul143
I need to do text mining on PDF files. I understand there is a readPDF command in tm that can be used. Have read the 2008 posts on converting PDF files to text by Tony Breyal and others. Wondering if the procedure has been standardized in any tutorial or otherwise? Being new to R, I was

[R] Problems with nls

2012-12-02 Thread rahul143
I'm trying to fit the Bass Diffusion Model using the nls function in R but I'm running into a strange problem. The model has either two or three parameters, depending on how it's parameterized, p (coefficient of innovation), q (coefficient of immitation), and sometimes m (maximum market

[R] postForm() in RCurl and library RHTMLForms

2012-12-02 Thread rahul143
Hi RUsers, Suppose I want to see the data on the website url - http://www.nseindia.com/content/indices/ind_histvalues.htm; for the index SP CNX NIFTY for dates FromDate=01-11-2010,ToDate=02-11-2010 then read the html table from the page using readHTMLtable() I am using this code webpage

[R] Sullivan, Timmerman and White 1999: TA rules, and R

2012-12-02 Thread rahul143
Friends I am trying to save myself some tedious work. I am processing a paper from The Journal Of Finance * Vol. LIV, No. 5 October 1999 by Sullivan, Timmerman and White. Data-Snooping, Technical Trading Rule Performance, and the Bootstrap I am aiming to reproduce their results using

[R] e1071 SVM: Cross-validation error confusion matrix

2012-12-02 Thread rahul143
Hi, I ran two svm models in R e1071 package: the first without cross-validation and the second with 10-fold cross-validation. I used the following syntax: #Model 1: Without cross-validation: svm.model - svm(Response ~ ., data=data.df, type=C-classification, kernel=linear, cost=1)

[R] Problems of metafile plots when converting word to pdf file

2012-12-02 Thread rahul143
Hi all, I copy metafile boxplot from R to Word. Then save as Pdf file. But I found there are some unexpected black lines in some plots within this PDF file. Please give me some advice how to get rid of these black lines in PDF files. Thanks a lot, Jim - TO GET MORE DETAILS CLICK

[R] RMySQL install on windows

2012-12-02 Thread rahul143
I have been trying to install RMySQL on Windows 7 following the procedure at: http://biostat.mc.vanderbilt.edu/wiki/Main/RMySQL I think I have properly installed RTools and created a proper Renviron.site file saying: MYSQL_HOME=C:/Program Files/MySQL/MySQL Server 5.5 When I try to install

[R] Object Browser

2012-12-02 Thread rahul143
What's the best object browser? Dear all, I have tried all the popular R IDE or editors like Eclipse, Komodo, JGR, Revolution... They all have fancy fucntions like auto completion, syntax highlight BUT, I JUST WANT A OBJECT BROWSER! The easiest way to view objects in R console is

[R] Linking C/C++ GUI to R.dll

2012-12-02 Thread rahul143
Hi, I am in the process of linking a C/C++ application to the R.dll directly. I have obtained the R source code and compiled it successfully. I have also successfully linked the R.dll directly into our application and have made calls successfully into the R.dll that are included in the

[R] Sullivan, Timmerman and White 1999: TA rules, and R

2012-12-02 Thread rahul143
Friends I am trying to save myself some tedious work. I am processing a paper from The Journal Of Finance * Vol. LIV, No. 5 October 1999 by Sullivan, Timmerman and White. Data-Snooping, Technical Trading Rule Performance, and the Bootstrap I am aiming to reproduce their results using

[R] residual and null deviance of an lme object with correlation structure

2012-12-02 Thread rahul143
Hello, I am attempting to calculate the residual and null deviance of an lme object that includes a corAR1 correlation structure. I tried deviance(lme.object) and it only returned NULL. Can anyone help? Thank you. - TO GET MORE DETAILS CLICK HERE -- View this message in

[R] finding index of maximum value in vector

2012-12-02 Thread rahul143
I found: max.col(matrix(c(1,3,2),nrow=1)) Is there a more concise/elegant way? Thanks, - TO GET MORE DETAILS CLICK HERE -- View this message in context: http://r.789695.n4.nabble.com/finding-index-of-maximum-value-in-vector-tp4651663.html Sent from the R help mailing list

[R] Fitting binomial lmer-model, high deviance and low logLik

2012-12-02 Thread rahul143
Hello I have a problem when fitting a mixed generalised linear model with the lmer-function in the Matrix package, version 0.98-7. I have a respons variable (sfox) that is 1 or 0, whether a roe deer fawn is killed or not by red fox. This is expected to be related to e.g. the density of red

[R] Hmisc latex cell background color

2012-12-02 Thread rahul143
Dear latex/R-Sweavers, Using the codel below, I can color text in individual cells for latex output. Is there a similar way to get a background shading? My attempts failed because I did not get the closing brace at the right place with Hmisc/latex. library(Hmisc) x -

[R] Kalman Filter Forecast using 'SSPIR'

2012-12-02 Thread rahul143
Dear R Users, I am new to state-space modeling. I am using SSPIR package for Kalman Filter. I have a data set containing one dependent variable and 7 independent variables with 250 data points. I want to use Kalman Filter for forecast the future values of the

[R] Fitting binomial lmer-model, high deviance and low logLik

2012-12-02 Thread rahul143
Hello I have a problem when fitting a mixed generalised linear model with the lmer-function in the Matrix package, version 0.98-7. I have a respons variable (sfox) that is 1 or 0, whether a roe deer fawn is killed or not by red fox. This is expected to be related to e.g. the density of red

[R] How to simulate correlated data

2012-12-02 Thread rahul143
Hello there, I would like to simulate X --Normal (20, 5) Y-- Normal (40, 10) and the correlation between X and Y is 0.6. How do I do it in R? Thank you very much - TO GET MORE DETAILS CLICK HERE -- View this message in context:

[R] concatenating expressions and standard text

2012-12-02 Thread rahul143
Hi all, is it possible to concatenate expressions and basic text when for instance labeling axis of a plot? I would like to see something like the concatenation of expression(C[0]) and for case 1 on my x axis. Obviously a plot(x, y, xlab=paste(expression(C[0]), in case1)) will not work.

Re: [R] superimpose density line over hist

2012-12-02 Thread rahul143
Romain == Romain Francois [hidden email] on Tue, 13 Dec 2005 15:40:59 +0100 writes: Romain A few comments : Romain - your code should be reproductible, otherwise it is useless. (that Romain recommandation is on the posting guide) Romain - that

Re: [R] Labeling a range of bars in barplot?

2012-12-02 Thread rahul143
Marc Schwartz (via MN) wrote: On Tue, 2005-12-13 at 10:53 +, Dan Bolser wrote: Hi, I am plotting a distribution of (ordered) values as a barplot. I would like to label groups of bars together to highlight aspects of the distribution. The label for the group should be the range of values

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

2012-12-02 Thread Adam Zeilinger
Dear Ravi, Thank you so much for the help. I switched to using the optimx function but I continue to use the spg method (for the most part) because I found that only spg consistently converges give different datasets. I also decided to use AIC rather that a likelihood ratio test. I have a

[R] R beginner

2012-12-02 Thread Akrem Zoghlami
Dear R-help group, I'm a R beginner and I find difficulty in manipulating R. Could you send to me a support that helps me to be familiar with R. Thank you in advance Sincerely, [[alternative HTML version deleted]] __ R-help@r-project.org

Re: [R] problem with convergence in mle2/optim function

2012-12-02 Thread Adam Zeilinger
Dear Berend, Thank you so much for your help! I was able to write the gradient function for my NLL function. For you're and others' possible interest, here is my final gradient function: Following from my description of the problem below, gr.p1, gr.p2, gr.mu1, and gr.mu2 are the (very

Re: [R] Warning message: In scan(file, what, nmax...)

2012-12-02 Thread F86
In addition, I tried Data1-read.table(/Users/kama/Analysis/GDP10.csv, sep=\t, fill=TRUE) the result: *** caught segfault *** address 0xd020, cause 'memory not mapped' Possible actions: 1: abort (with core dump, if enabled) 2: normal R exit 3: exit R without saving workspace 4: exit R

Re: [R] Warning message: In scan(file, what, nmax...)

2012-12-02 Thread F86
Also tries this: x-count.fields(/Users/kama/Analysis/GDP10.csv, sep=\t) x [1] 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 [38] 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 [75] 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1

Re: [R] simple subset question

2012-12-02 Thread Felipe Carrillo
I finally see what you and David W mean...I was under the assumption that by subsetting Year==2012 it would also extract the max for that year.I guess I have to review 'subset' again :) Thank you much Felipe D. Carrillo Supervisory Fishery Biologist Department of the Interior US Fish

Re: [R] Problem with figures

2012-12-02 Thread Duncan Murdoch
On 12-12-02 5:42 PM, Shige Song wrote: I am having problem making ggplot2, tikzDevice, and knitr working together. I used a very simple example: I don't use knitr so I can't really help, but you didn't tell us how you passed this file to knitr, so maybe nobody can. However, if you were

Re: [R] Can you turn a string into a (working) symbol?

2012-12-02 Thread R. Michael Weylandt
I believe I answered this a few weeks ago at this link: https://stat.ethz.ch/pipermail/r-help/2012-November/328053.html and following. Michael On Sun, Dec 2, 2012 at 2:43 PM, rahul143 rk204...@gmail.com wrote: Dear folks-- Suppose I have an expression that evaluates to a string, and that that

Re: [R] RMySQL install on windows

2012-12-02 Thread Gabor Grothendieck
On Sun, Dec 2, 2012 at 11:35 AM, rahul143 rk204...@gmail.com wrote: I have been trying to install RMySQL on Windows 7 following the procedure at: http://biostat.mc.vanderbilt.edu/wiki/Main/RMySQL I think I have properly installed RTools and created a proper Renviron.site file saying:

Re: [R] Problem with figures

2012-12-02 Thread Shige Song
Easiest way: copy and paste the code into Rstudio and hit compile pdf. From the command line, I believe you can do knit2pdf example.Rnw. Shige On Sun, Dec 2, 2012 at 6:12 PM, Duncan Murdoch murdoch.dun...@gmail.comwrote: On 12-12-02 5:42 PM, Shige Song wrote: I am having problem making

Re: [R] finding index of maximum value in vector

2012-12-02 Thread Jorge I Velez
?which.max On Mon, Dec 3, 2012 at 3:44 AM, rahul143 rk204...@gmail.com wrote: I found: max.col(matrix(c(1,3,2),nrow=1)) Is there a more concise/elegant way? Thanks, - TO GET MORE DETAILS CLICK HERE -- View this message in context:

  1   2   >