R/used.cars.R

#' @name used.cars
#' @aliases used.cars usa48.nb
#' @title US 1960 used car prices
#'
#' @description The \code{used.cars} data frame has 48 rows and 2 columns. The data set includes a neighbours list for the 48 states excluding DC from poly2nb().
#' 
#' @format This data frame contains the following columns:
#' \itemize{
#' \item{tax.charges} {taxes and delivery charges for 1955-9 new cars}
#' \item{price.1960} {1960 used car prices by state}
#' }
#' 
#' @source Hanna, F. A. 1966 Effects of regional differences in taxes and transport charges on automobile consumption, in Ostry, S., Rhymes, J. K. (eds) Papers on regional statistical studies, Toronto: Toronto University Press, pp. 199-223.
#' @references Hepple, L. W. 1976 A maximum likelihood model for econometric estimation with spatial series, in Masser, I (ed) Theory and practice in regional science, London: Pion, pp. 90-104.
#' @docType data
#' @keywords datasets spdep
#' 
#' @examples 
#' if (requireNamespace("spdep", quietly = TRUE)) {
#'   library(spdep)
#'   data(used.cars)
#'   moran.test(used.cars$price.1960, nb2listw(usa48.nb))
#'   moran.plot(used.cars$price.1960, nb2listw(usa48.nb),
#'            labels=rownames(used.cars))
#'   uc.lm <- lm(price.1960 ~ tax.charges, data=used.cars)
#'   summary(uc.lm)
#' 
#'   lm.morantest(uc.lm, nb2listw(usa48.nb))
#'   lm.morantest.sad(uc.lm, nb2listw(usa48.nb))
#'   lm.LMtests(uc.lm, nb2listw(usa48.nb))
#'   #if (requireNamespace("spatialreg", quietly = TRUE)) {
#'   #  library(spatialreg)
#'   #  uc.err <- errorsarlm(price.1960 ~ tax.charges, data=used.cars,
#'   #                     nb2listw(usa48.nb), tol.solve=1.0e-13, 
#'   #                     control=list(tol.opt=.Machine$double.eps^0.3))
#'   #  summary(uc.err)
#'   #  uc.lag <- lagsarlm(price.1960 ~ tax.charges, data=used.cars,
#'   #                   nb2listw(usa48.nb), tol.solve=1.0e-13, 
#'   #                   control=list(tol.opt=.Machine$double.eps^0.3))
#'   #  summary(uc.lag)
#'   #  uc.lag1 <- lagsarlm(price.1960 ~ 1, data=used.cars,
#'   #                    nb2listw(usa48.nb), tol.solve=1.0e-13, 
#'   #                    control=list(tol.opt=.Machine$double.eps^0.3))
#'   #  summary(uc.lag1)
#'   #  uc.err1 <- errorsarlm(price.1960 ~ 1, data=used.cars,
#'   #                      nb2listw(usa48.nb), tol.solve=1.0e-13, 
#'   #                      control=list(tol.opt=.Machine$double.eps^0.3))
#'   #  summary(uc.err1)
#'   #}
#' }
"used.cars"

Try the spData package in your browser

Any scripts or data that you put into this service are public.

spData documentation built on July 9, 2023, 6:24 p.m.