similar to: number of matches when using Match()

Hello! I need to compare 2 datasets whether they come from the same distribution. I use function ks.boot{Matching}. And what is the confidence level of the p-value, returned by ks.boot function? The code is: set=read.table("http://stella.sai.msu.ru:8080/~gala/data/testsets.csv", header=T,sep=',') set1=set[!is.na(set$set1),'set1']

Does anyone else find that using the Var.calc option (for heteroscedasticity consistent std. errors) in Match() (from the Matching library) slows down computation of the matching estimator by a lot? I don't really understand why when I use this option it slows down so much, but for me it does significantly. I want to use the heteroscedasticity consistent std. errors in my project, but as long

Dear R Users, I am a novice learner of R software. I am working with Genetic Matching - GenMatch(), but I am getting an Error message as follows:  Error in GenMatch(Tr = Tr, X = X.binarynp, BalanceMatrix = BalanceMatrix.binarynp,  :    GenMatch(): input includes NAs Could you please suggest me correcting the above problem? My GenMatch command is, > gen1 <- GenMatch(Tr = Tr, X =

Perhaps someone will have a solution to my more general problem, but here is the specific one: I used the round() function to round some estimates to 3 decimal places. I then sent put the rounded estimates in a matrix and used latex() to make a LaTeX table from them. However, in my table, there are estimtes which only have 2 decimal places. I assume that the third decimal place in these numbers

hello, I have been trying to use the predict function in the vars package to forecast from a seasonal VAR model. The following code sample illustrates what I am trying to do and the error that I get when trying to do it. I run the following code, that results in the following error: data(Canada) endoC <- Canada[1:72,1:3] exoC <- Canada[1:72,4] var.2c <- VAR(endoC, p = 2,

Hi I'm getting an error message with the roots() function in the vars package. Even the example in the help file comes up with an error: > data(Canada) > var.2c <- VAR(Canada, p = 2, type = "const") > roots(var.2c) Error in UseMethod("roots") : no applicable method for 'roots' applied to an object of class "varest" The error is odd, for

Thanks to the help of many on this list, I am now an R user and have been able to write some functioning code to do matching estimation. I have two for loops (i in 1:3, and j in 0:2). Within the loops, I had been creating matrices of relevant estimation coefficents in order to make lots of LaTeX tables. Well, now I want to be able to combine the results of many different estimations from within

Dear All I am using rgenoud to solve the following maximization problem: myfunc <- function(x) { x1 <- x[1] x2 <- x[2] if (x1^2+x2^2 > 1) return(-9999999) else x1+x2 } genoud(myfunc, nvars=2, Domains=rbind(c(0,1),c(0,1)),max=TRUE,boundary.enforcement=2,solution.tolerance=0.000001) How can one increase the precision of the solution $par [1] 0.7072442 0.7069694 ? I

Dear All I am a biginner of R. I have difficulty with reading the code of a method. I am using the vars package to estimate a VAR model and I want to view the code of "predict" method for objects with class attribute "varest". I thougt I could just type the name "predict" without anything to display the code of the method as I often do with generic function.

Hello! I need to purchase a new box, which I would like to optimize for good R performance. For the record, I will run Fedora Core 5 as and OS, and I wanted to know if anyone has experience with how the following affects R performance: - Is there a big advantage to having a 64-bit CPU over having a 32-bit? - Does an Opteron offer any advantages over an Athlon, and if yes, does it justify an

Dear all Off topic to some extent but Ted Harding has just called me to say he is (a) seeing large delays in the list messages when it is sent internally from hypatia.math.ethz.ch. I get only a delay of an hour which is reasonable in the scrap below but Ted reports delays of the order 24 hours, (b) messages he has sent to the list don't appear on it for a long time. Is anyone else

Hi guys. I have a dataset with 4 columns. In the first and second column I have the same qualitative variable referred to different teams of people. There are 10 teams in total and they compete against each other to perform a certain task whose result is stored in the third column for the team recorded in the first column, and in the fourth column for the team in the second column. For example,

I am trying to determine whether two samples are identical or not. I'm aware that somebody can use the Kolmogorov-Smirnoff test to compare empirical distributions, but since my samples have ties I'm not sure if I'm getting the right p-values for the comparison. Can the Kolmogorov-Smirnoff test be adjusted for the case when ties exists and are there any functions that already

I have been reviewing the arch test of vars package that is based on Engle's paper. R-metrics has this as a wants/todo item. I need an arch test and thought I might accomplish two things at once or at least try. The arch test requires a varest object and I am trying to write one that will use R-metrics arima, arch, Garch objects, or at least a vector. Also the arch function has the

Hello, I am running R version 2.4.1 on Debian. I am using some R functions which produce several graphics displayed successively with an interactive prompt between each graph. For example, the plot.varest (vars) function produces such output (one graph per time series on which the VAR model is fitted). Plots on GARCH models also produce such outputs. Moreover, the frontend I am using for R

I'm interested in doing multinomial logistic regression with equality constraints on some of the parameter values. For example, with categorical outcomes Y_1 (baseline), Y_2, and Y_3, and covariates X_1 and X_2, I might want to impose the equality constraint that \beta_{2,1} = \beta_{3,2} that is, that the effect of X_1 on the logit of Y_2 is the same as the effect of X_2 on the

Hi R users, "genoud" function of "rgenoud" package will optimize my function. If opt = genoud(fn,2,max=TRUE,starting.value=c(1,10),........) opt$value will give the optimized value of the function, "fn". My problem is from the same opt, can I get the value of the function at the initial parameter values? I need the initial value of the function for

I am trying to use the optimizing function genoud() with the snow package on a couple of i686 machines running Redhat Linux WS4 . I don't know anything about PVM or MPI, so I just followed the directions in snow and rgenoud for the simplest method and started a socket cluster. My function fn for genoud involves updating an input file for a separate numerical model with the latest parameter

Dear everyone, I am trying to run rgenoud on several chips simultaneusly. I used the instructions provided on Jasjeet Sekhon's Homepage (http://sekhon.berkeley.edu/rgenoud/multiple_cpus.html). However, I have the newer version of R (R 2.12) installed - for a 64-bit machine. So, when I tried to install the special version of "snow" from a zip file provided by Jasjeet on his page, R

Hi, If i have smbpasswd with: user1:1416:803A317873C24BBDAAD3B435B51404EE:2DF2CB1538FE718DE034707A521AA893: [U ]:LCT-1221415636: and i do: pdbedit -i smbpasswd:file And SID for domain SERVER is: S-1-5-21-1454471165-2146950999-672003340 How does samba generate the rid part of that users sid? [root@server samba]# pdbedit -L -v user1 User SID: