search for: phat

Hi, This code: n=40 x=17 phat=x/n SE=sqrt(phat*(1-phat)/n) zstar=qnorm(0.995) E=zstar*SE phat+c(-E,E) Gives this result: [1] 0.2236668 0.6263332 The TI Graphing calculator gives the same result. Whereas this test: prop.test(x,n,conf.level=0.99,correct=FALSE) Give this result: 0.2489036 0.6224374 I'm wondering why t...

...). I would like to evaluate these models not on their capability of adjustment but on their capability of prediction by calculating the (log)likelihood between predicted and observed values for each type of model. I found the following formula for Bernouilli model : -2 log lik = -2 sum (y*log phat + (1-y)*log(1-phat) ), with "phat" is the probaility (between 0 and 1) and "y" is the observed values (0 or 1). 1) Is anybody can tell me if this formula is statistically true? 2) Can someone tell me what is the formula of the likelihood between observed and predicted values...

Hi all, I'd like to do cross-validation on lm and get the resulting lift curve/table (or, alternatively, the estimates on 100% of my data with which I can get lift). If such a thing doesn't exist, could it be derived using cv.lm, or would we need to start from scratch? Thanks! -- Eric Siegel, Ph.D. President Prediction Impact, Inc. Predictive Analytics World Conference More info:

...***? 0.001 ?**? 0.01 ?*? 0.05 ?.? 0.1 ? ? 1 (Dispersion parameter for binomial family taken to be 1) Null deviance: 145934 on 105857 degrees of freedom Residual deviance: 145750 on 105856 degrees of freedom AIC: 145754 Number of Fisher Scoring iterations: 3 > > #PREDICTED SEs > phat.l.se.logit.nosvy <- predict(model.l.nosvy,se=TRUE) > as.matrix(table(phat.l.se.logit.nosvy$se.fit)) [,1] 0.00632408017609573 14456 0.00633130215261306 15188 0.00741988836010757 12896 0.00743834214717549 10392 0.00923404822144662 13207 0.00925875968615561 15864 0.011429466...

Hi r-users,   I would like to use prop.test code and I also calculate manually to test the proportions for 2 groups.  The problem is the answer for the p-value calculated manually are different from prop.test.  Here are the results:   ## Manually   z value: z= (phat-p)/sqrt(pq/n) = (.084-.081)/sqrt(.081(1-.081)/691)=0.289, pvalue=0.7718   ## Using prop.test code > low <- c(56,58) > tot <- c(691,691) > prop.test(low, tot, p = c(56/691,58/691), alternative = "two.sided", conf.level = 0.95, correct = TRUE)         2-sample test for give...

...dels not on their capability of > adjustment but on their capability of prediction by calculating the > (log)likelihood between predicted and observed values for each type of > model. > > I found the following formula for Bernouilli model : > > -2 log lik = -2 sum (y*log phat + (1-y)*log(1-phat) ), with "phat" is > the probaility (between 0 and 1) and "y" is the observed values (0 or 1). > > 1) Is anybody can tell me if this formula is statistically true? This looks correct. > 2) Can someone tell me what is the formula of the lik...

