Nothing
# Copyright (C) 2006, 2008, 2010 Heather Turner
#
# 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 or 3 of the License
# (at your option).
#
# 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.
#
# A copy of the GNU General Public License is available at
# http://www.r-project.org/Licenses/
asGnm.lm <- function(object, ...) {
lmExtra <- match(c("effects", "assign", "qr", "contrasts"), names(object))
modelData <- model.frame(object)
object[lmExtra] <- NULL
object$call[[1]] <- as.name("gnm")
constrain <- which(is.na(coef(object)))
object <- c(list(formula = formula(object), eliminate = NULL,
ofInterest = NULL, na.action = na.action(modelData),
constrain = constrain,
constrainTo = numeric(length(constrain)),
family = gaussian(), predictors = fitted.values(object),
deviance = deviance(object),
y = model.response(modelData)), object)
object$terms <- gnmTerms(object$formula, data = modelData)
object$weights <- object$prior.weights <- rep.int(1, length(object$y))
object$aic <- 2 * object$rank +
object$family$aic(object$y, object$weights, object$fitted.values,
object$weights, object$deviance)
if (is.null(object$offset))
object$offset <- rep.int(0, length(coef(object)))
object$tolerance <- object$iterStart <- object$iterMax <- object$iter <-
object$converged <- "Not available - model fitted by lm()"
class(object) <- c("gnm", "glm", "lm")
object
}
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.