search for: tirthankar

Hi, I have a simple beginner's question on removing a list of objects. Say I have objects C243.Daily1, C243.Daily2...C243.Daily5 in my workspace. I'd like to remove these without using rm five times. So I write. > a <- list(paste("C243.Daily",sep="",1:5)) > rm(a) Obviously this wouldn't work, as it would only remove the object a. But is there any way

...lt;- as.integer(substr(yyyymm,5,6)) if (month==12){ date <- as.Date(ISOdate(year+1,1,1))-1 }else{ date <- as.Date(ISOdate(year,month+1,1))-1 } print(date) } x$Date <- as.vector(sapply(x$YearEnd,GetEOM)) str(x$Date) Is there a simpler way to do this please? TIA and best, -Tir Tirthankar Patnaik India Strategy Citigroup Investment Research +91-22-6631 9887

On 04/11/2017 10:20 PM, Daniel Nordlund wrote: > Tirthankar, > > "random number generators" do not produce random numbers. Any given > generator produces a fixed sequence of numbers that appear to meet > various tests of randomness. By picking a seed you enter that sequence > in a particular place and subsequent numbers in the se...

...V5 5 $`10` V1 V5 1 5 $`11` integer(0) $`12` V5 5 > But now I'm stuck. I don't how to use this list of indices at each row to exclude my columns. Any pointers please? Would such an exercise be easier if I use a time-series based object, like a zoo. TIA and best, -Tir Tirthankar Patnaik India Strategy Citigroup Investment Research +91-22-6631 9887

...uth-TAOCP-2002");runif(1,17,26)>25.99}, > 1:10000) > [1] 280 415 826 1372 2224 2544 3270 3594 3809 4116 4236 5018 5692 7043 > 7212 7364 7747 9256 9491 9568 9886 > > > > Bill Dunlap > TIBCO Software > wdunlap tibco.com > > On Fri, Nov 3, 2017 at 10:31 AM, Tirthankar Chakravarty < > tirthankar.lists at gmail.com> wrote: > >> >> Bill, >> >> I have clarified this on SO, and I will copy that clarification in here: >> >> "Sure, we tested them on other 8-digit numbers as well & we could not >> replicate....

Hi, I have some confusion in applying a function over a column. Here's my function. I just need to shift non-March month-ends to March month-ends. Initially I tried seq.dates, but one cannot give a negative increment (decrement) here. return(as.Date(seq.dates(format(xdate,"%m/%d/%Y"),by="months",len=4)[4]) ) Hence this simple function: > mydate <-

Hi, A quick beginner's question. I have two time series, A with daily data, and another B with data at varying frequencies, but mostly annual. Both the series are sorted ascending. I need to merge these two series together in the following way: For any entry of A, the lookup should match with B until we find an entry of B that's larger than A's. For all A[i], i = 1,...,N and B[j],

Tirthankar, "random number generators" do not produce random numbers. Any given generator produces a fixed sequence of numbers that appear to meet various tests of randomness. By picking a seed you enter that sequence in a particular place and subsequent numbers in the sequence appear to be u...

Hi, Quick question: Say I have a date variable in a data frame or matrix, and I'd like to preserve the date format when using write.table. However, when I export the data, I get the generic number underlying the date, not the date per se, and a number such as 11323, 11324, etc are not meaningful in Excel. Is there any way I can preserve the format of a date on writing into a text-file? TIA

