similar to: p-value from GLM

Full_Name: Finn Knudsen Version: 1.8.0 OS: windows 2000 Submission from: (NULL) (194.192.22.33) The problem seems to happen when running the GLM. When both multiplicative effects and an offset is present. I experienced this problem on my own dataset when using af Poisson familiy with log link function but the behaviour can be reproduced with the following code. I do not know if it is a bug, but

A querent on StackOverflow asked about the with() function http://stackoverflow.com/questions/42283479/why-when-to-use-with-function#42283479 and asked about the example in ?with library(MASS) with(anorexia, { anorex.1 <- glm(Postwt ~ Prewt + Treat + offset(Prewt), family = gaussian) summary(anorex.1) }) which saves little or no typing

Dear Listers, I got this message when compiling a package: * creating pgirmess-manual.tex ... OK * checking pgirmess-manual.text ... ERROR LaTex errors when creating DVI version. This typically indicates Rd problems. The message is quite explicit but I struggled a lot before understanding that the trouble comes from a single file "selMod.rd" among 44 topics. Even though I have

Hello, I'm a starting user of R. I have installed R 1.9.1 and winedt 5.4 If I run the example from written with winedt. The summary command does not produce any output. It does when I repeat the command manualy in R. Can someone explain me what can be the problem? library(MASS) data(anorexia) anorex.1 <- glm(Postwt ~ Prewt + Treat + offset(Prewt), family =

Hello all, I've 2 strings that representing the start and end values of a date and time. For example, time1 <- c("21/04/2005","23/05/2005","11/04/2005") time2 <- c("15/07/2009", "03/06/2008", "15/10/2005") as.difftime(time1,time2) Time differences in secs [1] NA NA NA attr(,"tzone") [1] "" How can i

Hello, I am trying to compute the mdeian of the survival time from the function survfit: > fit <- survfit(Surv(time, status) ~ 1) > fit Call: survfit(formula = Surv(time, status) ~ 1) records n.max n.start events median 0.95LCL 0.95UCL 111 111 111 20 NA NA NA The results is NA? the fit$surv gives values between 1 and 0.749! Am I doing this correct?

I would like to perform partial mantel tests on only within group values, with "between group" values assigned to NA. This is possible in package ncf partial.mantel.test, however this sues a different permutation to that used in ecodist.ecodist will not accept data with NA values, returning a "matrix is not square error. is it possible to perform this test in ecodist? many thanks

Dear all, After running a glm, I use the summary ( ) function to extract its coefficients and related statistics for further use. Unfortunately, the screen only displays a small (last) part of the results. I tried to overcome the problem by creating/saving an object "coef" for coefficients of the model and export/save it e.g. as a cvs document. While I succed with this operatiion, I do

Dear listers, I was used to package pgirmess under Windows with everything OK, but, for the first time, I had a trial this afternoon on Ubuntu 7.10 gutsy (I have a double boot computer and work more and more under unix) and R 2.7.0. Everything went OK except this: sudo R CMD check pgirmess ..... * checking examples ... ERROR Running examples in 'pgirmess-Ex.R' failed. The error most

Dear Dr.Goslee and anyone may intrested in matrix manipulate, I am using your ecodist to do mantel and partial mantel test, I have locality data and shape variation data, and the two distance matrixs are given as belowings. When I run the analysis, it is always report that the matrix is not square, but I didn't know what's wrong with my data. Would you please help me on this. I am quite

With new "strict" settings in R-devel, the following glm() example > data(anorexia, package = "MASS") > fm <- glm(Postwt ~ Prewt + Treat + offset(Prewt), family = gaussian, data = anorexia) > methods(class = class(fm)) Warning in grep(name, row.names(info)) : argument 'pattern' has length > 1 and only the first element will be used Warning in

I think this would be a good change. I think most users use the 'methods(class = <...>)' function to answer the question, "what methods can I call on objects with these classes?", and in that context I think it would be sensible for the function to accept more than one class. Kevin On Wed, Oct 17, 2018 at 7:15 AM Martin Maechler <maechler at stat.math.ethz.ch>

Hi, I looking for a method that use F or Chisq test instead of AIC in a stepwise modelo selection. I try the grasp package using the grasp.step.anova, but It dont work. > library(grasp) Carregando pacotes exigidos: gam Carregando pacotes exigidos: splines Carregando pacotes exigidos: mda Carregando pacotes exigidos: class > data(anorexia,package="MASS") > > m1 <-

R-1.7.0 built on NetBSD 1.6, but the validation test suite failed: Machinetype: Intel Pentium III (600 MHz); NetBSD 1.6 (GENERIC) Remote gcc version: gcc (GCC) 3.2.2 Remote g++ version: g++ (GCC) 3.2.2 Configure environment: CC=gcc CXX=g++ LDFLAGS=-Wl,-rpath,/usr/local/lib make[5]: Entering directory `/local/build/R-1.7.0/src/library' >>> Building/Updating

Enviado desde mi smartphone BlackBerry Z10 4G Lte.   Mensaje original   De: admin at pr.copextel.com.cu Enviado: sábado, 7 de mayo de 2016 12:30 a.m. Para: admin at pr.copextel.com.cu Asunto: Daily mail report for 2016-05-06 Grand Totals ------------ messages 409 received 4135 delivered 0 forwarded 10 deferred (114 deferrals) 14 bounced 90 rejected (2%) 0 reject warnings 0 held 0 discarded (0%)

Hi, I was trying to evaluate the pearson correlation and the p-values for an nxm matrix, where each row represents a vector. One way to do it would be to iterate through each row, and find its correlation value( and the p-value) with respect to the other rows. Is there some function by which I can use the matrix as input? Ideally, the output would be an nxn matrix, containing the p-values

Hellow all, i need help to fix this bug: ACPI BIOS Bug: Error: A valid RSDP was not found (20120711/tbxfroot-219) Because of this first errors i get this after, and it causes don''t work USB ports and some PCI cards on the system: can''t find IRQ for PCI INT A; please try using pci=biosirq Kernel: 3.6.1-xen self compiled (work perfect without xen multiboot) Xen Version: 4.2

> >Error copying file or folder - - - - - - - - - - - - - - - - - - - - - - I have not - do not wrk with linux - however, I have seen similar behaviour working with Mac (OS9 and OSX) This usually has been the result of illegal character(s) in the file or folder name, or illegal placement of legal characters. Either on the server, or on the source/destination machine. the following legal

Hi all, I am using RODB to export multiple data.frames into one excel file with multiple sheets. One thing I can't seem to unserdatand is way the command adds the character: ' To the beginning of each of my (non numeric) values. Here is an example code: library(RODBC) library(MASS) data(anorexia) save2excel <- function(x, t.name) sqlSave(xlsFile, x, tablename = t.name, rownames =

Dear all, I have installed 1.7.0 on my Win XP and want to use MASS. But I get this error message: > library(MASS) Error in namespaceExport(ns, exports) : undefined exports: abbey, accdeaths, Aids2, Animals, anorexia, austres, bacteria, beav1, beav2, biopsy, birthwt, ... Error in library(MASS) : package/namespace load failed I therefore again downloaded VR.zip. When I extracted the files I