similar to: prop.test vs hand calculated confidence interval

Hello, I'm a first semester statistics student<http://r.789695.n4.nabble.com/Question-regarding-error-quot-x-and-y-lengths-differ-quot-td4665773.html#>and I am using R for roughly the third time ever. I am following a tutorial and yet I still get the error "x and y lengths differ." I am very new to this program, and I have searched for solutions, but because I do not understand

Dear R-users: The following code produces two cones in two panels. What I would like to have is to have them in one, and to meet in the origin. Does anyone have any good ideas how to do this? Thanks for your help Jaakko library(lattice) A<-matrix(ncol=2, nrow=64) for(i in 0:63) { A[i+1,1]<-sin(i/10) A[i+1,2]<-cos(i/10) }

Hi, I've two fitted models, one binomial model with presence-absence data that predicts probability of presence and one gaussian model (normal or log-normal abundances). 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

Hello , i am working in the functional time series using the multivariate time series data(hourly time series data). Sir? i am using FAR model more than one order for which no statistical package is available in R, so for this i convert my data into functional form and obtained the functional principle component and from those FPCA i extract their corresponding? FPCscores. Know i use the VAR model

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

This came up on r-help but indicates a bug. dnbinom(x,n,p) calls dbinom_raw(n-1,...) which returns 0 for n<1. -thomas ---------- Forwarded message ---------- Date: Thu, 08 Feb 2001 17:10:23 +0000 From: Yudi Pawitan <yudi@stat.ucc.ie> To: Mark Myatt <mark@myatt.demon.co.uk> Cc: R-Help <r-help@stat.math.ethz.ch> Subject: Re: [R] Goodness of fit to Poisson / NegBinomial

Hi, there I am a graduate student new to coding in S who is hitting a bit of a wall at present using an "optim" function. I am running into some troubles, and was hoping someone might be able to recognize where I am going wrong. As background: I have constructed a loop that carries out a 365-day calculation for a mass-balance model. Basically, the model depends on 2 variables (p,

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:

I still don't think the exp(lp)/(1+exp(lp)) gonna work. Since this is conditional logit model, 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

Hello , i am working in the functional time series using themultivariate time series data(hourly time series data). Sir? i am usingFAR model more than one order for which no statistical package is available inR, so for this i convert my data into functional form and obtained thefunctional principle component and from those FPCA i extract theircorresponding? FPCscores. Know i use the VAR model on

Hello R-list, I'm attempting to migrate from Stata to R for my complex survey work. It has been straight-forward so far except for the following problem: I have some code below, but first I'll describe the problem. When I compute predicted logits from a logistic regression, the standard errors of the predicted logits are way off (but the predicted logits are fine). Furthermore, the

Thank you so much for your help. The function "dbinom" seems to work very well. However, I'm a bit lost with the "dnorm" function. Apparently, I have to compute the mean "mu" and the standard deviation "sd" but what does it mean exactly? I only have a vector of predicted response and a vector of observed response that I would like to compare! What

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

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

UNSUBSCRIBE AT THE BOTTOM --------------------------------------------------------------- Dear Subscriber/Member, You don't need to spend hundreds of dollars on a DVD burner to backup your DVD's! DVD Professional is the most technologically advanced method of DVD reproduction ever available and it's the only system you will need to create backups of your DVD's that will play in

Hi there I thought this would be of particular interest to people using 'optim' functions and perhaps people involved with R development. I've been beaten down by R trying to get it to perform an optimization on a mass-balance model. I've written the same program in excel, and using the 'solver' function, it comes up with an answer for my variables (p, ACT, which

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, >>

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

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

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,