R Dataset / Package lme4 / Pastes
Attachment | Size |
---|---|
dataset-28546.csv | 1.14 KB |
Documentation |
---|
On this Picostat.com statistics page, you will find information about the Pastes data set which pertains to Paste strength by batch and cask. The Pastes data set is found in the lme4 R package. Try to load the Pastes data set in R by issuing the following command at the console data("Pastes"). This may load the data into a variable called Pastes. If R says the Pastes data set is not found, you can try installing the package by issuing this command install.packages("lme4") and then attempt to reload the data with library("lme4") followed by data("Pastes"). 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 Pastes 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 Pastes R data set. The size of this file is about 1,163 bytes. Paste strength by batch and caskDescriptionStrength of a chemical paste product; its quality depending on the delivery batch, and the cask within the delivery. FormatA data frame with 60 observations on the following 4 variables.
DetailsThe data are described in Davies and Goldsmith (1972) as coming from “ deliveries of a chemical paste product contained in casks where, in addition to sampling and testing errors, there are variations in quality between deliveries ... As a routine, three casks selected at random from each delivery were sampled and the samples were kept for reference. ... Ten of the delivery batches were sampled at random and two analytical tests carried out on each of the 30 samples”. SourceO.L. Davies and P.L. Goldsmith (eds), Statistical Methods in Research and Production, 4th ed., Oliver and Boyd, (1972), section 6.5 Examplesstr(Pastes) require(lattice) dotplot(cask ~ strength | reorder(batch, strength), Pastes, strip = FALSE, strip.left = TRUE, layout = c(1, 10), ylab = "Cask within batch", xlab = "Paste strength", jitter.y = TRUE) ## Modifying the factors to enhance the plot Pastes <- within(Pastes, batch <- reorder(batch, strength)) Pastes <- within(Pastes, sample <- reorder(reorder(sample, strength), as.numeric(batch))) dotplot(sample ~ strength | batch, Pastes, strip = FALSE, strip.left = TRUE, layout = c(1, 10), scales = list(y = list(relation = "free")), ylab = "Sample within batch", xlab = "Paste strength", jitter.y = TRUE) ## Four equivalent models differing only in specification (fm1 <- lmer(strength ~ (1|batch) + (1|sample), Pastes)) (fm2 <- lmer(strength ~ (1|batch/cask), Pastes)) (fm3 <- lmer(strength ~ (1|batch) + (1|batch:cask), Pastes)) (fm4 <- lmer(strength ~ (1|batch/sample), Pastes)) ## fm4 results in redundant labels on the sample:batch interaction head(ranef(fm4)[[1]]) ## compare to fm1 head(ranef(fm1)[[1]]) ## This model is different and NOT appropriate for these data (fm5 <- lmer(strength ~ (1|batch) + (1|cask), Pastes))L <- getME(fm1, "L") Matrix::image(L, sub = "Structure of random effects interaction in pastes model") -- Dataset imported from https://www.r-project.org. |
Picostat Manual |
---|
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 |