search for: probablility

...t <h$density> holds the values # plotted in the probability histogram. If that's the # case, I'd expect that the sum of h$density for a histogram where freq=F would equal 1.0 ... However: > sum(h$density) #returns the value : #[1] 10 I would really like to plot values in the probablility histogram that sum to 1, not 10. Is there a switch in the hist(x) command that I'm missing? Or, is there another function that can do this? Thanks in advance, jlb -- ************************************ Joseph P. LeBouton Forest Ecology PhD Candidate Department of Forestry Michigan State...

..., data=age) svymean(~y, age.des3) ## gives mean 26.296, SE 0.9862 age.des4<-svydesign(ids=~id, strata=~stratum, data=age) svymean(~y, age.des4) ## gives mean 25.875, SE 0.9437 age.des3 is the only estimator I am able to compute per hand correctly. It is stratified random sampling with inverse probablility weighting with weight= nh/Nh ## sample size/ stratum size. Basically, I thought the option weight=~Nh as well as weight=~I(nh/Nh) would result in the same number, but it does not. I am reading Thompson(02), Cochran(77) and of course Lumley on his Survey package, but I can't find my mistake...

..., "l"), each=6, 5), species=rep(c("sl", "hn"), each=3, 10), damage=rep(c(0,.1,.25), 20) ) bm.surv$survival<-ifelse(bm.surv$light=="h", rbinom(60, 4, .9), rbinom(60, 4, .6)) # difference in probablility should ensure a light effect bm.surv$death<-4-bm.surv$survival # fitting the model m1<-lmer(cbind(survival, death)~light+species+damage+(1|house), data=bm.surv, family="quasibinomial") summary(m1) # suggests that light is very significant Generalized linear mixed model f...

...ulating this probability requires knowing one of the 3 baseline functions-the survivor function or the hazard function or the cumulative hazard function-estimated from the Cox model. I would like to know : 1) How do I extract the correct baseline hazard so that I can predict the individual probablility of failure at 2000 days, given a particular patient covariate pattern (in this case, a particular patient's age? Is the baseline hazard for the average patient covariates the ideal one? 2) Is there a better way of directly extracting these probabilities? Brant Inman [[alterna...

...up would have implicitly avoid the glare condition through signalling the intent-to-use to the resource to the remote end (a-la E&M wink start). So what seems to be happening in my case is my CPE is set for a 'Descending' b channel sequence and Asterisk is set for same, thus the high probablility of glare. I can change the direction of my CPE which will solve the problem in this case as I own both ends of the link, but in a telco situation where Asterisk is the CPE how would I set it to be be the opposite of the serving central office? PRI b-channel sequence ordering in Asterisk can't b...

Here is a patch that fixes distortions that appear during scenes with no motion. Tbe problem was that UpdateFrame() was not being called when a MotionScore of 0 was computed. Since UpdateFrame() was not called the buffer returned by theora_encode_packetout() would be the same buffer that was output for the last frame that had a non-zero MotionScore. This is obviously incorrect behavior. I've

...f the nth most frequent word, we have, roughly p_n = 0.1/n . Zipf has pointed out that this type of formula, p_n = k/n, gives a rather good approximation to the word probabilities in many different languages. The formula [above] clearly can not hold indefinitely since the total probablility SUM p_n must be unity, while SUM_1^inf 0.1/n is infinite. If we assume (in the absence of any better estimate) that the formula p_n = 0.1/n holds out to the n at which the total probability is unity, and that p_n = 0 for larger n; we find that the critical n is the word rank 8,727. The entrop...

Hi, I have some doubts on sip.conf. 1) Can I have two or more SIP phones acting as extensions in one Asterisk box, and at the same time, registered to a SIP proxy, say Free World Dialup? If yes, how? 2) Why we need a section in the sip.conf for the proxy, say, Free World Dialup's fwd.pulver.com? In the case of 1), how to assign the value to section [fwd.pulver.com], since there are more

R-2.10.1.tar.gz was built a short while ago. This is a maintenance release and fixes a number of mostly minor issues. See the full list of changes below. You can get it from http://cran.r-project.org/src/base/R-2/R-2.10.1.tar.gz or wait for it to be mirrored at a CRAN site nearer to you. Binaries for various platforms will appear in due course. For the R Core Team Peter

R-2.10.1.tar.gz was built a short while ago. This is a maintenance release and fixes a number of mostly minor issues. See the full list of changes below. You can get it from http://cran.r-project.org/src/base/R-2/R-2.10.1.tar.gz or wait for it to be mirrored at a CRAN site nearer to you. Binaries for various platforms will appear in due course. For the R Core Team Peter

Currently rsync has a bit of a problem with very large files. Dynamic block sizes were introduced to try handle that automatically if the user didn't specify a block size. Unfortunately that isn't enough and the block size would need to grow faster than the file. Besides, overly large block sizes mean large amounts of data need to be copied even for small changes. The maths indicate