R Dataset / Package COUNT / ships
On this Picostat.com statistics page, you will find information about the ships data set which pertains to ships. The ships data set is found in the COUNT R package. Try to load the ships data set in R by issuing the following command at the console data("ships"). This may load the data into a variable called ships. If R says the ships data set is not found, you can try installing the package by issuing this command install.packages("COUNT") and then attempt to reload the data with library("COUNT") followed by data("ships"). Perhaps strangley, if R gives you no output after entering a command, it means the command succeeded. If it succeeded you can see the data by typing ships at the command-line which should display the entire dataset.
If you need to download R, you can go to the R project website. You can download a CSV (comma separated values) version of the ships R data set. The size of this file is about 738 bytes.
Data set used in McCullagh & Nelder (1989), Hardin & Hilbe (2003), and other sources. The data contains values on the number of reported accidents for ships belonging to a company over a given time period. When a ship was constructed is also recorded.
A data frame with 40 observations on the following 7 variables.
ships is saved as a data frame. Count models use accident as the response variable, with log(service) as the offset. ship can be used as a panel identifier.
McCullagh and Nelder, 1989.
Hilbe, Joseph M (2007, 2011), Negative Binomial Regression, Cambridge University Press Hilbe, Joseph M (2009), Logistic Regression Models, Chapman & Hall/CRC Hardin, JW and JM Hilbe (2001, 2007), Generalized Linear Models and Extensions, Stata Press McCullagh, P.A, and J. Nelder (1989), Generalized Linear Models, Chapman & Hall
data(ships) glmshp <- glm(accident ~ op + co.70.74 + co.75.79 + offset(log(service)), family=poisson, data=ships) summary(glmshp) exp(coef(glmshp)) library(MASS) glmshnb <- glm.nb(accident ~ op + co.70.74 + co.75.79 + offset(log(service)), data=ships) summary(glmshnb) exp(coef(glmshnb)) ## Not run: library(gee) shipgee <- gee(accident ~ op + co.70.74 + co.75.79 + offset(log(service)), data=ships, family=poisson, corstr="exchangeable", id=ship) summary(shipgee)## End(Not run)
Dataset imported from https://www.r-project.org.
How To Register With a Username
How To Register With Google Single Sign On (SSO)
How To Login With a Username and Password
How To Login With Google Single Sign On (SSO)
How To Import a Dataset
How To Perform Statistical Analysis with Picostat
How To Use Educational Applications with Picostat
|Recent Queries For This Dataset|
No queries made on this dataset yet.
|Title||Authored on||Content type|
|OpenIntro Statistics Dataset - unemploy_pres||August 9, 2020 - 2:38 PM||Dataset|
|OpenIntro Statistics Dataset - house||August 9, 2020 - 12:25 PM||Dataset|
|R Dataset / Package evir / nidd.annual||March 9, 2018 - 1:06 PM||Dataset|
|R Dataset / Package Stat2Data / WeightLossIncentive4||March 9, 2018 - 1:06 PM||Dataset|
|R Dataset / Package HSAUR / clouds||March 9, 2018 - 1:06 PM||Dataset|