...om_ function(x,n,p){ gamma(x+n)/gamma(n)/gamma(x+1)*p^n * (1-p)^x } fn_ function(p){ r_ p[1] th_ p[2] ll _ -sum(x[1:max(kk)]*log(dnbinom(kk[-length(kk)],r,th))) - x[length(kk)]*log(1-pnbinom((max(kk)-1),r,th)) return(ll) } a_ nlm(fn,p=c(.9,.5),hes=T) est_ a$est phat_ c(dnbinom(kk[-length(kk)],est[1],est[2]), 1-pnbinom((max(kk)-1),est[1],est[2])) E_ n*phat chi2_ sum((x-E)^2/E) cat('Chi2 goodness-of-fit = ',chi2,'\n') cat('P-value=',1-pchisq(chi2,df=length(kk)-1-2),'\n') print(cbind(O=x,E=round(E,1),err=round((x-E)^2/E,1)...

...functional form and obtained the functional principle component and from those FPCA i extract their corresponding? FPCscores. Know i use the VAR model on those FPCscores for the forecasting of each 24 hours through the VAR model, but the VAR give me the forecasted value for all 23hours? when i put phat=23, but whenever i put phat=24 i.e want to predict each 24 hours its give the results in the form of NA. the code is given below fdata<- function(mat){? nb = 27 # number of basis functions for the data? fbf = create.fourier.basis(rangeval=c(0,1), nbasis=nb) # basis for data? args=seq(0,1,length=...

...gt; EA <- 0.002938 > EB <- -0.2 > EQ <- 0.066 > a <- 0.8 > > #Specifying sex-specific parameters > > GSI<- NULL > > if (Sex==1) GSI<-0.05 else + if (Sex==2) GSI<-0.17 > > # Define margin of error functions > #merror <- function(phat,M,alpha) # (1-alpha)*100% merror for a proportion > # { > # z <- qnorm(1-alpha/2) > # merror <- z * sqrt(phat*(1-phat)/M) # M is (Monte Carlo) sample size > # merror > # } > > #Bring in temp file > > temper <- scan("temp.dat", n...

Hello, I'm fairly new to R so please excuse me if I am asking something obvious. I have looked in the FAQ, Introduction, and help pages, and searched the archives, but I don't know much about graphics yet. I'm running Red Hat Linux 2.14.18 on a machine blessed with dual 1.5 Xeon processors and 3.7GB of RAM. I have a very large dataset with 27 variables, and in exploring the data

...;s the only system you will need to create backups of your DVD's that will play in your home DVD player. Get it Today and you get: ------------------------------ + 70% off the Regular Price + Full Technical Support + Free Updates For more information, please visit their Web Site: http://www.phat-x.com/20025/ ---- [ REMOVE ] ------------------------------------------------- You are receiving this email as a subscriber to our mailing list. To remove yourself from this and related email lists click here: http://www.phat-x.com/cgi-bin/remove.cgi ---- [ NOTE ] -----------------------------...

...del, while this formula is only used in unconditional ones. By using this, one neglects the information based on stratum. Though I don't know how to solve it to. I am also working on a project on this and I do hope there's someone explaining this problem. Will that be a possibility that the phat can never be estimated as we never know the individual intercept? Disclaimer: This email (including any attachments) is for the use of the intended recipient only and may contain confidential information and/or copyright material. If you are not the intended recipient, please notify the sender...

...into functional form and obtained thefunctional principle component and from those FPCA i extract theircorresponding? FPCscores. Know i use the VAR model on those FPCscores forthe forecasting of each 24 hours through the VAR model, but the VAR give me theforecasted value for all 23hours? when i put phat=23, but whenever i putphat=24 i.e want to predict each 24 hours its give the results in the form ofNA. the code is given below ? fdata<- function(mat){ ? nb = 27 # number of basis functions for the data ? fbf = create.fourier.basis(rangeval=c(0,1), nbasis=nb) #basis for data ? args=seq(0,1,...

I am yet another person who is having a problem with Win98 clients needing to make several attempts before successfully logging in to our Samba box. I get the "broken pipe" message in the smb.log just like everyone else. This did not seem to be a problem until we got a phat broadband connection. All the Windows boxes are set to point to our ISP's DNS. I am not running DNS on the local network at all. Not even a cacheing DNS. I don't think this would matter anyway, as several people have pointed out that the issue is with reverse DNS, and this wasn't even...

Hi Aki, I'm not sure if this bug was worked on?? Sadly I ran into it too recently, on version 2.3.4.1. Regards, Aragon On 2018/02/07 13:23, Aki Tuomi wrote: > Hi! > > Thank you for reporting this issue. We'll look into it. Maybe you can > use sqlite3 instead as workaround? > > Aki > > > On 07.02.2018 14:22, Marco Giunta wrote: >> Hi, >>

...A<- 0.0253 UB<- 0.58 UG<- -0.299 #Mass balance model parameters EA <- 0.002938 EB <- -0.2 EQ <- 0.066 a <- 0.8 #Specifying sex-specific parameters GSI<- NULL if (Sex==1) GSI<-0.05 else if (Sex==2) GSI<-0.17 # Define margin of error functions #merror <- function(phat,M,alpha) # (1-alpha)*100% merror for a proportion # { # z <- qnorm(1-alpha/2) # merror <- z * sqrt(phat*(1-phat)/M) # M is (Monte Carlo) sample size # merror # } #Bring in temp file temper <- scan("temp.dat", na.strings = ".", list(Day=0, jday=0, Temp=...

It's not very high priority I'm afraid. Did you try using sqlite instead? Aki > On 13/12/2019 17:40 Aragon Gouveia <aragon at phat.za.net> wrote: > > > Hi Aki, > > I'm not sure if this bug was worked on?? Sadly I ran into it too > recently, on version 2.3.4.1. > > > Regards, > > Aragon > > > > On 2018/02/07 13:23, Aki Tuomi wrote: > > Hi! > > > >...

This is a formal announcement of Lguest64. Most are aware of the little puppies (lguest32, or simply lguest, or in some circles "rustyvisor"). But this time the puppies ate a bit too much. No more lean and mean puppies, now we got big fat lazy ones. Running on the hardware that's too lazy to do full virutalization. Yes, lguest now runs on x86_64! As you know, puppies are young,

The latest rsync (2.6.5) seems to invoke the transport specified by $RSYNC_RSH differently. I have this set to point at an ssh wrapper script (which is now I noticed), and an strace shows this: [archives/silva]#root@jazzbo*> strace -e trace=process -f /opt/bin/rsync -avHP /mnt/phat/archives/silva/. root@ali:/a/ali/home/............. [...snip...] execve("/u/cameron/scripts/rsync_rsh", ["rsync_rsh", "ali", "-l", "root", "rsync", "--server", "-vlHogDtpr", "--partial", ".", &qu...

This is a formal announcement of Lguest64. Most are aware of the little puppies (lguest32, or simply lguest, or in some circles "rustyvisor"). But this time the puppies ate a bit too much. No more lean and mean puppies, now we got big fat lazy ones. Running on the hardware that's too lazy to do full virutalization. Yes, lguest now runs on x86_64! As you know, puppies are young,