gee/0000755000176000001440000000000014531052253011042 5ustar ripleyusersgee/NAMESPACE0000644000176000001440000000030013013621652012251 0ustar ripleyusersuseDynLib(gee, .registration = TRUE) importFrom(stats, gaussian, model.extract, model.matrix, quantile) export(gee) S3method(print, gee) S3method(print, summary.gee) S3method(summary, gee) gee/README0000644000176000001440000000647313013621652011733 0ustar ripleyusersgee -- solve generalized estimating equations -- an Splus library implementation README VERSION 4.15 98/01/28 ----- 1998 Jan 26/27 -- "unification" -- condense to a single C source file, use 8.3 filenames to ease porting to windows. Change type of received as.integer data in Cgee from "int*" to "long*" so that code runs on DEC alpha. Subsequent to the "unification", Dr Brian Ripley made a number of code improvements, and these improvements are incorporated into this release (nominally 4.13 of the C source codes). The file set is now README (this file, version 4.15) COPYING (the FSF GPL for distribution control) README.too (addendum to README) blurb (for statlib indexing) ugee.c (version 4.13 of C engine) ugee.h (version 4.1 of unified header) convert.h f2c.h (headers related to f2c conversion of LINPACK) gee.d (version 4.8 of user doc) ugee.q (version 4.13 of Splus support functions) testgee.dpt tests.S testout.413 (output of tests.S) 1997 Dec 1 maintenance change A: repair offset handling per Brian Ripley 1997 (September 08) maintenance change: separate offset parameter eliminated (supply through model formula); clean up certain segments of gee() logic. Thanks to Douglas Grove for comments. 1997 (August 13) maintenance change: problem with working correlation estimation in presence of singletons resolved. Thanks to Donna Spiegelman for comments. 1998 contact: carey@biostat.harvard.edu ----- Modifications relative to last version (1994-1996) a) memory leaks diminished or eliminated b) cloglog link support added c) installation tasks standardized SHORT INSTALL/TEST DIRECTIONS FOR MOST POPULAR UNIX PLATFORMS RUNNING Splus 3.4 0)place all source code from the shar in the directory that you intend to use as a library component. you should name this directory "gee" or you may have to make some modifications to the "dyn.load" commands used in the test script. 1) Note that two versions of .First.lib are in the Splus source code file "ugee.q". By default, a version that uses "dyn.load" will be loaded. If you have a system that supports dyn.load.lib, you can comment out or remove the default .First.lib near the top of "ugee.q" and use the alternative (thanks to Brian Ripley): .First.lib <- function(...) dyn.load.lib(...) You will need to do so on SGI and DEC Alpha systems. After confirming or changing ugee.q to reflect the required .First.lib source, run the following commands Splus CHAPTER *.c *.d *.q make install make dyn.load # dyn.load.shared on SGI and DEC Alpha 2)perform the tests. The file "tests.S" is written so as to work from the source directory. If you run it from some other directory, you will have to modify the library() call. Splus < tests.S > testout.new diff testout.new testout.413 the diff should yield no differences beyond platform or S-plus version identifiers. ----- The file README.too contains miscellaneous archival information. Completion of this revision was aided by the following "beta testers": E Altpeter -- contributor of new Windows archive P Catalano D Dunson M Maechler J Maindonald -- contributor of new cloglog support P Ramsey B Ripley -- much helpful advice on portability A Rossini P Wollan The revision was stimulated by Bill Dunlap (Statsci) who created patches to rectify memory leaks. Sincere apologies to any testers who are not acknowledged here. gee/GPL-20000644000176000001440000004307613013621652011557 0ustar ripleyusers GNU GENERAL PUBLIC LICENSE Version 2, June 1991 Copyright (C) 1989, 1991 Free Software Foundation, Inc. 675 Mass Ave, Cambridge, MA 02139, USA Everyone is permitted to copy and distribute verbatim copies of this license document, but changing it is not allowed. Preamble The licenses for most software are designed to take away your freedom to share and change it. By contrast, the GNU General Public License is intended to guarantee your freedom to share and change free software--to make sure the software is free for all its users. This General Public License applies to most of the Free Software Foundation's software and to any other program whose authors commit to using it. (Some other Free Software Foundation software is covered by the GNU Library General Public License instead.) You can apply it to your programs, too. When we speak of free software, we are referring to freedom, not price. Our General Public Licenses are designed to make sure that you have the freedom to distribute copies of free software (and charge for this service if you wish), that you receive source code or can get it if you want it, that you can change the software or use pieces of it in new free programs; and that you know you can do these things. To protect your rights, we need to make restrictions that forbid anyone to deny you these rights or to ask you to surrender the rights. These restrictions translate to certain responsibilities for you if you distribute copies of the software, or if you modify it. For example, if you distribute copies of such a program, whether gratis or for a fee, you must give the recipients all the rights that you have. You must make sure that they, too, receive or can get the source code. And you must show them these terms so they know their rights. We protect your rights with two steps: (1) copyright the software, and (2) offer you this license which gives you legal permission to copy, distribute and/or modify the software. Also, for each author's protection and ours, we want to make certain that everyone understands that there is no warranty for this free software. If the software is modified by someone else and passed on, we want its recipients to know that what they have is not the original, so that any problems introduced by others will not reflect on the original authors' reputations. Finally, any free program is threatened constantly by software patents. We wish to avoid the danger that redistributors of a free program will individually obtain patent licenses, in effect making the program proprietary. To prevent this, we have made it clear that any patent must be licensed for everyone's free use or not licensed at all. The precise terms and conditions for copying, distribution and modification follow. GNU GENERAL PUBLIC LICENSE TERMS AND CONDITIONS FOR COPYING, DISTRIBUTION AND MODIFICATION 0. This License applies to any program or other work which contains a notice placed by the copyright holder saying it may be distributed under the terms of this General Public License. The "Program", below, refers to any such program or work, and a "work based on the Program" means either the Program or any derivative work under copyright law: that is to say, a work containing the Program or a portion of it, either verbatim or with modifications and/or translated into another language. (Hereinafter, translation is included without limitation in the term "modification".) Each licensee is addressed as "you". Activities other than copying, distribution and modification are not covered by this License; they are outside its scope. The act of running the Program is not restricted, and the output from the Program is covered only if its contents constitute a work based on the Program (independent of having been made by running the Program). Whether that is true depends on what the Program does. 1. You may copy and distribute verbatim copies of the Program's source code as you receive it, in any medium, provided that you conspicuously and appropriately publish on each copy an appropriate copyright notice and disclaimer of warranty; keep intact all the notices that refer to this License and to the absence of any warranty; and give any other recipients of the Program a copy of this License along with the Program. You may charge a fee for the physical act of transferring a copy, and you may at your option offer warranty protection in exchange for a fee. 2. You may modify your copy or copies of the Program or any portion of it, thus forming a work based on the Program, and copy and distribute such modifications or work under the terms of Section 1 above, provided that you also meet all of these conditions: a) You must cause the modified files to carry prominent notices stating that you changed the files and the date of any change. b) You must cause any work that you distribute or publish, that in whole or in part contains or is derived from the Program or any part thereof, to be licensed as a whole at no charge to all third parties under the terms of this License. c) If the modified program normally reads commands interactively when run, you must cause it, when started running for such interactive use in the most ordinary way, to print or display an announcement including an appropriate copyright notice and a notice that there is no warranty (or else, saying that you provide a warranty) and that users may redistribute the program under these conditions, and telling the user how to view a copy of this License. (Exception: if the Program itself is interactive but does not normally print such an announcement, your work based on the Program is not required to print an announcement.) These requirements apply to the modified work as a whole. If identifiable sections of that work are not derived from the Program, and can be reasonably considered independent and separate works in themselves, then this License, and its terms, do not apply to those sections when you distribute them as separate works. But when you distribute the same sections as part of a whole which is a work based on the Program, the distribution of the whole must be on the terms of this License, whose permissions for other licensees extend to the entire whole, and thus to each and every part regardless of who wrote it. Thus, it is not the intent of this section to claim rights or contest your rights to work written entirely by you; rather, the intent is to exercise the right to control the distribution of derivative or collective works based on the Program. In addition, mere aggregation of another work not based on the Program with the Program (or with a work based on the Program) on a volume of a storage or distribution medium does not bring the other work under the scope of this License. 3. You may copy and distribute the Program (or a work based on it, under Section 2) in object code or executable form under the terms of Sections 1 and 2 above provided that you also do one of the following: a) Accompany it with the complete corresponding machine-readable source code, which must be distributed under the terms of Sections 1 and 2 above on a medium customarily used for software interchange; or, b) Accompany it with a written offer, valid for at least three years, to give any third party, for a charge no more than your cost of physically performing source distribution, a complete machine-readable copy of the corresponding source code, to be distributed under the terms of Sections 1 and 2 above on a medium customarily used for software interchange; or, c) Accompany it with the information you received as to the offer to distribute corresponding source code. (This alternative is allowed only for noncommercial distribution and only if you received the program in object code or executable form with such an offer, in accord with Subsection b above.) The source code for a work means the preferred form of the work for making modifications to it. For an executable work, complete source code means all the source code for all modules it contains, plus any associated interface definition files, plus the scripts used to control compilation and installation of the executable. However, as a special exception, the source code distributed need not include anything that is normally distributed (in either source or binary form) with the major components (compiler, kernel, and so on) of the operating system on which the executable runs, unless that component itself accompanies the executable. If distribution of executable or object code is made by offering access to copy from a designated place, then offering equivalent access to copy the source code from the same place counts as distribution of the source code, even though third parties are not compelled to copy the source along with the object code. 4. You may not copy, modify, sublicense, or distribute the Program except as expressly provided under this License. Any attempt otherwise to copy, modify, sublicense or distribute the Program is void, and will automatically terminate your rights under this License. However, parties who have received copies, or rights, from you under this License will not have their licenses terminated so long as such parties remain in full compliance. 5. You are not required to accept this License, since you have not signed it. However, nothing else grants you permission to modify or distribute the Program or its derivative works. These actions are prohibited by law if you do not accept this License. Therefore, by modifying or distributing the Program (or any work based on the Program), you indicate your acceptance of this License to do so, and all its terms and conditions for copying, distributing or modifying the Program or works based on it. 6. Each time you redistribute the Program (or any work based on the Program), the recipient automatically receives a license from the original licensor to copy, distribute or modify the Program subject to these terms and conditions. You may not impose any further restrictions on the recipients' exercise of the rights granted herein. You are not responsible for enforcing compliance by third parties to this License. 7. If, as a consequence of a court judgment or allegation of patent infringement or for any other reason (not limited to patent issues), conditions are imposed on you (whether by court order, agreement or otherwise) that contradict the conditions of this License, they do not excuse you from the conditions of this License. If you cannot distribute so as to satisfy simultaneously your obligations under this License and any other pertinent obligations, then as a consequence you may not distribute the Program at all. For example, if a patent license would not permit royalty-free redistribution of the Program by all those who receive copies directly or indirectly through you, then the only way you could satisfy both it and this License would be to refrain entirely from distribution of the Program. If any portion of this section is held invalid or unenforceable under any particular circumstance, the balance of the section is intended to apply and the section as a whole is intended to apply in other circumstances. It is not the purpose of this section to induce you to infringe any patents or other property right claims or to contest validity of any such claims; this section has the sole purpose of protecting the integrity of the free software distribution system, which is implemented by public license practices. Many people have made generous contributions to the wide range of software distributed through that system in reliance on consistent application of that system; it is up to the author/donor to decide if he or she is willing to distribute software through any other system and a licensee cannot impose that choice. This section is intended to make thoroughly clear what is believed to be a consequence of the rest of this License. 8. If the distribution and/or use of the Program is restricted in certain countries either by patents or by copyrighted interfaces, the original copyright holder who places the Program under this License may add an explicit geographical distribution limitation excluding those countries, so that distribution is permitted only in or among countries not thus excluded. In such case, this License incorporates the limitation as if written in the body of this License. 9. The Free Software Foundation may publish revised and/or new versions of the General Public License from time to time. Such new versions will be similar in spirit to the present version, but may differ in detail to address new problems or concerns. Each version is given a distinguishing version number. If the Program specifies a version number of this License which applies to it and "any later version", you have the option of following the terms and conditions either of that version or of any later version published by the Free Software Foundation. If the Program does not specify a version number of this License, you may choose any version ever published by the Free Software Foundation. 10. If you wish to incorporate parts of the Program into other free programs whose distribution conditions are different, write to the author to ask for permission. For software which is copyrighted by the Free Software Foundation, write to the Free Software Foundation; we sometimes make exceptions for this. Our decision will be guided by the two goals of preserving the free status of all derivatives of our free software and of promoting the sharing and reuse of software generally. NO WARRANTY 11. BECAUSE THE PROGRAM IS LICENSED FREE OF CHARGE, THERE IS NO WARRANTY FOR THE PROGRAM, TO THE EXTENT PERMITTED BY APPLICABLE LAW. EXCEPT WHEN OTHERWISE STATED IN WRITING THE COPYRIGHT HOLDERS AND/OR OTHER PARTIES PROVIDE THE PROGRAM "AS IS" WITHOUT WARRANTY OF ANY KIND, EITHER EXPRESSED OR IMPLIED, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE. THE ENTIRE RISK AS TO THE QUALITY AND PERFORMANCE OF THE PROGRAM IS WITH YOU. SHOULD THE PROGRAM PROVE DEFECTIVE, YOU ASSUME THE COST OF ALL NECESSARY SERVICING, REPAIR OR CORRECTION. 12. IN NO EVENT UNLESS REQUIRED BY APPLICABLE LAW OR AGREED TO IN WRITING WILL ANY COPYRIGHT HOLDER, OR ANY OTHER PARTY WHO MAY MODIFY AND/OR REDISTRIBUTE THE PROGRAM AS PERMITTED ABOVE, BE LIABLE TO YOU FOR DAMAGES, INCLUDING ANY GENERAL, SPECIAL, INCIDENTAL OR CONSEQUENTIAL DAMAGES ARISING OUT OF THE USE OR INABILITY TO USE THE PROGRAM (INCLUDING BUT NOT LIMITED TO LOSS OF DATA OR DATA BEING RENDERED INACCURATE OR LOSSES SUSTAINED BY YOU OR THIRD PARTIES OR A FAILURE OF THE PROGRAM TO OPERATE WITH ANY OTHER PROGRAMS), EVEN IF SUCH HOLDER OR OTHER PARTY HAS BEEN ADVISED OF THE POSSIBILITY OF SUCH DAMAGES. END OF TERMS AND CONDITIONS Appendix: How to Apply These Terms to Your New Programs If you develop a new program, and you want it to be of the greatest possible use to the public, the best way to achieve this is to make it free software which everyone can redistribute and change under these terms. To do so, attach the following notices to the program. It is safest to attach them to the start of each source file to most effectively convey the exclusion of warranty; and each file should have at least the "copyright" line and a pointer to where the full notice is found. Copyright (C) 19yy This program is free software; you can redistribute it and/or modify it under the terms of the GNU General Public License as published by the Free Software Foundation; either version 2 of the License, or (at your option) any later version. This program is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License for more details. You should have received a copy of the GNU General Public License along with this program; if not, write to the Free Software Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA. Also add information on how to contact you by electronic and paper mail. If the program is interactive, make it output a short notice like this when it starts in an interactive mode: Gnomovision version 69, Copyright (C) 19yy name of author Gnomovision comes with ABSOLUTELY NO WARRANTY; for details type `show w'. This is free software, and you are welcome to redistribute it under certain conditions; type `show c' for details. The hypothetical commands `show w' and `show c' should show the appropriate parts of the General Public License. Of course, the commands you use may be called something other than `show w' and `show c'; they could even be mouse-clicks or menu items--whatever suits your program. You should also get your employer (if you work as a programmer) or your school, if any, to sign a "copyright disclaimer" for the program, if necessary. Here is a sample; alter the names: Yoyodyne, Inc., hereby disclaims all copyright interest in the program `Gnomovision' (which makes passes at compilers) written by James Hacker. , 1 April 1989 Ty Coon, President of Vice This General Public License does not permit incorporating your program into proprietary programs. If your program is a subroutine library, you may consider it more useful to permit linking proprietary applications with the library. If this is what you want to do, use the GNU Library General Public License instead of this License. gee/ChangeLog0000644000176000001440000000302713561110142012610 0ustar ripleyusersVersion 4.13-21 2019-11-07 Use Authors@R in DESCRIPTION file. Version 4.13-20 2019-11-07 Note authorship of LINPACK routines in Authors field. Version 4.13-19 2015-06-29 register symbols remove unused assignments tweak NAMESPACE imports Version 4.13-18 2012-03-14 tidy up Version 4.13-17 2011-07-11 rename COPYING to GPL-2 Version 4.13-16 2010-11-25 print.gee assumed >= 4 groups Version 4.13-15 2010-07-15 another divergence check contributed by Emmanuel Paradis Version 4.13-13 2009-07-28 add reference examples spelling, markup tweaks standardize licence. Version 4.13-13 2007-06-19 add NAMESPACE Version 4.13-12 2007-02-20 tidy up code, fix rank-deficient case Version 4.13-11 2006-09-05 use NULL not NA for R, b. use model.extract(, "response") for codetools. add back $(FLIBS) as per 'Writing R Extensions'. Version 4.13-10 2002-10-09 use (*) not (1) in Fortran use stdlib.h not malloc.h Version 4.13-9 2002-07-03 add PACKAGE= Version 4.13-8 2002-05-21 remove _, unneeded FLIBS T/F -> TRUE/FALSE in tests/tests.R Version 4.13-7 2002-04-19 replace Mathlib.h by Rmath.h remove unnecessary convert.h, zzz.R use BLAS routines, not the translations to C in ugee.c add LINPACK routines dgefa and dgedi Version 4.13-6 2001-12-08 add ... to S3 methods parent.frame not sys.parent. Version 4.13-5 2001-04-09 add maintainer field Version 4.13-4 2000-12-20 T/F -> TRUE/FALSE Re-format documentation Version 4.13-3 2000-02-07 Taken over by BDR (previously maintained by TSL) gee/man/0000755000176000001440000000000013013621652011614 5ustar ripleyusersgee/man/gee.Rd0000644000176000001440000001307613013621652012652 0ustar ripleyusers\name{gee} \alias{gee} \alias{print.gee} \alias{summary.gee} \alias{print.summary.gee} \title{ Function to solve a Generalized Estimation Equation Model } \description{ Produces an object of class \code{"gee"} which is a Generalized Estimation Equation fit of the data. } \usage{ gee(formula, id, data, subset, na.action, R = NULL, b = NULL, tol = 0.001, maxiter = 25, family = gaussian, corstr = "independence", Mv = 1, silent = TRUE, contrasts = NULL, scale.fix = FALSE, scale.value = 1, v4.4compat = FALSE) } \arguments{ \item{formula}{ a formula expression as for other regression models, of the form \code{response ~ predictors}. See the documentation of \code{\link{lm}} and \code{\link{formula}} for details. } \item{id}{ a vector which identifies the clusters. The length of \code{id} should be the same as the number of observations. Data are assumed to be sorted so that observations on a cluster are contiguous rows for all entities in the formula. } \item{data}{ an optional data frame in which to interpret the variables occurring in the \code{formula}, along with the \code{id} and \code{n} variables. } \item{subset}{ expression saying which subset of the rows of the data should be used in the fit. This can be a logical vector (which is replicated to have length equal to the number of observations), or a numeric vector indicating which observation numbers are to be included, or a character vector of the row names to be included. All observations are included by default. } \item{na.action}{ a function to filter missing data. For \code{gee} only \code{na.omit} should be used here. } \item{R}{ a square matrix of dimension maximum cluster size containing the user specified correlation. This is only appropriate if \code{corstr = "fixed"}. } \item{b}{ an initial estimate for the parameters. } \item{tol}{ the tolerance used in the fitting algorithm. } \item{maxiter}{ the maximum number of iterations. } \item{family}{ a \code{family} object: a list of functions and expressions for defining link and variance functions. Families supported in \code{gee} are \code{gaussian}, \code{binomial}, \code{poisson}, \code{Gamma}, and \code{quasi}; see the \code{\link{glm}} and \code{\link{family}} documentation. Some links are not currently available: \code{1/mu^2} and \code{sqrt} have not been hard-coded in the \samp{cgee} engine at present. The inverse gaussian variance function is not available. All combinations of remaining functions can be obtained either by family selection or by the use of \code{quasi}. } \item{corstr}{ a character string specifying the correlation structure. The following are permitted: \code{"independence"}, \code{"fixed"}, \code{"stat_M_dep"}, \code{"non_stat_M_dep"}, \code{"exchangeable"}, \code{"AR-M"} and \code{"unstructured"} } \item{Mv}{ When \code{corstr} is \code{"stat_M_dep"}, \code{"non_stat_M_dep"}, or \code{"AR-M"} then \code{Mv} must be specified. } \item{silent}{ a logical variable controlling whether parameter estimates at each iteration are printed. } \item{contrasts}{ a list giving contrasts for some or all of the factors appearing in the model formula. The elements of the list should have the same name as the variable and should be either a contrast matrix (specifically, any full-rank matrix with as many rows as there are levels in the factor), or else a function to compute such a matrix given the number of levels. } \item{scale.fix}{ a logical variable; if true, the scale parameter is fixed at the value of \code{scale.value}. } \item{scale.value}{ numeric variable giving the value to which the scale parameter should be fixed; used only if \code{scale.fix == TRUE}. } \item{v4.4compat}{ logical variable requesting compatibility of correlation parameter estimates with previous versions; the current version revises to be more faithful to the Liang and Zeger (1986) proposals (compatible with the Groemping SAS macro, version 2.03) }} \value{ An object of class \code{"gee"} representing the fit. } \section{Side Effects}{ Offsets must be specified in the model formula, as in \code{\link{glm}}. } \details{ Though input data need not be sorted by the variable named \code{"id"}, the program will interpret physically contiguous records possessing the same value of \code{id} as members of the same cluster. Thus it is possible to use the following vector as an \code{id} vector to discriminate 4 clusters of size 4: \code{c(0,0,0,0,1,1,1,1,0,0,0,0,1,1,1,1)}. } \note{ This is version 4.8 of this user documentation file, revised 98/01/27. The assistance of Dr B Ripley is gratefully acknowledged. } \references{ Liang, K.Y. and Zeger, S.L. (1986) Longitudinal data analysis using generalized linear models. \emph{Biometrika}, \bold{73} 13--22. Zeger, S.L. and Liang, K.Y. (1986) Longitudinal data analysis for discrete and continuous outcomes. \emph{Biometrics}, \bold{42} 121--130. } \seealso{ \code{\link{glm}}, \code{\link{lm}}, \code{\link{formula}}. } \examples{ data(warpbreaks) ## marginal analysis of random effects model for wool summary(gee(breaks ~ tension, id=wool, data=warpbreaks, corstr="exchangeable")) ## test for serial correlation in blocks summary(gee(breaks ~ tension, id=wool, data=warpbreaks, corstr="AR-M", Mv=1)) if(require(MASS)) { data(OME) ## not fully appropriate link for these data. (fm <- gee(cbind(Correct, Trials-Correct) ~ Loud + Age + OME, id = ID, data = OME, family = binomial, corstr = "exchangeable")) summary(fm) }} \keyword{nonlinear} gee/DESCRIPTION0000644000176000001440000000204214531052253012546 0ustar ripleyusersPackage: gee Title: Generalized Estimation Equation Solver Version: 4.13-26 Depends: stats Suggests: MASS Date: 2023-11-27 DateNote: Gee version 1998-01-27 Authors@R: c(person(c("Vincent", "J"), "Carey", role = "aut"), person(c("Thomas", "S"), "Lumley", role = "trl", comment = "R port of versions 3.13 and 4.4"), person("Cleve", "Moler", role = "ctb", comment = "LINPACK routines in src/d*"), person("Brian", "Ripley", role = c("trl", "cre", "ctb"), email = "ripley@stats.ox.ac.uk", comment = "R port of version 4.13 and updates")) Description: Generalized Estimation Equation solver. License: GPL-2 NeedsCompilation: yes Packaged: 2023-11-27 08:10:53 UTC; ripley Author: Vincent J Carey [aut], Thomas S Lumley [trl] (R port of versions 3.13 and 4.4), Cleve Moler [ctb] (LINPACK routines in src/d*), Brian Ripley [trl, cre, ctb] (R port of version 4.13 and updates) Maintainer: Brian Ripley Repository: CRAN Date/Publication: 2023-11-27 08:34:51 UTC gee/tests/0000755000176000001440000000000013013621652012203 5ustar ripleyusersgee/tests/Examples/0000755000176000001440000000000013013621652013761 5ustar ripleyusersgee/tests/Examples/gee-Ex.Rout.save0000644000176000001440000011115313013621652016705 0ustar ripleyusers R Under development (unstable) (2014-10-19 r66802) -- "Unsuffered Consequences" Copyright (C) 2014 The R Foundation for Statistical Computing Platform: x86_64-apple-darwin13.3.0 (64-bit) R is free software and comes with ABSOLUTELY NO WARRANTY. You are welcome to redistribute it under certain conditions. Type 'license()' or 'licence()' for distribution details. Natural language support but running in an English locale R is a collaborative project with many contributors. Type 'contributors()' for more information and 'citation()' on how to cite R or R packages in publications. Type 'demo()' for some demos, 'help()' for on-line help, or 'help.start()' for an HTML browser interface to help. Type 'q()' to quit R. > pkgname <- "gee" > source(file.path(R.home("share"), "R", "examples-header.R")) > options(warn = 1) > library('gee') > > base::assign(".oldSearch", base::search(), pos = 'CheckExEnv') > cleanEx() > nameEx("gee") > ### * gee > > flush(stderr()); flush(stdout()) > > ### Name: gee > ### Title: Function to solve a Generalized Estimation Equation Model > ### Aliases: gee print.gee summary.gee print.summary.gee > ### Keywords: nonlinear > > ### ** Examples > > data(warpbreaks) > ## marginal analysis of random effects model for wool > summary(gee(breaks ~ tension, id=wool, data=warpbreaks, corstr="exchangeable")) Beginning Cgee S-function, @(#) geeformula.q 4.13 98/01/27 running glm to get initial regression estimate (Intercept) tensionM tensionH 36.38889 -10.00000 -14.72222 GEE: GENERALIZED LINEAR MODELS FOR DEPENDENT DATA gee S-function, version 4.13 modified 98/01/27 (1998) Model: Link: Identity Variance to Mean Relation: Gaussian Correlation Structure: Exchangeable Call: gee(formula = breaks ~ tension, id = wool, data = warpbreaks, corstr = "exchangeable") Summary of Residuals: Min 1Q Median 3Q Max -22.388889 -8.138889 -2.666667 6.333333 33.611111 Coefficients: Estimate Naive S.E. Naive z Robust S.E. Robust z (Intercept) 36.38889 3.069434 11.855246 5.774705 6.301428 tensionM -10.00000 3.910008 -2.557539 7.463905 -1.339781 tensionH -14.72222 3.910008 -3.765266 3.731952 -3.944912 Estimated Scale Parameter: 141.1481 Number of Iterations: 1 Working Correlation [,1] [,2] [,3] [,4] [,5] [,6] [1,] 1.00000000 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 [2,] 0.02518404 1.00000000 0.02518404 0.02518404 0.02518404 0.02518404 [3,] 0.02518404 0.02518404 1.00000000 0.02518404 0.02518404 0.02518404 [4,] 0.02518404 0.02518404 0.02518404 1.00000000 0.02518404 0.02518404 [5,] 0.02518404 0.02518404 0.02518404 0.02518404 1.00000000 0.02518404 [6,] 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 1.00000000 [7,] 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 [8,] 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 [9,] 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 [10,] 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 [11,] 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 [12,] 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 [13,] 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 [14,] 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 [15,] 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 [16,] 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 [17,] 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 [18,] 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 [19,] 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 [20,] 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 [21,] 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 [22,] 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 [23,] 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 [24,] 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 [25,] 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 [26,] 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 [27,] 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 [,7] [,8] [,9] [,10] [,11] [,12] [1,] 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 [2,] 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 [3,] 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 [4,] 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 [5,] 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 [6,] 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 [7,] 1.00000000 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 [8,] 0.02518404 1.00000000 0.02518404 0.02518404 0.02518404 0.02518404 [9,] 0.02518404 0.02518404 1.00000000 0.02518404 0.02518404 0.02518404 [10,] 0.02518404 0.02518404 0.02518404 1.00000000 0.02518404 0.02518404 [11,] 0.02518404 0.02518404 0.02518404 0.02518404 1.00000000 0.02518404 [12,] 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 1.00000000 [13,] 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 [14,] 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 [15,] 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 [16,] 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 [17,] 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 [18,] 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 [19,] 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 [20,] 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 [21,] 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 [22,] 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 [23,] 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 [24,] 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 [25,] 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 [26,] 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 [27,] 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 [,13] [,14] [,15] [,16] [,17] [,18] [1,] 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 [2,] 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 [3,] 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 [4,] 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 [5,] 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 [6,] 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 [7,] 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 [8,] 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 [9,] 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 [10,] 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 [11,] 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 [12,] 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 [13,] 1.00000000 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 [14,] 0.02518404 1.00000000 0.02518404 0.02518404 0.02518404 0.02518404 [15,] 0.02518404 0.02518404 1.00000000 0.02518404 0.02518404 0.02518404 [16,] 0.02518404 0.02518404 0.02518404 1.00000000 0.02518404 0.02518404 [17,] 0.02518404 0.02518404 0.02518404 0.02518404 1.00000000 0.02518404 [18,] 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 1.00000000 [19,] 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 [20,] 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 [21,] 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 [22,] 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 [23,] 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 [24,] 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 [25,] 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 [26,] 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 [27,] 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 [,19] [,20] [,21] [,22] [,23] [,24] [1,] 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 [2,] 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 [3,] 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 [4,] 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 [5,] 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 [6,] 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 [7,] 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 [8,] 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 [9,] 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 [10,] 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 [11,] 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 [12,] 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 [13,] 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 [14,] 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 [15,] 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 [16,] 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 [17,] 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 [18,] 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 [19,] 1.00000000 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 [20,] 0.02518404 1.00000000 0.02518404 0.02518404 0.02518404 0.02518404 [21,] 0.02518404 0.02518404 1.00000000 0.02518404 0.02518404 0.02518404 [22,] 0.02518404 0.02518404 0.02518404 1.00000000 0.02518404 0.02518404 [23,] 0.02518404 0.02518404 0.02518404 0.02518404 1.00000000 0.02518404 [24,] 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 1.00000000 [25,] 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 [26,] 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 [27,] 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 0.02518404 [,25] [,26] [,27] [1,] 0.02518404 0.02518404 0.02518404 [2,] 0.02518404 0.02518404 0.02518404 [3,] 0.02518404 0.02518404 0.02518404 [4,] 0.02518404 0.02518404 0.02518404 [5,] 0.02518404 0.02518404 0.02518404 [6,] 0.02518404 0.02518404 0.02518404 [7,] 0.02518404 0.02518404 0.02518404 [8,] 0.02518404 0.02518404 0.02518404 [9,] 0.02518404 0.02518404 0.02518404 [10,] 0.02518404 0.02518404 0.02518404 [11,] 0.02518404 0.02518404 0.02518404 [12,] 0.02518404 0.02518404 0.02518404 [13,] 0.02518404 0.02518404 0.02518404 [14,] 0.02518404 0.02518404 0.02518404 [15,] 0.02518404 0.02518404 0.02518404 [16,] 0.02518404 0.02518404 0.02518404 [17,] 0.02518404 0.02518404 0.02518404 [18,] 0.02518404 0.02518404 0.02518404 [19,] 0.02518404 0.02518404 0.02518404 [20,] 0.02518404 0.02518404 0.02518404 [21,] 0.02518404 0.02518404 0.02518404 [22,] 0.02518404 0.02518404 0.02518404 [23,] 0.02518404 0.02518404 0.02518404 [24,] 0.02518404 0.02518404 0.02518404 [25,] 1.00000000 0.02518404 0.02518404 [26,] 0.02518404 1.00000000 0.02518404 [27,] 0.02518404 0.02518404 1.00000000 > ## test for serial correlation in blocks > summary(gee(breaks ~ tension, id=wool, data=warpbreaks, corstr="AR-M", Mv=1)) Beginning Cgee S-function, @(#) geeformula.q 4.13 98/01/27 running glm to get initial regression estimate (Intercept) tensionM tensionH 36.38889 -10.00000 -14.72222 GEE: GENERALIZED LINEAR MODELS FOR DEPENDENT DATA gee S-function, version 4.13 modified 98/01/27 (1998) Model: Link: Identity Variance to Mean Relation: Gaussian Correlation Structure: AR-M , M = 1 Call: gee(formula = breaks ~ tension, id = wool, data = warpbreaks, corstr = "AR-M", Mv = 1) Summary of Residuals: Min 1Q Median 3Q Max -22.413608 -8.095783 -2.688893 6.311107 33.586392 Coefficients: Estimate Naive S.E. Naive z Robust S.E. Robust z (Intercept) 36.41361 2.889306 12.602890 5.788146 6.291065 tensionM -10.09043 4.077388 -2.474730 7.513413 -1.342989 tensionH -14.72471 4.086060 -3.603646 3.758972 -3.917218 Estimated Scale Parameter: 141.1501 Number of Iterations: 2 Working Correlation [,1] [,2] [,3] [,4] [,5] [1,] 1.000000e+00 3.540253e-02 1.253339e-03 4.437136e-05 1.570858e-06 [2,] 3.540253e-02 1.000000e+00 3.540253e-02 1.253339e-03 4.437136e-05 [3,] 1.253339e-03 3.540253e-02 1.000000e+00 3.540253e-02 1.253339e-03 [4,] 4.437136e-05 1.253339e-03 3.540253e-02 1.000000e+00 3.540253e-02 [5,] 1.570858e-06 4.437136e-05 1.253339e-03 3.540253e-02 1.000000e+00 [6,] 5.561236e-08 1.570858e-06 4.437136e-05 1.253339e-03 3.540253e-02 [7,] 1.968818e-09 5.561236e-08 1.570858e-06 4.437136e-05 1.253339e-03 [8,] 6.970113e-11 1.968818e-09 5.561236e-08 1.570858e-06 4.437136e-05 [9,] 2.467596e-12 6.970113e-11 1.968818e-09 5.561236e-08 1.570858e-06 [10,] 8.735914e-14 2.467596e-12 6.970113e-11 1.968818e-09 5.561236e-08 [11,] 3.092734e-15 8.735914e-14 2.467596e-12 6.970113e-11 1.968818e-09 [12,] 1.094906e-16 3.092734e-15 8.735914e-14 2.467596e-12 6.970113e-11 [13,] 3.876244e-18 1.094906e-16 3.092734e-15 8.735914e-14 2.467596e-12 [14,] 1.372288e-19 3.876244e-18 1.094906e-16 3.092734e-15 8.735914e-14 [15,] 4.858248e-21 1.372288e-19 3.876244e-18 1.094906e-16 3.092734e-15 [16,] 1.719942e-22 4.858248e-21 1.372288e-19 3.876244e-18 1.094906e-16 [17,] 6.089031e-24 1.719942e-22 4.858248e-21 1.372288e-19 3.876244e-18 [18,] 2.155671e-25 6.089031e-24 1.719942e-22 4.858248e-21 1.372288e-19 [19,] 7.631619e-27 2.155671e-25 6.089031e-24 1.719942e-22 4.858248e-21 [20,] 2.701786e-28 7.631619e-27 2.155671e-25 6.089031e-24 1.719942e-22 [21,] 9.565005e-30 2.701786e-28 7.631619e-27 2.155671e-25 6.089031e-24 [22,] 3.386253e-31 9.565005e-30 2.701786e-28 7.631619e-27 2.155671e-25 [23,] 1.198819e-32 3.386253e-31 9.565005e-30 2.701786e-28 7.631619e-27 [24,] 4.244123e-34 1.198819e-32 3.386253e-31 9.565005e-30 2.701786e-28 [25,] 1.502527e-35 4.244123e-34 1.198819e-32 3.386253e-31 9.565005e-30 [26,] 5.319325e-37 1.502527e-35 4.244123e-34 1.198819e-32 3.386253e-31 [27,] 1.883175e-38 5.319325e-37 1.502527e-35 4.244123e-34 1.198819e-32 [,6] [,7] [,8] [,9] [,10] [1,] 5.561236e-08 1.968818e-09 6.970113e-11 2.467596e-12 8.735914e-14 [2,] 1.570858e-06 5.561236e-08 1.968818e-09 6.970113e-11 2.467596e-12 [3,] 4.437136e-05 1.570858e-06 5.561236e-08 1.968818e-09 6.970113e-11 [4,] 1.253339e-03 4.437136e-05 1.570858e-06 5.561236e-08 1.968818e-09 [5,] 3.540253e-02 1.253339e-03 4.437136e-05 1.570858e-06 5.561236e-08 [6,] 1.000000e+00 3.540253e-02 1.253339e-03 4.437136e-05 1.570858e-06 [7,] 3.540253e-02 1.000000e+00 3.540253e-02 1.253339e-03 4.437136e-05 [8,] 1.253339e-03 3.540253e-02 1.000000e+00 3.540253e-02 1.253339e-03 [9,] 4.437136e-05 1.253339e-03 3.540253e-02 1.000000e+00 3.540253e-02 [10,] 1.570858e-06 4.437136e-05 1.253339e-03 3.540253e-02 1.000000e+00 [11,] 5.561236e-08 1.570858e-06 4.437136e-05 1.253339e-03 3.540253e-02 [12,] 1.968818e-09 5.561236e-08 1.570858e-06 4.437136e-05 1.253339e-03 [13,] 6.970113e-11 1.968818e-09 5.561236e-08 1.570858e-06 4.437136e-05 [14,] 2.467596e-12 6.970113e-11 1.968818e-09 5.561236e-08 1.570858e-06 [15,] 8.735914e-14 2.467596e-12 6.970113e-11 1.968818e-09 5.561236e-08 [16,] 3.092734e-15 8.735914e-14 2.467596e-12 6.970113e-11 1.968818e-09 [17,] 1.094906e-16 3.092734e-15 8.735914e-14 2.467596e-12 6.970113e-11 [18,] 3.876244e-18 1.094906e-16 3.092734e-15 8.735914e-14 2.467596e-12 [19,] 1.372288e-19 3.876244e-18 1.094906e-16 3.092734e-15 8.735914e-14 [20,] 4.858248e-21 1.372288e-19 3.876244e-18 1.094906e-16 3.092734e-15 [21,] 1.719942e-22 4.858248e-21 1.372288e-19 3.876244e-18 1.094906e-16 [22,] 6.089031e-24 1.719942e-22 4.858248e-21 1.372288e-19 3.876244e-18 [23,] 2.155671e-25 6.089031e-24 1.719942e-22 4.858248e-21 1.372288e-19 [24,] 7.631619e-27 2.155671e-25 6.089031e-24 1.719942e-22 4.858248e-21 [25,] 2.701786e-28 7.631619e-27 2.155671e-25 6.089031e-24 1.719942e-22 [26,] 9.565005e-30 2.701786e-28 7.631619e-27 2.155671e-25 6.089031e-24 [27,] 3.386253e-31 9.565005e-30 2.701786e-28 7.631619e-27 2.155671e-25 [,11] [,12] [,13] [,14] [,15] [1,] 3.092734e-15 1.094906e-16 3.876244e-18 1.372288e-19 4.858248e-21 [2,] 8.735914e-14 3.092734e-15 1.094906e-16 3.876244e-18 1.372288e-19 [3,] 2.467596e-12 8.735914e-14 3.092734e-15 1.094906e-16 3.876244e-18 [4,] 6.970113e-11 2.467596e-12 8.735914e-14 3.092734e-15 1.094906e-16 [5,] 1.968818e-09 6.970113e-11 2.467596e-12 8.735914e-14 3.092734e-15 [6,] 5.561236e-08 1.968818e-09 6.970113e-11 2.467596e-12 8.735914e-14 [7,] 1.570858e-06 5.561236e-08 1.968818e-09 6.970113e-11 2.467596e-12 [8,] 4.437136e-05 1.570858e-06 5.561236e-08 1.968818e-09 6.970113e-11 [9,] 1.253339e-03 4.437136e-05 1.570858e-06 5.561236e-08 1.968818e-09 [10,] 3.540253e-02 1.253339e-03 4.437136e-05 1.570858e-06 5.561236e-08 [11,] 1.000000e+00 3.540253e-02 1.253339e-03 4.437136e-05 1.570858e-06 [12,] 3.540253e-02 1.000000e+00 3.540253e-02 1.253339e-03 4.437136e-05 [13,] 1.253339e-03 3.540253e-02 1.000000e+00 3.540253e-02 1.253339e-03 [14,] 4.437136e-05 1.253339e-03 3.540253e-02 1.000000e+00 3.540253e-02 [15,] 1.570858e-06 4.437136e-05 1.253339e-03 3.540253e-02 1.000000e+00 [16,] 5.561236e-08 1.570858e-06 4.437136e-05 1.253339e-03 3.540253e-02 [17,] 1.968818e-09 5.561236e-08 1.570858e-06 4.437136e-05 1.253339e-03 [18,] 6.970113e-11 1.968818e-09 5.561236e-08 1.570858e-06 4.437136e-05 [19,] 2.467596e-12 6.970113e-11 1.968818e-09 5.561236e-08 1.570858e-06 [20,] 8.735914e-14 2.467596e-12 6.970113e-11 1.968818e-09 5.561236e-08 [21,] 3.092734e-15 8.735914e-14 2.467596e-12 6.970113e-11 1.968818e-09 [22,] 1.094906e-16 3.092734e-15 8.735914e-14 2.467596e-12 6.970113e-11 [23,] 3.876244e-18 1.094906e-16 3.092734e-15 8.735914e-14 2.467596e-12 [24,] 1.372288e-19 3.876244e-18 1.094906e-16 3.092734e-15 8.735914e-14 [25,] 4.858248e-21 1.372288e-19 3.876244e-18 1.094906e-16 3.092734e-15 [26,] 1.719942e-22 4.858248e-21 1.372288e-19 3.876244e-18 1.094906e-16 [27,] 6.089031e-24 1.719942e-22 4.858248e-21 1.372288e-19 3.876244e-18 [,16] [,17] [,18] [,19] [,20] [1,] 1.719942e-22 6.089031e-24 2.155671e-25 7.631619e-27 2.701786e-28 [2,] 4.858248e-21 1.719942e-22 6.089031e-24 2.155671e-25 7.631619e-27 [3,] 1.372288e-19 4.858248e-21 1.719942e-22 6.089031e-24 2.155671e-25 [4,] 3.876244e-18 1.372288e-19 4.858248e-21 1.719942e-22 6.089031e-24 [5,] 1.094906e-16 3.876244e-18 1.372288e-19 4.858248e-21 1.719942e-22 [6,] 3.092734e-15 1.094906e-16 3.876244e-18 1.372288e-19 4.858248e-21 [7,] 8.735914e-14 3.092734e-15 1.094906e-16 3.876244e-18 1.372288e-19 [8,] 2.467596e-12 8.735914e-14 3.092734e-15 1.094906e-16 3.876244e-18 [9,] 6.970113e-11 2.467596e-12 8.735914e-14 3.092734e-15 1.094906e-16 [10,] 1.968818e-09 6.970113e-11 2.467596e-12 8.735914e-14 3.092734e-15 [11,] 5.561236e-08 1.968818e-09 6.970113e-11 2.467596e-12 8.735914e-14 [12,] 1.570858e-06 5.561236e-08 1.968818e-09 6.970113e-11 2.467596e-12 [13,] 4.437136e-05 1.570858e-06 5.561236e-08 1.968818e-09 6.970113e-11 [14,] 1.253339e-03 4.437136e-05 1.570858e-06 5.561236e-08 1.968818e-09 [15,] 3.540253e-02 1.253339e-03 4.437136e-05 1.570858e-06 5.561236e-08 [16,] 1.000000e+00 3.540253e-02 1.253339e-03 4.437136e-05 1.570858e-06 [17,] 3.540253e-02 1.000000e+00 3.540253e-02 1.253339e-03 4.437136e-05 [18,] 1.253339e-03 3.540253e-02 1.000000e+00 3.540253e-02 1.253339e-03 [19,] 4.437136e-05 1.253339e-03 3.540253e-02 1.000000e+00 3.540253e-02 [20,] 1.570858e-06 4.437136e-05 1.253339e-03 3.540253e-02 1.000000e+00 [21,] 5.561236e-08 1.570858e-06 4.437136e-05 1.253339e-03 3.540253e-02 [22,] 1.968818e-09 5.561236e-08 1.570858e-06 4.437136e-05 1.253339e-03 [23,] 6.970113e-11 1.968818e-09 5.561236e-08 1.570858e-06 4.437136e-05 [24,] 2.467596e-12 6.970113e-11 1.968818e-09 5.561236e-08 1.570858e-06 [25,] 8.735914e-14 2.467596e-12 6.970113e-11 1.968818e-09 5.561236e-08 [26,] 3.092734e-15 8.735914e-14 2.467596e-12 6.970113e-11 1.968818e-09 [27,] 1.094906e-16 3.092734e-15 8.735914e-14 2.467596e-12 6.970113e-11 [,21] [,22] [,23] [,24] [,25] [1,] 9.565005e-30 3.386253e-31 1.198819e-32 4.244123e-34 1.502527e-35 [2,] 2.701786e-28 9.565005e-30 3.386253e-31 1.198819e-32 4.244123e-34 [3,] 7.631619e-27 2.701786e-28 9.565005e-30 3.386253e-31 1.198819e-32 [4,] 2.155671e-25 7.631619e-27 2.701786e-28 9.565005e-30 3.386253e-31 [5,] 6.089031e-24 2.155671e-25 7.631619e-27 2.701786e-28 9.565005e-30 [6,] 1.719942e-22 6.089031e-24 2.155671e-25 7.631619e-27 2.701786e-28 [7,] 4.858248e-21 1.719942e-22 6.089031e-24 2.155671e-25 7.631619e-27 [8,] 1.372288e-19 4.858248e-21 1.719942e-22 6.089031e-24 2.155671e-25 [9,] 3.876244e-18 1.372288e-19 4.858248e-21 1.719942e-22 6.089031e-24 [10,] 1.094906e-16 3.876244e-18 1.372288e-19 4.858248e-21 1.719942e-22 [11,] 3.092734e-15 1.094906e-16 3.876244e-18 1.372288e-19 4.858248e-21 [12,] 8.735914e-14 3.092734e-15 1.094906e-16 3.876244e-18 1.372288e-19 [13,] 2.467596e-12 8.735914e-14 3.092734e-15 1.094906e-16 3.876244e-18 [14,] 6.970113e-11 2.467596e-12 8.735914e-14 3.092734e-15 1.094906e-16 [15,] 1.968818e-09 6.970113e-11 2.467596e-12 8.735914e-14 3.092734e-15 [16,] 5.561236e-08 1.968818e-09 6.970113e-11 2.467596e-12 8.735914e-14 [17,] 1.570858e-06 5.561236e-08 1.968818e-09 6.970113e-11 2.467596e-12 [18,] 4.437136e-05 1.570858e-06 5.561236e-08 1.968818e-09 6.970113e-11 [19,] 1.253339e-03 4.437136e-05 1.570858e-06 5.561236e-08 1.968818e-09 [20,] 3.540253e-02 1.253339e-03 4.437136e-05 1.570858e-06 5.561236e-08 [21,] 1.000000e+00 3.540253e-02 1.253339e-03 4.437136e-05 1.570858e-06 [22,] 3.540253e-02 1.000000e+00 3.540253e-02 1.253339e-03 4.437136e-05 [23,] 1.253339e-03 3.540253e-02 1.000000e+00 3.540253e-02 1.253339e-03 [24,] 4.437136e-05 1.253339e-03 3.540253e-02 1.000000e+00 3.540253e-02 [25,] 1.570858e-06 4.437136e-05 1.253339e-03 3.540253e-02 1.000000e+00 [26,] 5.561236e-08 1.570858e-06 4.437136e-05 1.253339e-03 3.540253e-02 [27,] 1.968818e-09 5.561236e-08 1.570858e-06 4.437136e-05 1.253339e-03 [,26] [,27] [1,] 5.319325e-37 1.883175e-38 [2,] 1.502527e-35 5.319325e-37 [3,] 4.244123e-34 1.502527e-35 [4,] 1.198819e-32 4.244123e-34 [5,] 3.386253e-31 1.198819e-32 [6,] 9.565005e-30 3.386253e-31 [7,] 2.701786e-28 9.565005e-30 [8,] 7.631619e-27 2.701786e-28 [9,] 2.155671e-25 7.631619e-27 [10,] 6.089031e-24 2.155671e-25 [11,] 1.719942e-22 6.089031e-24 [12,] 4.858248e-21 1.719942e-22 [13,] 1.372288e-19 4.858248e-21 [14,] 3.876244e-18 1.372288e-19 [15,] 1.094906e-16 3.876244e-18 [16,] 3.092734e-15 1.094906e-16 [17,] 8.735914e-14 3.092734e-15 [18,] 2.467596e-12 8.735914e-14 [19,] 6.970113e-11 2.467596e-12 [20,] 1.968818e-09 6.970113e-11 [21,] 5.561236e-08 1.968818e-09 [22,] 1.570858e-06 5.561236e-08 [23,] 4.437136e-05 1.570858e-06 [24,] 1.253339e-03 4.437136e-05 [25,] 3.540253e-02 1.253339e-03 [26,] 1.000000e+00 3.540253e-02 [27,] 3.540253e-02 1.000000e+00 > > if(require(MASS)) { + data(OME) + ## not fully appropriate link for these data. + (fm <- gee(cbind(Correct, Trials-Correct) ~ Loud + Age + OME, id = ID, + data = OME, family = binomial, corstr = "exchangeable")) + summary(fm) + } Loading required package: MASS Beginning Cgee S-function, @(#) geeformula.q 4.13 98/01/27 running glm to get initial regression estimate (Intercept) Loud Age OMEhigh OMElow -5.90984676 0.15516993 0.01876291 -0.07569691 -0.31660542 GEE: GENERALIZED LINEAR MODELS FOR DEPENDENT DATA gee S-function, version 4.13 modified 98/01/27 (1998) Model: Link: Logit Variance to Mean Relation: Binomial Correlation Structure: Exchangeable Call: gee(formula = cbind(Correct, Trials - Correct) ~ Loud + Age + OME, id = ID, data = OME, family = binomial, corstr = "exchangeable") Summary of Residuals: Min 1Q Median 3Q Max -0.9585811 0.5507230 2.0414189 3.1210663 12.0241951 Coefficients: Estimate Naive S.E. Naive z Robust S.E. Robust z (Intercept) -5.90100336 0.336455218 -17.5387482 0.231463368 -25.4943295 Loud 0.15507565 0.007513314 20.6401142 0.005321880 29.1392592 Age 0.01851319 0.003561640 5.1979385 0.003305525 5.6006793 OMEhigh -0.04183516 0.160288285 -0.2609995 0.152182352 -0.2749015 OMElow -0.28563252 0.131496124 -2.1721745 0.117510247 -2.4307031 Estimated Scale Parameter: 1.250466 Number of Iterations: 4 Working Correlation [,1] [,2] [,3] [,4] [,5] [,6] [1,] 1.00000000 -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [2,] -0.02344803 1.00000000 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [3,] -0.02344803 -0.02344803 1.00000000 -0.02344803 -0.02344803 -0.02344803 [4,] -0.02344803 -0.02344803 -0.02344803 1.00000000 -0.02344803 -0.02344803 [5,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 1.00000000 -0.02344803 [6,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 1.00000000 [7,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [8,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [9,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [10,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [11,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [12,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [13,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [14,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [15,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [16,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [17,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [18,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [19,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [20,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [21,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [22,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [23,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [24,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [25,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [26,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [27,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [28,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [29,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [30,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [,7] [,8] [,9] [,10] [,11] [,12] [1,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [2,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [3,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [4,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [5,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [6,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [7,] 1.00000000 -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [8,] -0.02344803 1.00000000 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [9,] -0.02344803 -0.02344803 1.00000000 -0.02344803 -0.02344803 -0.02344803 [10,] -0.02344803 -0.02344803 -0.02344803 1.00000000 -0.02344803 -0.02344803 [11,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 1.00000000 -0.02344803 [12,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 1.00000000 [13,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [14,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [15,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [16,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [17,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [18,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [19,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [20,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [21,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [22,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [23,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [24,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [25,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [26,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [27,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [28,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [29,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [30,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [,13] [,14] [,15] [,16] [,17] [,18] [1,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [2,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [3,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [4,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [5,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [6,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [7,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [8,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [9,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [10,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [11,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [12,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [13,] 1.00000000 -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [14,] -0.02344803 1.00000000 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [15,] -0.02344803 -0.02344803 1.00000000 -0.02344803 -0.02344803 -0.02344803 [16,] -0.02344803 -0.02344803 -0.02344803 1.00000000 -0.02344803 -0.02344803 [17,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 1.00000000 -0.02344803 [18,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 1.00000000 [19,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [20,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [21,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [22,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [23,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [24,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [25,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [26,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [27,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [28,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [29,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [30,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [,19] [,20] [,21] [,22] [,23] [,24] [1,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [2,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [3,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [4,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [5,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [6,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [7,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [8,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [9,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [10,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [11,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [12,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [13,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [14,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [15,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [16,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [17,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [18,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [19,] 1.00000000 -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [20,] -0.02344803 1.00000000 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [21,] -0.02344803 -0.02344803 1.00000000 -0.02344803 -0.02344803 -0.02344803 [22,] -0.02344803 -0.02344803 -0.02344803 1.00000000 -0.02344803 -0.02344803 [23,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 1.00000000 -0.02344803 [24,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 1.00000000 [25,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [26,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [27,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [28,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [29,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [30,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [,25] [,26] [,27] [,28] [,29] [,30] [1,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [2,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [3,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [4,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [5,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [6,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [7,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [8,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [9,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [10,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [11,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [12,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [13,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [14,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [15,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [16,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [17,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [18,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [19,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [20,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [21,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [22,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [23,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [24,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [25,] 1.00000000 -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [26,] -0.02344803 1.00000000 -0.02344803 -0.02344803 -0.02344803 -0.02344803 [27,] -0.02344803 -0.02344803 1.00000000 -0.02344803 -0.02344803 -0.02344803 [28,] -0.02344803 -0.02344803 -0.02344803 1.00000000 -0.02344803 -0.02344803 [29,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 1.00000000 -0.02344803 [30,] -0.02344803 -0.02344803 -0.02344803 -0.02344803 -0.02344803 1.00000000 > > > ### *