similar to: Likelihood between observed and predicted response

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

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:

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

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

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

On 18/11/2020 10:42, Gregory ROCHER via samba wrote: > Hi all > > I'm looking for some help on winbind/idmap for a new host > > > 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 group(s) You are using AD now, so the primary group

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

Hi, In rw-FAQ, there is a paragraph about optimizing R by using a fast ATLAS library in place of the standard BLAS. It says: "BLAS support is supplied by the single DLL `R_HOME\bin\Rblas.dll', and you can add a fast BLAS just by replacing that." Do I need to recompile the whole R with the option `USE_ATLAS = YES', or is it possible to compile only Rblas.dll on each separate

Hi, Is it any way to convert .Rd files to CHTML Splus 6 files under Windows? Thanks, Philippe Grosjean ...........]<(({?<...............<?}))><............................... ( ( ( ( ( ) ) ) ) ) Philippe Grosjean ( ( ( ( ( ) ) ) ) ) IFREMER Nantes - DEL/AO ( ( ( ( ( rue de l'Ile d'Yeu, BP 21105, 44311 Nantes Cedex 3 ) ) ) ) ) tel: (33) 02.40.37.42.29, fax: (33)

dear R users, please, consider the following code you may run: setClass("my.class1", representation( my.list = "list", my.array1 = "array"), prototype=prototype( my.list =list(), my.array1 =array(0,c(3,3,3)) )) #------------------

Many Thanks Rowland to have read and respond to this very long mail Le 18/11/2020 ? 12:46, Rowland penny via samba a ?crit?: >>> NT_STATUS_NO_LOGON_SERVERS > Are the Unix clients joined to the domain and do they use a DC as their > first nameserver ? No linux clients aren't binded by any mean to the AD Their name servers in /etc/resolv.conf are standards bind resolvers