...at give about the same first value: > > > Filter(function(s){set.seed(s);runif(1,17,26)>25.99}, 1:10000) > [1] 514 532 1951 2631 3974 4068 4229 6092 6432 7264 9090 > > > > Bill Dunlap > TIBCO Software > wdunlap tibco.com > > On Fri, Nov 3, 2017 at 12:49 AM, Tirthankar Chakravarty < > tirthankar.lists at gmail.com> wrote: > >> This is cross-posted from SO (https://stackoverflow.com/q/4 >> 7079702/1414455), >> but I now feel that this needs someone from R-Devel to help understand why >> this is happening. >> >> We ar...

...> Filter(function(s){set.seed(s, kind="Knuth-TAOCP-2002");runif(1,17,26)>25.99}, 1:10000) [1] 280 415 826 1372 2224 2544 3270 3594 3809 4116 4236 5018 5692 7043 7212 7364 7747 9256 9491 9568 9886 Bill Dunlap TIBCO Software wdunlap tibco.com On Fri, Nov 3, 2017 at 10:31 AM, Tirthankar Chakravarty < tirthankar.lists at gmail.com> wrote: > > Bill, > > I have clarified this on SO, and I will copy that clarification in here: > > "Sure, we tested them on other 8-digit numbers as well & we could not > replicate. However, these are honest-to-goodnes...

Hi, I have a history dataset, a matrix with about 1590 obs, and 242 cols, and I need to update this matrix with an 'update' matrix that has about 30 rows, and roughly similar number of columns as the history ds (but not necessarily equal). The update dataset is read from an Excel ODBC connection. When I try and merge these datasets, I get counter-intuitive results. library(RODBC) chn

Dear list members, Is there a package somewhere for jointly estimating two poisson processes? I think the closest I've come is using the "SUR" option in the Zelig package (see below), but when I try the "poisson" option instead of the "SUR" optioin I get an error (error given below, and indeed, reading the documentation of the Zelig package, I get the impression

> On 5 Nov 2017, at 15:17 , Duncan Murdoch <murdoch.duncan at gmail.com> wrote: > > On 04/11/2017 10:20 PM, Daniel Nordlund wrote: >> Tirthankar, >> "random number generators" do not produce random numbers. Any given >> generator produces a fixed sequence of numbers that appear to meet >> various tests of randomness. By picking a seed you enter that sequence >> in a particular place and subsequent numbers...

This is cross-posted from SO (https://stackoverflow.com/q/47079702/1414455), but I now feel that this needs someone from R-Devel to help understand why this is happening. We are facing a weird situation in our code when using R's [`runif`][1] and setting seed with `set.seed` with the `kind = NULL` option (which resolves, unless I am mistaken, to `kind = "default"`; the default being

>>>>> Tirthankar Chakravarty <tirthankar.lists at gmail.com> >>>>> on Fri, 3 Nov 2017 13:19:12 +0530 writes: > This is cross-posted from SO > (https://stackoverflow.com/q/47079702/1414455), but I now > feel that this needs someone from R-Devel to help > understa...

...' was, but it is easy to come up with seeds that give about the same first value: > Filter(function(s){set.seed(s);runif(1,17,26)>25.99}, 1:10000) [1] 514 532 1951 2631 3974 4068 4229 6092 6432 7264 9090 Bill Dunlap TIBCO Software wdunlap tibco.com On Fri, Nov 3, 2017 at 12:49 AM, Tirthankar Chakravarty < tirthankar.lists at gmail.com> wrote: > This is cross-posted from SO (https://stackoverflow.com/q/47079702/1414455 > ), > but I now feel that this needs someone from R-Devel to help understand why > this is happening. > > We are facing a weird situation in our...

Can R invert a 160000x160000 matrix with all positive numbers? Thanks a lot!

...atterplot of two groups side-by-side? (nonunah at yahoo.de) > 8. Re: Nomogram with stratified cph in Design package (reneepark) > 9. Stochastic Gradient Ascent for logistic regression (Tim LIU) > 10. Re: Question of "Quantile Regression for Longitudinal Data" > (Tirthankar Chakravarty) > 11. Re: THE EQUIVALENT OF SQL INNER TABLE JOIN IN R > (Gabor Grothendieck) > 12. Re: THE EQUIVALENT OF SQL INNER TABLE JOIN IN R (Peter Dalgaard) > 13. Re: 3 questions regarding matrix copy/shuffle/compares (Esmail) > 14. Re: Conditional plot labels (bapt...

Martin, Thanks for the helpful reply. Alas I had forgotten that (implied) unfavorable comparisons of *nix systems with Windows systems would likely draw irate (but always substantive) responses on the R-devel list -- poor phrasing on my part. :) Regardless, let me try to address some of the concerns related to the construction of the MRE itself and try to see if we can clean away the shrubbery