search for: emper

Hello all, I would like to plot the emperical CDF, normal CDF and pareto CDF in the same graph and I amusing the following codes. "z" is a vector and I just need the part when z between 1.6 and 3. plot(ecdf(z), do.points=FALSE, verticals=TRUE, xlim=c(1.6,3),ylim=c(1-sum(z>1.6)/length(z), 1)) x <- seq(1.6, 3, 0.1) lines(x,...

Dear all, I am trying to get the emperical bayes estimates together with their standard errors out of lme4. Up to now I have used MLwiN to get these estimates. I have fitted the following - very simple - model, just to find out how this works. test<-lmer(y~(1|subject),data,REML=F) ranef(test,postVar=T) str(ranef(test,postVar=T) I...

Dear friends. In Carlin and Louis "Bayes and emperical Bayes methods.." 1996 the classical example of 12 independent tosses of a fair coin producing 9 heads and 3 tails is given. If the situation is seen as a fixed sample of 12, a binomial lieklihood is used, and Carlin et al reports a probability of 0.075. Using sum(dbinom(9:12,12,.5)) I...

Hi all, I would like to plot part of the emperical CDF. Suppose the variable is x, I just need the part when x>1,therefore, I am using the following codes. tail <- x>1 plot(ecdf(x[tail]), do.points=FALSE, verticals=TRUE) The "x" value starts from 1, but the yaxs still begins from 0, not the corresponding value when "x...

Hi, I would like to use the following codes to plot the CDF for pareto distribution. Before doing this, I have plot the emperical one. x <- seq(1.6, 3, 0.1) lines(x,pgpd(x, 1.544,0.4477557,), col="red") Could anyone give me some advice whether the above codes are correct? Many thanks. -- View this message in context: http://www.nabble.com/CDF-for-pareto-distribution-tf4061253.html#a11538272 Sent from t...

Hi all. Belove is the example from the cluster-help page wtih the output. I simply cannot figure out how to relate the estimate and robust Std. Err to the p-value. I am aware this a marginal model applying the sandwich estimator using (here I guess) an emperical (unstructered/exchangeable?) ICC. Shouldent it be, at least to some extend, comparable to the robust z-test, for rx : 2*pnorm(-0.239/0.0816)=0.0034 ? Any help/hints are appreciated. Steen and R 2.2.1, survival 2.21 on win XP. library(survival) data(rats) marginal.model <- survreg(Surv(ti...

Hi Ewan, I''m almost ready to integrate some PPC-specific stuff into xend, and I was wondering if you had a plan for how that should work. First example: the device tree data structure we talked about a few weeks ago. We will need to pass the config data to PPC code, probably in XendDomainInfo.initDomain(), and then pass the resulting data structure into libxc''s xc_linux_load()

Dear R experts, density estimation on a 2 dimensional bounded domain --------------------------------------------------------------------- I am currently trying to estimate the probability density (PD) of cancers within the breast using the sm library with the routine sm.density Of course a practical PD must be limited by the curve of the breast outline. I don't have a clue after perusing

Sorry, for the re-post, but I''m new to the mailing list and wanted to bring back up and old topic I saw in the archives. http://rubyforge.org/pipermail/mongrel-users/2008-February/004991.html I think a patch to delay garbage collection and run it later is pretty important for high performance web applications. I do understand the trade-offs of having explicit vs. implicit garbage

Hi folks. I just converted my site from UW+mbox to dovecot+maildir (stable-test59, 2005/02/25). For the most part, everything is great. The server seems faster and less memory intense than the UW server (this part has made me very happy). And it's definitely more user friendly to me than Courier. There is one thing however... A few of my users are using the web based mail reader

1. HFSC have 4 curve such sc, rc, ls, ul and 1.1 In leaf class can specify rc for guarantee service (bandwidth and delay) and If want to sharing fairness exceess service, we must specify ls and ul curve too (ls curve with paramater m2 specify at lease sharing bandwidth in that class will receive and ul curve mean maximum bandwidth in that class will receive) so i''m doubt .. about if i

We are running samba 2.2.3a ( for about 3000 users ) on HP/UX 11.11 fileservers. Our customers are complaining about slow response on the opening/closing of files in office/excel etc. In the samba log file we see the following messages : [2002/10/29 09:03:27, 0] ../source/smbd/oplock.c:(761) oplock_break: receive_smb timed out after 30 seconds. oplock_break failed for file My

Hi, ok I know zfs-fuse is still incomplete and performance has not been considered, but still, before I''m going to use it for my /home I wanted a rough estimate. Another benchmark already asserted that zfs by itself, on Solaris, is a very fast beast (http://cmynhier.blogspot.com/2006/05/zfs-benchmarking.html), faster than ext3 on Linux and also faster than UFS. Interestingly UFS seems