similar to: prediction from glm...

Hello r-help, I try to fit birds counts over years using glm. I have done (with Estate and year as factors): Model1 <- glm(Females~Estate+Year+offset = log(area)), family = quasipoisson(link = log), na.action = "na.exclude") After I have calculated the prediction using: Pred1 <- predict(Model1, type = "response", na.action = "na.exclude") My question

madame, monsieur, help me please, Je debute avec R. J'ai r??alis?? une analyse factorielle discriminante avec R. J'obtient les resultat graphiques et les tables assici??es ?? l'analyse. Cependant je souhaiterais savoir si les distances separant les groupes (decades) sont statistiquement significatives. pour cela il me faut utiliser la distance de mahalanobis et la statistique de

Un texte encapsul? et encod? dans un jeu de caract?res inconnu a ?t? nettoy?... Nom : non disponible Url : https://stat.ethz.ch/pipermail/r-help/attachments/20060609/ab46b028/attachment.pl

Hi, I have found quite a few posts on normality checking of response variables, but I am still in doubt about that. As it is easy to understand I'm not a statistician so be patient please. I want to estimate the possible effects of some predictors on my response variable that is nº of males and nº of females (cbind(males,females)), so, it would be:

could someone explain why the name of FPVAL gets " .value" concatenated onto it when the code below is run and temp is returned. I've been trying to figure this out for too long. It doesn't matter when I put the FPVAL in the return statement. It happens regardless of whether it's first or last. Thanks. f.lmmultenhanced <- function(response, pred1, pred2) {

Hi, I''m lost and my mysql connection as well... I got this error in my production.log: ActiveRecord::StatementInvalid (Mysql::Error: Lost connection to MySQL server during query: SET NAMES ''utf8''): Config: Debian Sarge (last updated), Rails 1.1.2, Webrick, mysqld Ver 4.1.11-Debian_4sarge2 for pc-linux-gnu on i386 apache2 with mod_rewrite and webrick as proxy on 8080.

hi, i'm working with mgcv packages and specially gam. My exemple is: >test<-gam(B~s(pred1)+s(pred2)) >plot(test,pages=1) when ploting test, you can view pred1 vs s(pred1, edf[1] ) & pred2 vs s(pred2, edf[2] ) I would like to know if there is a way to access to those terms (s(pred1) & s(pred2)). Does someone know how? the purpose is to access to equation of smooths terms

https://bugs.freedesktop.org/show_bug.cgi?id=105629 Bug ID: 105629 Summary: WebGL shader regular update crashes nouveau driver 17.3.3 Product: Mesa Version: 17.3 Hardware: x86-64 (AMD64) OS: Linux (All) Status: NEW Severity: normal Priority: medium Component:

Dear R Helpers I am having difficulty understanding how to use the penalty matrix for the nomROC function in package 'nonbinROC'. The documentation says that the values of the penalty matrix code the penalty function L[i,j] in which 0 <= L[i,j] <= 1 for j > i. It gives an example that if we have an ordered response with 4 categories, then we might wish to penalise larger

THANKS so very much for your help (previous and future!). I have a two follow-up questions. 1) You say that dispersion = 1 by definition ....dispersion changes from 1 to 13.5 when I go from binomial to quasibinomial....does this suggest that I should use the binomial? i.e., is the dispersion factor more important that the 2) Is there a cutoff for too much overdispersion - mine seems to be

Hi. I'm looking for samba 2.2.5 or 2.2.6 binary package to use with squid. So I need a version compile with option --with-winbind-auth-challenge and the one provide with red hat 8 doesn't have it. I tried to compile it but it doesn't work. Can someone help me? Thanks in advance, Hugues

Say I fit a logistic model and want to calculate an odds ratio between 2 sets of predictors. It is easy to obtain the difference in the predicted logodds using the predict() function, and thus get a point-estimate OR. But I can't see how to obtain the confidence interval for such an OR. For example: model <- glm(chd ~age.cat + male + lowed, family=binomial(logit)) pred1 <-

I plan on analyzing some bird point count data. I need to develop a re-sampling regime for data like this. I used scan(file="Sample.txt",what=list(0,0,0,0,0,0),sep="\t") [[1]] [1] 1 2 3 4 5 [[2]] [1] 6 7 8 9 NA [[3]] [1] 10 11 12 13 NA [[4]] [1] 14 15 16 NA NA [[5]] [1] 17 18 19 20 21 Each row represents one site and the numbers (1-21) are visits to that site.

Hi R users I am trying to get cross-tabulation tables using tables. All variables used are binary ones (0 and 1). Each time I constructed cross-tabluation table using a different variable pair (e.g., variable 1 and variable 2, variable 1 and variable 3 etc) In doing so, I ran into some problems i.V2 i.V1 1 0 17 1 33 For variable 2 (i.V2) there was no one belonging to ZERO

I'm new to R and am trying to extract the factors of a dataframe using numeric indices (e.g. df[1]) that are input to a function definition instead of the other types of references (e.g. df$out). df[1] is a list(?) whose class is "dataframe". These indexed lists can be printed successfuly but are not agreeable to the plot() and lm() functions shown below as are their df$out

I have a function that, among other things, runs a linear model and returns r2. But, the number of predictor variables passed to the function changes from 1 to 3. How can I change the formula inside the function depending on the number of variables passed in? An example: get.model.fit <- function(response.dat, pred1.dat, pred2.dat = NULL, pred3.dat = NULL) { res <- lm(response.dat ~

Hello I got this error message yesterday evening for a user I've checked LDA error log it happened 6 times in 3 days for 5 differents users, not much but I wonder what's going on as actually the server is moderately loaded ( Easter holidays ) we have currently 250/300 IMAP processes , in "normal" production time we have approx 1000 IMAP processes running. Apr 14 19:24:37

Hello, When using rsync 2.4.6 on a linux machine to copy a complete web directory (rescursively) on an other linux machine on my local network. For one file I got : skipping non-regular file www/svg_mysql/mysql.sock This file got the "s" flag set (and I suppose that is the problem). Anybody can tell me which option should I use to also in the syntax of the rsync command to copy this

I am an R beginner and trying to run a market model using event study in R framework. First, I run a market model, that is lm(stock security~SP500 index, subset=Obs[197, 396]) ->result1 Then I get predict results for a new dataset using predict (result1, newdata=Obs[397,399]) ->pred1 Pred1 should have three numbers. Now I need to calculate abnormal return by the formula stock

Hi - How can I 'manually' reproduce the results in 'pred1' below? My attempt is pred_manual, but is not correct. Any help is much appreciated. library(splines) set.seed(12345) y <- rgamma(1000, shape =0.5) age <- rnorm(1000, 45, 10) glm1 <- glm(y ~ ns(age, 4), family=Gamma(link=log)) dd <- data.frame(age = 16:80) mm <- model.matrix( ~ ns(dd$age, 4)) pred1 <-