search for: sekhon

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

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 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 did not allow me to insta...

...ters and it can handle datasets in which the number of categories varies by observation. DESCRIPTION: Package: multinomRob Version: 1.0 Date: 2004/02/18 Title: Robust Estimation of Overdispersed Multinomial Regression Models Author: Walter R. Mebane, Jr. <wrm1 at cornell.edu>, Jasjeet Singh Sekhon <jasjeet_sekhon at harvard.edu> Maintainer: Jasjeet Singh Sekhon <jasjeet_sekhon at harvard.edu> Description: overdispersed multinomial regression using robust (LQD and tanh) estimation Depends: R (>= 1.7.0), rgenoud (>= 1.22), MASS (>= 7.1-8), mvtnorm (>= 0.6-3) License: GP...

...ters and it can handle datasets in which the number of categories varies by observation. DESCRIPTION: Package: multinomRob Version: 1.0 Date: 2004/02/18 Title: Robust Estimation of Overdispersed Multinomial Regression Models Author: Walter R. Mebane, Jr. <wrm1 at cornell.edu>, Jasjeet Singh Sekhon <jasjeet_sekhon at harvard.edu> Maintainer: Jasjeet Singh Sekhon <jasjeet_sekhon at harvard.edu> Description: overdispersed multinomial regression using robust (LQD and tanh) estimation Depends: R (>= 1.7.0), rgenoud (>= 1.22), MASS (>= 7.1-8), mvtnorm (>= 0.6-3) License: GP...

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']

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

To anyone who uses the Match() function in the Matching library... How do you go about deciding how many matches you will use? With my data, my standard errors generally get smaller if I use more matches. Speaking of standard errors, when correcting for heteroscedasticity, how many matches do you use (this is the Var.cal option). It seems to me that it might make sense to use the same number

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

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

...offer any advice on what to do about this? I am running 3.0.24 on Debian Etch and 3.0.28a on Gentoo, for which those are the latest stable versions packaged for the systems. I have tried 3.0.32 and the problem seems to occur there too. Is this a bug that has crept in after 3.0.24? -h -- Hari Sekhon

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:

I see something new and unexpected here. > update.packages() trying URL `http://cran.r-project.org/src/contrib/PACKAGES' Content type `text/plain; charset=iso-8859-1' length 163467 bytes opened URL .......... .......... .......... .......... .......... .......... .......... .......... .......... .......... .......... .......... .......... .......... .......... ......... downloaded

Hi, I am trying the following: library(ismev) library(evd) fit <- gev.fit(x,show=FALSE) ks.test(x,pgev,fit$mle[1],fit$mle[2],fit$mle[3]) but I am getting: Warning message: cannot compute correct p-values with ties in: ks.test(x, pgev, fit$mle[1], fit$mle[2], fit$mle[3]) where x is: [1] 239 38 1 43 22 1 5 9 15 6 1 9 156 25 3 100 6 [18] 5 100

...ow and rgenoud packages. The function I want to run in task-parallel fashion across multiple machines is one that pre- and post-processes data and runs an external model code. My problem is that external file I/O is happening only in the master node and not in the slaves. I have followed Jasjeet Sekhon's suggestion to test the cluster setup, and that is fine: > library(snow) > > #pick two machines > cl <- makeCluster(c("moab","escalante")) > > clusterCall(cl, sin, 2) > The output should be: > > clusterCall(cl, sin, 2) > [[1]] > [1]...

...t; y <- rnorm(100000) > z <- rbinom(100000,1,0.002) > logit.reg <- glm(z~x+y,family=binomial(link='logit')) > match <- > Match(Y=NULL,Tr=z,X=logit.reg$fitted,version='fast',ties=TRUE,M=1,distance.tolerance=1e-5) According to the function definition (http://sekhon.berkeley.edu/matching/Match.html): "distance.tolerance: This is a scalar which is used to determine if distances between two observations are different from zero. Values less than distance.tolerance are deemed to be equal to zero. This option can be used to perform a type of optimal matching&...

Hello list, I am dealing with a noisy function (gradient,hessian not available) with simple boundary constraints (x_i>0). I've tried constrOptim() using nelder mead to minimize it but it is way too slow and the returned results are not satisfying. simulated annealing is so hard to tune and it always crashes R program in my case. I wonder if there are any packages or functions can do

I'm looking to buy a new desktop which will primarily be used for analyses of large datasets (100s of MB). I've seen postings from several years back re the 'optimal' platform for running R, but nothing more recently. Specifically, I want to know: 1) if I run R under Windows, does having a dual-processor machine help speed things up? And 2) is it still true that R performs