similar to: Extracting values from linear models

Does anyone know R is to be cited in a publication? I've looked everywhere on the website and cannot find this. Thanks, Heather -- Heather Maughan Department of Ecology and Evolutionary Biology Biosciences West 310 University of Arizona Tucson, AZ 85701 Phone: 520-626-5108 Fax: 520-621-9190 hmaughan at u.arizona.edu

Dear R-users: I am doing multiple regressions using the "lm" function and would like to force the intercept to be equal to a specific value (such as 4.3). I was able to find out how to force it through the origin but this does not work for other values. I am also interested in forcing the regression parameters obtained from one regression in another regression with a subset of the

Hello, I am basically familiar with R and am trying to import a module that someone else has written. I know that it must go into the R library but even after I place the file there R doesn't recognize it. The module is maanova, available from the Churchill lab group for analysis of microarray data, if anyone is familiar with it. Any ideas/help? Thanks, Heather --

Dear R-users: I am using the package MAANOVA to analyze microarray data and have encountered problems when trying to plot data. I have tried emailing a MAANOVA discussion group, as well as the author of the package, and have not yet received a response so I am hoping that someone on this listserv can be of assistance. There are several functions in MAANOVA (riplot, resiplot) which call the

Dear all, Using: library(vegan) data(BCI) mod <- radfit(BCI[1,]) mod RAD models, family poisson No. of species 93, total abundance 448 par1 par2 par3 Deviance AIC BIC Null 39.5261 315.4362 315.4362 Preemption 0.042797 21.8939 299.8041 302.3367 Lognormal 1.0687 1.0186 25.1528 305.0629 310.1281

hi. why choose 1.4? correctly 1.3 FreeBSD hq.local 9.1-PRERELEASE FreeBSD 9.1-PRERELEASE #1: Tue Jul 31 19:48:36 MSK 2012 root at hq.local:/sys/amd64/compile/HQ amd64 ugen1.3: <DELL> at usbus1 uhid0: <DELL Dell UPS Rack 1920W HV, class 0/0, rev 1.10/1.00, addr 3> on usbus1 root at hq: # usbconfig ugen0.1: <EHCI root HUB Intel> at usbus0, cfg=0 md=HOST spd=HIGH (480Mbps)

Dear list, I want to perform a logistic regression analysis with multiple categorical predictors (i.e., a logit) on some data where there is a very definite relationship between one predicator and the response/independent variable. The problem I have is that in such a case the p value goes very high (while I as a naive newbie would expect it to crash towards 0). I'll illustrate my problem

Hi all, I've compiled R (version 2.8.1) from the source code with "--with- blas=" option in order to use external multi-threaded blas package. However, it's very hard to verify that compilation is correct and R is using multi-threads correctly or not. Which command or operation in R will run with multi-threads? I assumed that matrix multiplication %*% may call multi-threaded

Dear Rusers, I have met a difficult problem on explaining the differences of principal component analysis(PCA) between R,S-PLUS and SAS/STATA/SPSS, which wasn't met before. Althought they have got the same eigenvalues, their coeffiecients were different. First, I list my results from R,S-PLUS and SAS/STATA/SPSS, and then show the original dataset, hoping sb. to try and explain it.

Dear all R-users, I wanted to fit a Garch(1,1) model to a dataset by: >garch1 = garch(na.omit(dat)) But I got a warning message while executing, which is: >Warning message: >NaNs produced in: sqrt(pred$e) The garch parameters that I got are: > garch1 Call: garch(x = na.omit(dat)) Coefficient(s): a0 a1 b1 1.212e-04 1.001e+00 1.111e-14 Can any one

Recently I had the chance to look at smb service response times over OpenVPN using Wireshark/TShark. Here are 2 samples. In both samples a Win XP Pro machine is connecting to a server running samba-3.3.7 and openvpn-2.0.6_9 on FreeBSD 7.2-STABLE. Location #1 =================================================================== SMB RTT Statistics: Filter: Commands Calls Min

Dear all I'm not sure if I did the right analysis for my specific split splot design. We are studying biomass increase with different CO2 concentrations with four different functional plant groups (e.g. grasses, herbs, broad-leafed trees and conifers). Of each functional plant group we have four species. The design is orthogonal. The design is: Blocks: 2 (climate chambers, called

Dear all, I have run the confint() function according to below and I get the following error message: > confint(stepAIC.glm.spe.var.konn2.abund, level=0.95) Waiting for profiling to be done... Error: missing value where logical needed In addition: Warning message: NaNs produced in: sqrt((fm$deviance - OriginalDeviance)/DispersionParameter) or > confint(stepAIC.glm.spe.var.konn2.abund,

OS name and version, *CentOS release 6.10 (Final) Linux 2.6.32-754.3.5.el6.centos.plus.x86_64 #1 SMP Wed Aug 15 02:28:38 UTC 2018 x86_64 x86_64 x86_64 GNU/Linux* exact NUT version Name : nut Arch : x86_64 Version : 2.6.5 Release : 2.el6 Size : 4.2 M Repo : installed >From repo : epel Summary : Network UPS Tools URL :

Broken! [ 0.000000] Initializing cgroup subsys cpuset [ 0.000000] Initializing cgroup subsys cpu [ 0.000000] Initializing cgroup subsys cpuacct [ 0.000000] Linux version 3.10.0-693.11.6.el7.x86_64 (builder at kbuilder.dev.centos.org) (gcc version 4.8.5 20150623 (Red Hat 4.8.5-16) (GCC) ) #1 SMP Thu Jan 4 01:06:37 UTC 2018 [ 0.000000] Command line: root=/dev/xvda ro

Hello R users, For the last few days I am struggling with the following task: my data.frame: A1 A2 A3 B1 B2 B3 58.81 53.292 54.501 13.013 17.39 19.407 56.02 56.251 54.033 20.099 13.15 10.411 55.376 53.099 57.625 13.396 21.031 13.22 58.584 53.194 54.218 13.038 16.854 19.289 55.7 55.921 53.847 19.942 13.153 9.828 55.093 52.934

I have a two part question one about statistical theory and the other about implementations in R. Thank you for all help in advance. (1) Am I correct in understanding that Heteroscedasticity is a problem for Generalized Additive Models as it is for standard linear models? I am asking particularly about the GAMs as implemented in the mgcv package. Based upon my online search it seems that some