similar to: Likelihood between observed and predicted response

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

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 there is a difference. D. -- View this message in context:

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

Dear all, I am encountering some issues with my data and need some help. I am trying to run glm analysis with a presence/absence variable as response variable and several explanatory variable (time, location, presence/absence data, abundance data). First I tried to use the glm() function, however I was having 2 warnings concerning glm.fit () : # 1: glm.fit: algorithm did not converge # 2:

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

Hi all I'm looking for some help on winbind/idmap for a new host The debian version is new on this host : debian 10 buster so samba is > root at homedir10:~# samba --version > Version 4.9.5-Debian We want to use security=ads so we join this host to the domain No problem for windows clients : they can mount shares that are accessible to their primary unix group and secondary unix

Hi, I am using the densCols to draw a high volume scatter plot. Instead of using the default local density, I would like to take log of the local density and then map them to the colors. I could not figure out how to do that. For example: plot(x,y,col=densCols(x,y,"log")) ? Any help would be appreciated! Jiaxiu [[alternative HTML version deleted]]

Dear all, I am organizing a set of specific R code as package (to ease the documentation and deployment of it to users). Before doing so, I would like to know if there are written coding rules for R (with functions, objects naming convention for example). Thanks a lot, Thomas -- ------------------------------------------------------------- Thomas LOUBRIEU IFREMER IDM/ISI BP70 29280

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

I am trying to fit regression models type II with R, but it seems to me that most (all) of the linear model functions are for type I regressions. Does anyone knows whether type II regressions functions exist in R. Benjamin -- Benjamin Planque IFREMER Laboratoire d'Ecologie Halieutique BP 21105 44311 Nantes Cedex 03 France Tel: +33 (0)2 40 37 41 17 Fax: +33 (0)2 40 37 40 75 e-mail:

Hello, I don't manage to see if you have already focussed on this point in some previous messages so I post my question: I have a little problem with the S4 style of programming. I tried to formalize my question: please consider the following example that you can run I think: #------------------ setClass("my.class1", representation( my.list =

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

Hi, On a linux (suse 11p1 64 bits) and intel compiler (11.1 20100806), for R 2.14.1 I ran : ./configure make The config.log file is available at ftp://ftp.ifremer.fr/ifremer/divers_permanents/config.log Binary file R is created, all seems Ok. However : make check fails with the message : make[1]: Entering directory `/tmp/R/R-2.14.1/tests' make[2]: Entering directory

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