similar to: sensitivity tests fo causal inference

Hi all, does anyone know if there exists some library that implements the sensitivity tests for matched pairs and unmatched groups proposed in Rosenbaum's Observational Studies (2002: ch.4)? Cheers, Holger

SELECTION FOR CONSULTANTS BY THE WORLD BANK REQUEST FOR EXPRESSIONS OF INTEREST All Submissions must be made via the World Bank Group's "eConsult2" system. Please go to: https://wbgeconsult2.worldbank.org/wbgec/index.html Click on "Business Opportunities", and find the following assignment: Assignment Title: 1067014 - Statistical Analyst with Extensive Experience in

You might look at "eval" and "parse". Reid Huntsinger -----Original Message----- From: Spencer Graves [mailto:spencer.graves at pdf.com] Sent: Tuesday, March 11, 2003 4:03 PM To: janet rosenbaum Cc: r-help Subject: Re: [R] different file types The command "'source' causes R to accept its input from the named file". I don't know what you mean by

I have a question for our European readers: how common is it to use commas as decimal points in spread-sheet and statistics applications? Is it an inconvenience to require that all data use a period as decimal point? (i.e., 3.14159 rather than 3,14159). We're trying to make our program as foolproof as possible, and would prefer not to give users a chance to have commas as both decimal

I want to eliminate certain observations in a large dataframe (21000x100). I have written code which does this using a binary vector (0=delete obs, 1=keep), but it uses for loops, and so it's slow and in the extreme it causes R to hang for indefinite time periods. I'm looking for one of two things: 1. A document which discusses how to avoid for loops and situations in which it's

Dear all, I am using the Instrumental Variable approach to estimate the causal effects of TWO endogenous variables in a Mendelian Randomization study. As long as point estimation is concerned, I have no problem: both "ivreg" in library "AER" and "tsls" in library "sem" do the job perfectly. The problems begin when I try to obtain confidence intervals for

I have a semantics question, I am reading Bowerman and O'Connell and they state that forecasting models fall into two categories, univariate and causal. My question is whether a time series regression model that relates the time series of interest to functions of time such as the day of the week or month of the year could be considered a causal model? thanks, spencer [[alternative HTML

Windows (both XP and 98) seems to have a problem with our GUI, written using R's tcl/tk package. When we open another window, and then return to the GUI, the GUI has old pieces of text left from the previous windows and/or places which are cut off where a window had been. Also, when it computes something, the window will often decide to move itself over for no apparent reason, so it moves

Hi netters. Assume Y = {Y1,.....Yn} and X = {X1......Xm}, where Yi and Xi are random variables that can take on discrete values from V={0,1,2}. Each Yi in Y has some (0-k) parent variables in X, which means given the values of the parent variables (Xi0....Xik) the values of Yi is set. Yi = F(Xi0,...Xik), where F is a mapping function from parents to sons. Considering there are some noise in

Hi R-list, I am new to optimization in R and would appreciate help on the following question. I would like to minimize the following function using two constraints: ###### fn <- function(par,H,F){ fval <- 0.5 * t(par) %*% H %*% par + F%*% par fval } # matrix H is (n by k) # matrix F is (n by 1) # par is a (n by 1) set of weights # I need two constraints: # 1.

Hi. I want to use R with very large files, a couple hundred megabytes, but R crashes every time that I try. Reading help files seems to indicate that R ought to manage its memory itself. I know I have enough memory since stata handles these files perfectly well. I have a mac running OS 10.3 and am running RAqua 1.8.1. Is there anything I can do to make it deal with these files successfully?

Hi, I am looking for an elegant way to transform a vector into percentages of values that meet certain criteria. store<-c(1,1.4,3,1.1,0.3,0.6,4,5) # now I want to get the precentages of values # that fall into the categories <=M , >M & <=N , >N # let M <-.8 N <- 1.2 # In my real example I have many more of these cutoff-points # What I did is: out <- matrix(NA,1,3)

Dear R-list, i have 6 different sets of samples. Each sample has about 5000 observations, with each observation comprised of 150 baseline covariates (X), 125 of which are dichotomous. Roughly 20% of the observations in each sample are "treatment" and the rest are "control" units. i am doing propensity score matching, i have already estimated propensity scores(predicted

Hi, I would like to add a prefix to colnames in a matrix but I can't get the prefix option in colnames to work. What am I doing wrong? > X<-matrix(NA,3,4) > colnames(X)<-c("test","test","test","test") > colnames(X)<-colnames(X,prefix="PREFIX.") > X test test test test [1,] NA NA NA NA [2,] NA NA NA NA

Ilija Barukcic, Jever, Germany. GMT + 1h. Barukcic at t-online.de Hallo, you are receiving this message from me as someone who potentially has bought the first Edition of my book Causality. New statistical methods. ISBN 3 - 8 3 3 4 - 3 6 4 5 - X. Your email address was obtained either through my friends, your publication, a colleague or via correspondence with you. I am pleased to

Hi, i have been running R 2.11 compiles with the goto blas libraries successfully for about a year. I am running Ubuntu 10.10 64 bit. However, recently I have been periodically getting a segfault when running R (and numpy -- also compiled against goto blas). When I check the dmesg error log I get the following message R[28792]: segfault at 7f4183ed4000 ip 00007f4196f05f03 sp 00007f4192d55d20

Hi. I'm trying to install a new package. I'm a relatively long-time (though not advanced) R user and have never seen this error before. For the first example, I tried a few different CRAN mirrors. In the second example, the file does exist; I downloaded it from the CRAN website for the package and pasted in the name exactly a few different times to make sure it was right. >

Hi, With two data sets, one complete and another one partial, I would like to merge them and keep the unmatched lines. The problem is that merge() dosen't keep the unmatched lines. Is there another function that I could use to merge the data frames. Example: completedf <- expand.grid(alpha=letters[1:3],beta=1:3) partdf <- data.frame( alpha= c('a','a','c'),

Hi, I am running Linux (64 bit) R, compiled under gnu compilers. I am looking to compile a Fortran program with the ifort (the intel compiler) and then be able to import that as a library into my version of R. I am using the flags FC=ifort and SHLIB_FCLD=ifort, however ifort does not seem to recognize the libR.so file. Does any one have experience with this/know of a a way to load an ifort

Hi, I am trying to use the merge command in the data.tables package. However, when I run the command I am not sure if it is running the merge command from the base package or the merge command from data.tables. When I run "methods(generic.function="merge")' it informs me that 'merge.data.table" is "non-visible". I am just trying to run the merge command on two