Paul Bernal
2019-May-06 17:35 UTC
[R] Errors Generated by Function gof from LogisticDx Package
Dear Christopher and friends, Hope you are all doing great. I am currently using R version 3.6.0 and I have a Windows 8, 64-bit Operating System. When applying function gof on my glm model, I get the following errors:> gof(GLM1)Error in factor(G, labels = dx1[, format(max(P), digits = 3), by = G]$V1) : invalid 'labels'; length 11 should be 1 or 10> gof(glm(TRANSIT ~ Draft + TOTALCOST + BUNKER + CHARTERVALUE,family=binomial, data=Dataset), plot=FALSE) Warning: glm.fit: fitted probabilities numerically 0 or 1 occurred Error in factor(G, labels = dx1[, format(max(P), digits = 3), by = G]$V1) : invalid 'labels'; length 11 should be 1 or 10 I would like to know what I am doing wrong, or if there is any issue with the gof function that I am not aware of. I provide de dput() of my dataset as well as other details below. My code is as follows:> library(LogisticDx)Registered S3 methods overwritten by 'ggplot2': method from [.quosures rlang c.quosures rlang print.quosures rlang> GLM1 <- glm(TRANSIT ~ TOTALCOST + Draft + CHARTERVALUE + BUNKER,+ family=binomial(logit), data=Dataset)> dput(Dataset)structure(list(TRANSIT = c(1L, 1L, 1L, 0L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 0L, 0L, 0L, 1L, 0L, 1L, 1L, 0L, 1L, 0L, 0L, 1L, 1L, 1L, 0L, 1L, 1L, 1L, 0L, 0L, 0L, 1L, 1L, 0L, 0L, 1L, 0L, 0L, 1L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 1L, 1L, 0L, 1L, 1L, 0L, 0L, 0L, 1L, 1L, 1L, 1L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 1L, 1L, 1L, 1L, 1L, 1L, 0L, 1L, 0L, 0L, 0L, 0L, 1L, 1L, 1L, 0L, 1L, 1L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 1L, 1L, 0L, 1L, 1L, 0L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 1L, 0L, 0L, 1L, 1L, 1L, 0L, 1L, 0L, 1L, 1L, 0L, 0L, 1L, 1L, 1L, 1L, 1L, 1L, 0L, 1L, 1L, 1L, 0L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 0L, 1L, 1L, 1L, 1L, 0L, 0L, 0L, 0L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 0L, 0L, 0L, 0L, 1L, 1L, 1L, 0L, 1L, 1L, 1L, 0L, 0L, 0L, 0L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 1L, 0L, 0L, 1L, 1L, 1L, 0L, 0L, 0L, 1L, 1L, 0L, 0L, 1L, 0L, 0L, 1L, 0L, 1L, 0L, 1L, 1L, 1L, 1L, 0L, 0L, 1L, 0L, 1L, 1L, 0L, 1L, 0L, 0L, 0L, 1L, 1L, 1L, 0L, 1L, 0L, 0L, 0L, 1L, 1L, 1L, 1L, 0L, 1L, 1L, 1L, 0L, 0L, 1L, 1L, 0L, 1L, 0L, 0L, 1L, 1L, 1L, 0L, 1L, 1L, 1L, 1L, 1L, 0L, 0L, 0L, 1L, 0L, 0L, 1L, 1L, 1L, 1L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 1L, 1L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 1L, 0L, 0L, 1L, 0L, 1L, 1L, 1L, 0L, 0L, 0L, 0L, 1L, 0L, 1L, 0L, 1L, 1L, 1L, 1L, 1L, 0L, 1L, 0L, 0L, 0L, 1L, 1L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 1L, 1L, 0L, 1L, 0L, 1L, 1L, 1L, 0L, 1L, 0L, 0L, 1L, 1L, 0L, 1L, 1L, 1L, 1L, 1L, 0L, 1L, 1L, 1L, 0L, 0L, 1L, 1L, 0L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 0L, 1L, 0L, 0L, 0L, 0L, 0L, 1L, 1L, 1L, 1L, 1L, 0L, 1L, 0L, 0L, 1L, 1L, 1L, 1L, 1L, 0L, 1L, 0L, 0L, 0L, 0L, 1L, 1L, 0L, 0L, 1L, 0L, 1L, 1L, 1L, 1L, 1L, 1L, 0L, 1L, 1L, 1L, 1L, 1L, 0L, 1L, 1L, 1L, 0L, 1L, 1L, 1L, 1L, 0L, 1L, 0L, 0L, 1L, 1L, 0L, 0L, 1L, 1L, 1L, 0L, 1L, 0L, 0L, 0L, 1L, 0L, 0L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 0L, 1L, 1L, 1L, 0L, 0L, 0L, 0L, 1L, 0L, 1L, 1L, 1L, 1L, 0L, 1L, 0L, 1L, 1L, 1L, 0L, 1L, 1L, 0L, 0L, 1L, 0L, 1L, 0L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L), Draft = c(12L, 12L, 12L, 13L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 14L, 14L, 14L, 12L, 13L, 12L, 12L, 14L, 12L, 14L, 13L, 12L, 12L, 12L, 14L, 12L, 12L, 11L, 13L, 13L, 14L, 12L, 12L, 13L, 14L, 12L, 13L, 13L, 12L, 14L, 14L, 13L, 12L, 14L, 14L, 13L, 14L, 14L, 12L, 13L, 12L, 12L, 13L, 12L, 12L, 14L, 14L, 14L, 12L, 12L, 12L, 12L, 14L, 13L, 14L, 12L, 13L, 14L, 13L, 12L, 12L, 12L, 12L, 12L, 12L, 14L, 12L, 14L, 13L, 14L, 14L, 12L, 12L, 12L, 14L, 12L, 12L, 13L, 14L, 13L, 13L, 12L, 14L, 13L, 13L, 14L, 12L, 12L, 14L, 12L, 12L, 14L, 12L, 13L, 13L, 14L, 14L, 14L, 14L, 12L, 12L, 14L, 13L, 12L, 12L, 12L, 13L, 12L, 14L, 12L, 12L, 14L, 14L, 12L, 12L, 12L, 12L, 12L, 12L, 13L, 12L, 12L, 12L, 13L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 14L, 12L, 12L, 12L, 12L, 14L, 14L, 14L, 14L, 10L, 12L, 11L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 13L, 14L, 14L, 14L, 12L, 12L, 12L, 14L, 12L, 12L, 12L, 14L, 14L, 14L, 14L, 12L, 12L, 11L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 10L, 12L, 12L, 12L, 11L, 13L, 14L, 14L, 12L, 13L, 14L, 14L, 12L, 13L, 14L, 12L, 12L, 12L, 14L, 13L, 14L, 12L, 12L, 14L, 14L, 12L, 14L, 13L, 12L, 14L, 12L, 14L, 12L, 12L, 12L, 12L, 14L, 13L, 12L, 13L, 12L, 12L, 14L, 12L, 14L, 14L, 14L, 12L, 12L, 12L, 13L, 12L, 14L, 14L, 14L, 12L, 12L, 12L, 12L, 13L, 12L, 12L, 12L, 14L, 14L, 12L, 12L, 14L, 12L, 14L, 14L, 12L, 12L, 12L, 13L, 12L, 12L, 12L, 12L, 12L, 14L, 14L, 13L, 12L, 13L, 14L, 12L, 12L, 12L, 12L, 13L, 14L, 12L, 14L, 13L, 14L, 14L, 14L, 14L, 14L, 13L, 14L, 14L, 13L, 14L, 12L, 12L, 14L, 14L, 14L, 14L, 14L, 14L, 14L, 14L, 14L, 14L, 14L, 14L, 13L, 12L, 14L, 14L, 14L, 12L, 14L, 14L, 14L, 12L, 12L, 14L, 14L, 14L, 14L, 12L, 14L, 14L, 12L, 14L, 14L, 12L, 13L, 12L, 12L, 12L, 14L, 14L, 13L, 14L, 12L, 13L, 13L, 13L, 12L, 12L, 12L, 12L, 12L, 14L, 12L, 14L, 13L, 14L, 12L, 12L, 14L, 13L, 14L, 14L, 14L, 12L, 14L, 14L, 12L, 12L, 14L, 12L, 14L, 12L, 12L, 12L, 14L, 12L, 13L, 14L, 12L, 12L, 14L, 12L, 12L, 12L, 12L, 12L, 14L, 12L, 12L, 12L, 14L, 14L, 12L, 12L, 14L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 14L, 12L, 13L, 13L, 13L, 14L, 14L, 12L, 12L, 12L, 12L, 12L, 13L, 12L, 14L, 13L, 12L, 12L, 12L, 12L, 12L, 13L, 12L, 14L, 13L, 13L, 13L, 11L, 12L, 14L, 14L, 12L, 14L, 12L, 11L, 12L, 12L, 12L, 12L, 14L, 12L, 12L, 12L, 12L, 14L, 14L, 12L, 12L, 12L, 14L, 12L, 12L, 12L, 12L, 14L, 12L, 13L, 14L, 12L, 12L, 14L, 14L, 12L, 12L, 12L, 13L, 12L, 14L, 14L, 14L, 12L, 14L, 13L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 13L, 6L, 12L, 12L, 14L, 14L, 14L, 14L, 12L, 14L, 12L, 12L, 12L, 12L, 14L, 12L, 14L, 12L, 12L, 12L, 14L, 12L, 12L, 13L, 14L, 12L, 13L, 12L, 14L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L), TOTALCOST = c(194364L, 219364L, 198260L, 237456L, 197159L, 198992L, 194337L, 219337L, 199198L, 196604L, 230607L, 196604L, 196604L, 194496L, 238600L, 236936L, 237476L, 197220L, 236950L, 197300L, 182042L, 237938L, 199221L, 237475L, 239190L, 157406L, 157211L, 182211L, 237475L, 182475L, 181599L, 156599L, 238269L, 238402L, 238069L, 161436L, 225031L, 238180L, 237572L, 189861L, 239005L, 239049L, 163814L, 240064L, 239171L, 238410L, 200878L, 239019L, 239087L, 239350L, 239352L, 240275L, 164844L, 238400L, 225158L, 202495L, 239681L, 201791L, 226863L, 244092L, 244590L, 239171L, 189811L, 219412L, 228480L, 203650L, 237514L, 247451L, 244739L, 211770L, 244308L, 239197L, 238419L, 224977L, 157362L, 162434L, 162434L, 162434L, 162434L, 239681L, 163316L, 237265L, 243920L, 244088L, 244163L, 202256L, 159592L, 201346L, 239187L, 189800L, 191959L, 239476L, 239171L, 238087L, 238052L, 164169L, 245057L, 244215L, 240812L, 239156L, 156879L, 197853L, 245367L, 164710L, 164710L, 244192L, 211110L, 239156L, 244213L, 237504L, 239018L, 241150L, 244447L, 238506L, 210298L, 243482L, 239166L, 159489L, 184600L, 226439L, 239127L, 235243L, 244296L, 159696L, 189046L, 244355L, 244446L, 187595L, 162595L, 162595L, 162595L, 170604L, 188774L, 244103L, 188680L, 163611L, 200551L, 244055L, 170606L, 169154L, 194154L, 170905L, 200551L, 191412L, 166412L, 243969L, 170483L, 210719L, 168554L, 164016L, 245158L, 245131L, 245186L, 239166L, 116360L, 155698L, 155698L, 155698L, 155698L, 223827L, 191968L, 159650L, 189999L, 201193L, 201011L, 226218L, 201218L, 243970L, 244291L, 243993L, 243993L, 236035L, 236035L, 236035L, 244070L, 159692L, 194183L, 169110L, 241994L, 238216L, 238301L, 242948L, 169810L, 189280L, 164662L, 164156L, 189156L, 163989L, 163924L, 159577L, 159650L, 170566L, 170598L, 188975L, 189006L, 99983L, 191595L, 166907L, 228744L, 166621L, 243593L, 244001L, 239035L, 172934L, 238288L, 241665L, 241665L, 193991L, 238361L, 238361L, 164215L, 168867L, 194304L, 241732L, 237745L, 237911L, 195374L, 195374L, 244044L, 244044L, 169118L, 244040L, 244040L, 198518L, 244106L, 236206L, 244136L, 191390L, 164516L, 165137L, 232682L, 244021L, 244101L, 236136L, 244101L, 194181L, 169181L, 244058L, 212313L, 238240L, 242502L, 239175L, 166221L, 184500L, 170027L, 237701L, 211035L, 244050L, 243745L, 242782L, 164482L, 166341L, 189482L, 174552L, 244213L, 190960L, 184494L, 169116L, 239123L, 239121L, 165097L, 206396L, 241738L, 165622L, 242651L, 250331L, 178778L, 169133L, 238280L, 244044L, 193182L, 194156L, 194156L, 169156L, 196240L, 244060L, 244060L, 244060L, 196050L, 243546L, 243546L, 195500L, 195500L, 170389L, 195389L, 243549L, 243503L, 211398L, 243510L, 238436L, 238546L, 243907L, 243654L, 238709L, 238656L, 244171L, 244136L, 243215L, 243957L, 243957L, 164455L, 164455L, 243287L, 238203L, 243738L, 243266L, 243294L, 243548L, 243262L, 243262L, 237628L, 243266L, 243382L, 243927L, 243574L, 168364L, 243598L, 243596L, 243647L, 191094L, 243655L, 244550L, 243907L, 200636L, 210208L, 243632L, 243632L, 243367L, 243048L, 212125L, 244651L, 243357L, 202542L, 243778L, 243502L, 170036L, 237911L, 195234L, 195220L, 170220L, 239391L, 244397L, 244397L, 238631L, 225921L, 244034L, 244051L, 243310L, 189976L, 164976L, 164976L, 164999L, 165154L, 243439L, 211003L, 244034L, 243859L, 243859L, 170008L, 175602L, 238078L, 243484L, 243619L, 243333L, 243289L, 200618L, 243392L, 243376L, 164873L, 235797L, 243930L, 191502L, 243906L, 195351L, 170527L, 195307L, 243551L, 175551L, 244759L, 238122L, 178863L, 170249L, 243701L, 200549L, 236254L, 189982L, 163055L, 203863L, 243561L, 165089L, 164574L, 193750L, 238061L, 240569L, 175435L, 164313L, 243153L, 189825L, 189825L, 164825L, 189825L, 164340L, 203691L, 168483L, 243970L, 193608L, 243054L, 243115L, 243115L, 243043L, 243115L, 201917L, 204065L, 177917L, 178745L, 178735L, 243911L, 200920L, 242726L, 243042L, 204204L, 181109L, 179157L, 200093L, 179164L, 243676L, 235476L, 243862L, 243873L, 243945L, 243927L, 168102L, 168102L, 243734L, 243929L, 179053L, 246381L, 204130L, 200546L, 200301L, 174699L, 199699L, 178309L, 243549L, 204424L, 216428L, 203785L, 204101L, 245074L, 243224L, 163661L, 179036L, 199248L, 243458L, 199190L, 200330L, 200406L, 174754L, 243138L, 195257L, 244796L, 243069L, 179132L, 204171L, 243718L, 243719L, 200616L, 175749L, 179010L, 243037L, 178405L, 243953L, 243923L, 243485L, 200891L, 239635L, 243661L, 204041L, 179002L, 204070L, 206036L, 198896L, 164487L, 166891L, 246375L, 200217L, 179153L, 210112L, 243941L, 243052L, 243724L, 246328L, 164311L, 243736L, 154373L, 192956L, 237690L, 193282L, 244901L, 198985L, 246315L, 179272L, 204007L, 202386L, 246315L, 202386L, 178856L, 243704L, 243750L, 164533L, 246330L, 204082L, 243790L, 189359L, 164359L, 168286L, 168286L, 175262L, 164395L, 189395L, 164299L, 189299L, 189110L, 154953L, 166251L, 175373L, 235883L), BUNKER = c(350L, 405L, 276L, 350L, 373L, 355L, 370L, 343L, 345L, 288L, 313L, 358L, 440L, 292L, 318L, 360L, 318L, 288L, 350L, 349L, 350L, 318L, 345L, 313L, 313L, 378L, 298L, 363L, 315L, 435L, 423L, 440L, 343L, 355L, 313L, 318L, 435L, 313L, 345L, 318L, 349L, 353L, 368L, 362L, 348L, 345L, 296L, 313L, 365L, 355L, 368L, 362L, 378L, 348L, 313L, 418L, 348L, 418L, 345L, 362L, 318L, 350L, 300L, 343L, 348L, 349L, 298L, 313L, 303L, 388L, 370L, 360L, 362L, 338L, 313L, 350L, 313L, 423L, 313L, 343L, 353L, 313L, 318L, 360L, 292L, 423L, 298L, 343L, 313L, 367L, 368L, 303L, 355L, 353L, 370L, 296L, 303L, 355L, 343L, 313L, 353L, 370L, 313L, 303L, 418L, 373L, 353L, 349L, 349L, 363L, 367L, 355L, 365L, 443L, 440L, 350L, 363L, 318L, 423L, 364L, 313L, 422L, 358L, 430L, 358L, 343L, 370L, 298L, 362L, 378L, 419L, 445L, 362L, 313L, 432L, 373L, 355L, 318L, 353L, 283L, 338L, 255L, 276L, 276L, 430L, 313L, 367L, 276L, 300L, 313L, 283L, 350L, 313L, 313L, 362L, 288L, 425L, 313L, 348L, 426L, 345L, 313L, 353L, 355L, 443L, 355L, 423L, 343L, 355L, 348L, 303L, 298L, 318L, 367L, 313L, 435L, 313L, 425L, 355L, 318L, 368L, 370L, 343L, 430L, 348L, 300L, 313L, 423L, 350L, 443L, 338L, 276L, 292L, 358L, 378L, 313L, 443L, 313L, 348L, 338L, 370L, 313L, 318L, 360L, 363L, 358L, 345L, 353L, 318L, 313L, 338L, 345L, 345L, 313L, 355L, 348L, 313L, 422L, 363L, 313L, 276L, 318L, 350L, 363L, 313L, 292L, 350L, 368L, 418L, 298L, 375L, 313L, 315L, 353L, 313L, 288L, 348L, 360L, 413L, 318L, 345L, 365L, 292L, 348L, 318L, 362L, 426L, 313L, 365L, 367L, 315L, 368L, 425L, 276L, 345L, 360L, 350L, 405L, 362L, 313L, 350L, 343L, 360L, 313L, 355L, 303L, 358L, 419L, 350L, 298L, 367L, 313L, 343L, 405L, 419L, 345L, 303L, 367L, 265L, 378L, 345L, 318L, 432L, 350L, 445L, 303L, 364L, 296L, 418L, 365L, 370L, 313L, 362L, 318L, 313L, 353L, 373L, 360L, 345L, 313L, 353L, 422L, 365L, 315L, 365L, 313L, 313L, 360L, 413L, 345L, 318L, 338L, 355L, 313L, 349L, 418L, 360L, 303L, 313L, 355L, 313L, 318L, 367L, 425L, 270L, 318L, 349L, 353L, 318L, 349L, 345L, 368L, 318L, 313L, 362L, 338L, 303L, 296L, 345L, 364L, 283L, 368L, 368L, 343L, 423L, 367L, 368L, 313L, 298L, 355L, 405L, 292L, 368L, 355L, 440L, 313L, 313L, 313L, 438L, 358L, 313L, 292L, 338L, 313L, 313L, 373L, 360L, 345L, 423L, 348L, 370L, 292L, 303L, 345L, 265L, 364L, 315L, 338L, 350L, 368L, 313L, 318L, 370L, 303L, 423L, 388L, 343L, 362L, 355L, 426L, 350L, 365L, 345L, 355L, 343L, 443L, 313L, 270L, 360L, 350L, 435L, 445L, 313L, 348L, 355L, 430L, 362L, 349L, 349L, 298L, 313L, 292L, 375L, 367L, 318L, 315L, 368L, 296L, 300L, 318L, 296L, 425L, 355L, 288L, 353L, 370L, 362L, 355L, 318L, 313L, 435L, 343L, 435L, 292L, 355L, 440L, 338L, 313L, 355L, 288L, 440L, 435L, 303L, 360L, 270L, 435L, 283L, 373L, 353L, 265L, 265L, 425L, 367L, 353L, 367L, 448L, 368L, 283L, 350L, 343L, 353L, 303L, 355L, 368L, 373L, 343L, 375L, 348L, 413L, 362L, 303L, 298L, 313L, 300L, 440L, 349L, 355L, 318L, 355L, 388L, 363L, 440L, 292L, 373L, 349L, 300L, 315L, 338L, 373L, 353L, 348L, 370L, 362L, 338L, 440L, 440L, 350L, 296L, 343L, 368L, 349L, 423L, 364L, 348L, 349L, 423L, 353L, 345L, 370L, 292L, 355L, 349L, 355L, 276L, 440L, 283L, 358L, 375L, 348L, 440L, 355L, 423L, 445L, 368L, 348L, 355L, 367L), CHARTERVALUE = c(14000L, 12825L, 10475L, 11850L, 13250L, 12100L, 11875L, 14500L, 12500L, 10500L, 13375L, 14500L, 13400L, 11000L, 12750L, 11625L, 11875L, 10500L, 11850L, 11900L, 11850L, 12750L, 12500L, 12000L, 12250L, 12750L, 10450L, 12900L, 12425L, 13375L, 12075L, 13400L, 12625L, 11125L, 12000L, 11875L, 13400L, 12000L, 12500L, 12750L, 11900L, 13625L, 12750L, 11800L, 12500L, 12500L, 9850L, 12000L, 12350L, 11125L, 12750L, 11800L, 12750L, 12500L, 12250L, 13125L, 13125L, 13125L, 12500L, 11800L, 11875L, 11850L, 11500L, 12625L, 13125L, 11900L, 10425L, 12250L, 12375L, 12400L, 11875L, 11625L, 11800L, 12400L, 12000L, 14000L, 12000L, 13125L, 12250L, 12625L, 13875L, 12400L, 11875L, 11625L, 11000L, 12075L, 10450L, 12625L, 13375L, 12875L, 13125L, 12375L, 11125L, 13625L, 11875L, 9850L, 12375L, 12100L, 14500L, 12000L, 13875L, 11875L, 12400L, 12375L, 13125L, 13250L, 13875L, 11900L, 11900L, 12900L, 12875L, 12100L, 12350L, 12375L, 13125L, 11850L, 12900L, 12750L, 13125L, 13875L, 13375L, 13025L, 14500L, 13400L, 14500L, 12625L, 11875L, 10450L, 11800L, 12750L, 12625L, 12250L, 11800L, 12250L, 13250L, 13250L, 12100L, 12750L, 13625L, 11125L, 12400L, 10250L, 10475L, 10475L, 13400L, 13375L, 12875L, 10475L, 11500L, 12400L, 11125L, 11850L, 13375L, 12400L, 11800L, 10500L, 13375L, 13375L, 12500L, 12625L, 12500L, 12000L, 13875L, 11125L, 12375L, 11125L, 13125L, 12625L, 12100L, 12500L, 12375L, 10450L, 12750L, 12875L, 12250L, 13400L, 12250L, 13375L, 11125L, 12750L, 12750L, 11875L, 14500L, 13400L, 12500L, 11500L, 13375L, 13125L, 11850L, 12375L, 12400L, 10475L, 11000L, 14500L, 12750L, 12400L, 12375L, 12250L, 12500L, 12400L, 11875L, 12250L, 12750L, 11625L, 12900L, 14500L, 12500L, 13875L, 11875L, 13375L, 12400L, 12500L, 12500L, 13375L, 12100L, 12500L, 12400L, 13025L, 12900L, 12400L, 10475L, 12750L, 11850L, 12900L, 13375L, 11000L, 11850L, 13125L, 13125L, 10450L, 12500L, 12250L, 12425L, 13875L, 12000L, 10500L, 13125L, 11625L, 12975L, 12750L, 12500L, 12350L, 11000L, 13125L, 12750L, 11800L, 12625L, 13375L, 12350L, 12875L, 12425L, 12750L, 12675L, 10475L, 12500L, 11625L, 11850L, 12825L, 11800L, 13375L, 14000L, 12625L, 11625L, 12400L, 11125L, 12375L, 14500L, 12625L, 14000L, 10425L, 12875L, 13375L, 14500L, 12825L, 12625L, 12500L, 12375L, 12875L, 9875L, 12750L, 12500L, 12750L, 13250L, 11850L, 12250L, 12375L, 13875L, 9850L, 13125L, 12350L, 11875L, 12000L, 11800L, 11875L, 12000L, 13875L, 13250L, 11625L, 12500L, 12400L, 13625L, 13025L, 12350L, 12425L, 12350L, 12400L, 12400L, 11625L, 12975L, 12500L, 11875L, 12400L, 11125L, 12000L, 11900L, 13125L, 11625L, 12375L, 12250L, 12100L, 13375L, 12750L, 12875L, 12675L, 10000L, 11875L, 11900L, 13625L, 12750L, 11900L, 12500L, 12750L, 12750L, 12000L, 11800L, 12400L, 12375L, 9850L, 12500L, 13875L, 11125L, 12750L, 12750L, 12625L, 12075L, 12875L, 13125L, 12250L, 10450L, 11125L, 12825L, 11000L, 13125L, 11125L, 13125L, 12000L, 12000L, 13375L, 13375L, 14500L, 12000L, 11000L, 12400L, 12250L, 12000L, 13250L, 11625L, 12500L, 13125L, 13125L, 11875L, 11000L, 12375L, 12500L, 9875L, 13875L, 12425L, 12400L, 14000L, 12750L, 12400L, 11875L, 11875L, 12375L, 12075L, 12400L, 14500L, 11800L, 12100L, 12625L, 14000L, 12350L, 12500L, 12100L, 12625L, 12375L, 12250L, 10000L, 11625L, 14000L, 13375L, 12250L, 13375L, 12500L, 11125L, 13400L, 11800L, 11900L, 11900L, 10425L, 12400L, 11000L, 12500L, 12875L, 12750L, 12425L, 12750L, 9850L, 11500L, 12750L, 9850L, 13375L, 11125L, 10500L, 13875L, 11875L, 11800L, 11125L, 12750L, 12250L, 13375L, 12625L, 13375L, 11000L, 11125L, 13125L, 12400L, 13375L, 12100L, 10500L, 13075L, 13375L, 12375L, 11625L, 10000L, 13400L, 11125L, 13250L, 13875L, 9875L, 9875L, 13375L, 12875L, 13875L, 12875L, 13500L, 12750L, 11125L, 11850L, 12625L, 13875L, 12375L, 12100L, 13125L, 13250L, 12625L, 12500L, 13125L, 12975L, 11800L, 12375L, 10425L, 13375L, 11500L, 13075L, 11900L, 12100L, 11875L, 11125L, 12400L, 12900L, 13400L, 11000L, 13250L, 11900L, 11500L, 12425L, 12400L, 13250L, 13625L, 12500L, 11875L, 11800L, 12400L, 13125L, 13075L, 14000L, 9850L, 14500L, 13125L, 11900L, 13125L, 13875L, 13125L, 11900L, 13125L, 13875L, 12500L, 11875L, 11000L, 11125L, 11900L, 11125L, 10475L, 13075L, 11125L, 14500L, 12500L, 13125L, 13125L, 12100L, 13125L, 12250L, 13125L, 12500L, 11125L, 12875L)), class = "data.frame", row.names = c(NA, -527L)) Best regards, Paul [[alternative HTML version deleted]]
Bert Gunter
2019-May-06 17:48 UTC
[R] Errors Generated by Function gof from LogisticDx Package
The warning message should be your hint. You cannot have gof statistics with perfect separation. See here for some suggestions: https://stats.stackexchange.com/questions/11109/how-to-deal-with-perfect-separation-in-logistic-regression Bert Gunter "The trouble with having an open mind is that people keep coming along and sticking things into it." -- Opus (aka Berkeley Breathed in his "Bloom County" comic strip ) On Mon, May 6, 2019 at 10:36 AM Paul Bernal <paulbernal07 at gmail.com> wrote:> Dear Christopher and friends, > > Hope you are all doing great. I am currently using R version 3.6.0 and I > have a Windows 8, 64-bit Operating System. > > When applying function gof on my glm model, I get the following errors: > > > gof(GLM1) > Error in factor(G, labels = dx1[, format(max(P), digits = 3), by = G]$V1) : > invalid 'labels'; length 11 should be 1 or 10 > > > gof(glm(TRANSIT ~ Draft + TOTALCOST + BUNKER + CHARTERVALUE, > family=binomial, data=Dataset), plot=FALSE) > Warning: glm.fit: fitted probabilities numerically 0 or 1 occurred > Error in factor(G, labels = dx1[, format(max(P), digits = 3), by = G]$V1) : > invalid 'labels'; length 11 should be 1 or 10 > > I would like to know what I am doing wrong, or if there is any issue with > the gof function that I am not aware of. > > I provide de dput() of my dataset as well as other details below. > > > My code is as follows: > > > library(LogisticDx) > Registered S3 methods overwritten by 'ggplot2': > method from > [.quosures rlang > c.quosures rlang > print.quosures rlang > > > GLM1 <- glm(TRANSIT ~ TOTALCOST + Draft + CHARTERVALUE + BUNKER, > + family=binomial(logit), data=Dataset) > > > dput(Dataset) > structure(list(TRANSIT = c(1L, 1L, 1L, 0L, 1L, 1L, 1L, 1L, 1L, > 1L, 1L, 1L, 1L, 1L, 0L, 0L, 0L, 1L, 0L, 1L, 1L, 0L, 1L, 0L, 0L, > 1L, 1L, 1L, 0L, 1L, 1L, 1L, 0L, 0L, 0L, 1L, 1L, 0L, 0L, 1L, 0L, > 0L, 1L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 1L, 1L, 0L, > 1L, 1L, 0L, 0L, 0L, 1L, 1L, 1L, 1L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, > 1L, 1L, 1L, 1L, 1L, 1L, 0L, 1L, 0L, 0L, 0L, 0L, 1L, 1L, 1L, 0L, > 1L, 1L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 1L, 1L, 0L, 1L, 1L, > 0L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 1L, 0L, 0L, 1L, 1L, 1L, 0L, > 1L, 0L, 1L, 1L, 0L, 0L, 1L, 1L, 1L, 1L, 1L, 1L, 0L, 1L, 1L, 1L, > 0L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 0L, 1L, 1L, 1L, 1L, 0L, 0L, 0L, > 0L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 0L, 0L, > 0L, 0L, 1L, 1L, 1L, 0L, 1L, 1L, 1L, 0L, 0L, 0L, 0L, 1L, 1L, 1L, > 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 0L, > 0L, 0L, 1L, 0L, 0L, 0L, 1L, 0L, 0L, 1L, 1L, 1L, 0L, 0L, 0L, 1L, > 1L, 0L, 0L, 1L, 0L, 0L, 1L, 0L, 1L, 0L, 1L, 1L, 1L, 1L, 0L, 0L, > 1L, 0L, 1L, 1L, 0L, 1L, 0L, 0L, 0L, 1L, 1L, 1L, 0L, 1L, 0L, 0L, > 0L, 1L, 1L, 1L, 1L, 0L, 1L, 1L, 1L, 0L, 0L, 1L, 1L, 0L, 1L, 0L, > 0L, 1L, 1L, 1L, 0L, 1L, 1L, 1L, 1L, 1L, 0L, 0L, 0L, 1L, 0L, 0L, > 1L, 1L, 1L, 1L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, > 0L, 0L, 0L, 1L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, > 0L, 0L, 1L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 1L, 1L, 0L, 0L, 0L, 0L, > 1L, 0L, 0L, 1L, 0L, 0L, 1L, 0L, 1L, 1L, 1L, 0L, 0L, 0L, 0L, 1L, > 0L, 1L, 0L, 1L, 1L, 1L, 1L, 1L, 0L, 1L, 0L, 0L, 0L, 1L, 1L, 0L, > 0L, 0L, 0L, 0L, 1L, 0L, 0L, 1L, 1L, 0L, 1L, 0L, 1L, 1L, 1L, 0L, > 1L, 0L, 0L, 1L, 1L, 0L, 1L, 1L, 1L, 1L, 1L, 0L, 1L, 1L, 1L, 0L, > 0L, 1L, 1L, 0L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 0L, 1L, 0L, 0L, 0L, > 0L, 0L, 1L, 1L, 1L, 1L, 1L, 0L, 1L, 0L, 0L, 1L, 1L, 1L, 1L, 1L, > 0L, 1L, 0L, 0L, 0L, 0L, 1L, 1L, 0L, 0L, 1L, 0L, 1L, 1L, 1L, 1L, > 1L, 1L, 0L, 1L, 1L, 1L, 1L, 1L, 0L, 1L, 1L, 1L, 0L, 1L, 1L, 1L, > 1L, 0L, 1L, 0L, 0L, 1L, 1L, 0L, 0L, 1L, 1L, 1L, 0L, 1L, 0L, 0L, > 0L, 1L, 0L, 0L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 0L, 1L, 1L, 1L, 0L, > 0L, 0L, 0L, 1L, 0L, 1L, 1L, 1L, 1L, 0L, 1L, 0L, 1L, 1L, 1L, 0L, > 1L, 1L, 0L, 0L, 1L, 0L, 1L, 0L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, > 1L, 1L, 1L, 1L, 1L, 1L), Draft = c(12L, 12L, 12L, 13L, 12L, 12L, > 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 14L, 14L, 14L, 12L, 13L, > 12L, 12L, 14L, 12L, 14L, 13L, 12L, 12L, 12L, 14L, 12L, 12L, 11L, > 13L, 13L, 14L, 12L, 12L, 13L, 14L, 12L, 13L, 13L, 12L, 14L, 14L, > 13L, 12L, 14L, 14L, 13L, 14L, 14L, 12L, 13L, 12L, 12L, 13L, 12L, > 12L, 14L, 14L, 14L, 12L, 12L, 12L, 12L, 14L, 13L, 14L, 12L, 13L, > 14L, 13L, 12L, 12L, 12L, 12L, 12L, 12L, 14L, 12L, 14L, 13L, 14L, > 14L, 12L, 12L, 12L, 14L, 12L, 12L, 13L, 14L, 13L, 13L, 12L, 14L, > 13L, 13L, 14L, 12L, 12L, 14L, 12L, 12L, 14L, 12L, 13L, 13L, 14L, > 14L, 14L, 14L, 12L, 12L, 14L, 13L, 12L, 12L, 12L, 13L, 12L, 14L, > 12L, 12L, 14L, 14L, 12L, 12L, 12L, 12L, 12L, 12L, 13L, 12L, 12L, > 12L, 13L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 14L, 12L, 12L, 12L, > 12L, 14L, 14L, 14L, 14L, 10L, 12L, 11L, 12L, 12L, 12L, 12L, 12L, > 12L, 12L, 12L, 12L, 12L, 13L, 14L, 14L, 14L, 12L, 12L, 12L, 14L, > 12L, 12L, 12L, 14L, 14L, 14L, 14L, 12L, 12L, 11L, 12L, 12L, 12L, > 12L, 12L, 12L, 12L, 12L, 12L, 12L, 10L, 12L, 12L, 12L, 11L, 13L, > 14L, 14L, 12L, 13L, 14L, 14L, 12L, 13L, 14L, 12L, 12L, 12L, 14L, > 13L, 14L, 12L, 12L, 14L, 14L, 12L, 14L, 13L, 12L, 14L, 12L, 14L, > 12L, 12L, 12L, 12L, 14L, 13L, 12L, 13L, 12L, 12L, 14L, 12L, 14L, > 14L, 14L, 12L, 12L, 12L, 13L, 12L, 14L, 14L, 14L, 12L, 12L, 12L, > 12L, 13L, 12L, 12L, 12L, 14L, 14L, 12L, 12L, 14L, 12L, 14L, 14L, > 12L, 12L, 12L, 13L, 12L, 12L, 12L, 12L, 12L, 14L, 14L, 13L, 12L, > 13L, 14L, 12L, 12L, 12L, 12L, 13L, 14L, 12L, 14L, 13L, 14L, 14L, > 14L, 14L, 14L, 13L, 14L, 14L, 13L, 14L, 12L, 12L, 14L, 14L, 14L, > 14L, 14L, 14L, 14L, 14L, 14L, 14L, 14L, 14L, 13L, 12L, 14L, 14L, > 14L, 12L, 14L, 14L, 14L, 12L, 12L, 14L, 14L, 14L, 14L, 12L, 14L, > 14L, 12L, 14L, 14L, 12L, 13L, 12L, 12L, 12L, 14L, 14L, 13L, 14L, > 12L, 13L, 13L, 13L, 12L, 12L, 12L, 12L, 12L, 14L, 12L, 14L, 13L, > 14L, 12L, 12L, 14L, 13L, 14L, 14L, 14L, 12L, 14L, 14L, 12L, 12L, > 14L, 12L, 14L, 12L, 12L, 12L, 14L, 12L, 13L, 14L, 12L, 12L, 14L, > 12L, 12L, 12L, 12L, 12L, 14L, 12L, 12L, 12L, 14L, 14L, 12L, 12L, > 14L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 14L, 12L, 13L, 13L, 13L, > 14L, 14L, 12L, 12L, 12L, 12L, 12L, 13L, 12L, 14L, 13L, 12L, 12L, > 12L, 12L, 12L, 13L, 12L, 14L, 13L, 13L, 13L, 11L, 12L, 14L, 14L, > 12L, 14L, 12L, 11L, 12L, 12L, 12L, 12L, 14L, 12L, 12L, 12L, 12L, > 14L, 14L, 12L, 12L, 12L, 14L, 12L, 12L, 12L, 12L, 14L, 12L, 13L, > 14L, 12L, 12L, 14L, 14L, 12L, 12L, 12L, 13L, 12L, 14L, 14L, 14L, > 12L, 14L, 13L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 13L, 6L, 12L, > 12L, 14L, 14L, 14L, 14L, 12L, 14L, 12L, 12L, 12L, 12L, 14L, 12L, > 14L, 12L, 12L, 12L, 14L, 12L, 12L, 13L, 14L, 12L, 13L, 12L, 14L, > 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, > 12L), TOTALCOST = c(194364L, 219364L, 198260L, 237456L, 197159L, > 198992L, 194337L, 219337L, 199198L, 196604L, 230607L, 196604L, > 196604L, 194496L, 238600L, 236936L, 237476L, 197220L, 236950L, > 197300L, 182042L, 237938L, 199221L, 237475L, 239190L, 157406L, > 157211L, 182211L, 237475L, 182475L, 181599L, 156599L, 238269L, > 238402L, 238069L, 161436L, 225031L, 238180L, 237572L, 189861L, > 239005L, 239049L, 163814L, 240064L, 239171L, 238410L, 200878L, > 239019L, 239087L, 239350L, 239352L, 240275L, 164844L, 238400L, > 225158L, 202495L, 239681L, 201791L, 226863L, 244092L, 244590L, > 239171L, 189811L, 219412L, 228480L, 203650L, 237514L, 247451L, > 244739L, 211770L, 244308L, 239197L, 238419L, 224977L, 157362L, > 162434L, 162434L, 162434L, 162434L, 239681L, 163316L, 237265L, > 243920L, 244088L, 244163L, 202256L, 159592L, 201346L, 239187L, > 189800L, 191959L, 239476L, 239171L, 238087L, 238052L, 164169L, > 245057L, 244215L, 240812L, 239156L, 156879L, 197853L, 245367L, > 164710L, 164710L, 244192L, 211110L, 239156L, 244213L, 237504L, > 239018L, 241150L, 244447L, 238506L, 210298L, 243482L, 239166L, > 159489L, 184600L, 226439L, 239127L, 235243L, 244296L, 159696L, > 189046L, 244355L, 244446L, 187595L, 162595L, 162595L, 162595L, > 170604L, 188774L, 244103L, 188680L, 163611L, 200551L, 244055L, > 170606L, 169154L, 194154L, 170905L, 200551L, 191412L, 166412L, > 243969L, 170483L, 210719L, 168554L, 164016L, 245158L, 245131L, > 245186L, 239166L, 116360L, 155698L, 155698L, 155698L, 155698L, > 223827L, 191968L, 159650L, 189999L, 201193L, 201011L, 226218L, > 201218L, 243970L, 244291L, 243993L, 243993L, 236035L, 236035L, > 236035L, 244070L, 159692L, 194183L, 169110L, 241994L, 238216L, > 238301L, 242948L, 169810L, 189280L, 164662L, 164156L, 189156L, > 163989L, 163924L, 159577L, 159650L, 170566L, 170598L, 188975L, > 189006L, 99983L, 191595L, 166907L, 228744L, 166621L, 243593L, > 244001L, 239035L, 172934L, 238288L, 241665L, 241665L, 193991L, > 238361L, 238361L, 164215L, 168867L, 194304L, 241732L, 237745L, > 237911L, 195374L, 195374L, 244044L, 244044L, 169118L, 244040L, > 244040L, 198518L, 244106L, 236206L, 244136L, 191390L, 164516L, > 165137L, 232682L, 244021L, 244101L, 236136L, 244101L, 194181L, > 169181L, 244058L, 212313L, 238240L, 242502L, 239175L, 166221L, > 184500L, 170027L, 237701L, 211035L, 244050L, 243745L, 242782L, > 164482L, 166341L, 189482L, 174552L, 244213L, 190960L, 184494L, > 169116L, 239123L, 239121L, 165097L, 206396L, 241738L, 165622L, > 242651L, 250331L, 178778L, 169133L, 238280L, 244044L, 193182L, > 194156L, 194156L, 169156L, 196240L, 244060L, 244060L, 244060L, > 196050L, 243546L, 243546L, 195500L, 195500L, 170389L, 195389L, > 243549L, 243503L, 211398L, 243510L, 238436L, 238546L, 243907L, > 243654L, 238709L, 238656L, 244171L, 244136L, 243215L, 243957L, > 243957L, 164455L, 164455L, 243287L, 238203L, 243738L, 243266L, > 243294L, 243548L, 243262L, 243262L, 237628L, 243266L, 243382L, > 243927L, 243574L, 168364L, 243598L, 243596L, 243647L, 191094L, > 243655L, 244550L, 243907L, 200636L, 210208L, 243632L, 243632L, > 243367L, 243048L, 212125L, 244651L, 243357L, 202542L, 243778L, > 243502L, 170036L, 237911L, 195234L, 195220L, 170220L, 239391L, > 244397L, 244397L, 238631L, 225921L, 244034L, 244051L, 243310L, > 189976L, 164976L, 164976L, 164999L, 165154L, 243439L, 211003L, > 244034L, 243859L, 243859L, 170008L, 175602L, 238078L, 243484L, > 243619L, 243333L, 243289L, 200618L, 243392L, 243376L, 164873L, > 235797L, 243930L, 191502L, 243906L, 195351L, 170527L, 195307L, > 243551L, 175551L, 244759L, 238122L, 178863L, 170249L, 243701L, > 200549L, 236254L, 189982L, 163055L, 203863L, 243561L, 165089L, > 164574L, 193750L, 238061L, 240569L, 175435L, 164313L, 243153L, > 189825L, 189825L, 164825L, 189825L, 164340L, 203691L, 168483L, > 243970L, 193608L, 243054L, 243115L, 243115L, 243043L, 243115L, > 201917L, 204065L, 177917L, 178745L, 178735L, 243911L, 200920L, > 242726L, 243042L, 204204L, 181109L, 179157L, 200093L, 179164L, > 243676L, 235476L, 243862L, 243873L, 243945L, 243927L, 168102L, > 168102L, 243734L, 243929L, 179053L, 246381L, 204130L, 200546L, > 200301L, 174699L, 199699L, 178309L, 243549L, 204424L, 216428L, > 203785L, 204101L, 245074L, 243224L, 163661L, 179036L, 199248L, > 243458L, 199190L, 200330L, 200406L, 174754L, 243138L, 195257L, > 244796L, 243069L, 179132L, 204171L, 243718L, 243719L, 200616L, > 175749L, 179010L, 243037L, 178405L, 243953L, 243923L, 243485L, > 200891L, 239635L, 243661L, 204041L, 179002L, 204070L, 206036L, > 198896L, 164487L, 166891L, 246375L, 200217L, 179153L, 210112L, > 243941L, 243052L, 243724L, 246328L, 164311L, 243736L, 154373L, > 192956L, 237690L, 193282L, 244901L, 198985L, 246315L, 179272L, > 204007L, 202386L, 246315L, 202386L, 178856L, 243704L, 243750L, > 164533L, 246330L, 204082L, 243790L, 189359L, 164359L, 168286L, > 168286L, 175262L, 164395L, 189395L, 164299L, 189299L, 189110L, > 154953L, 166251L, 175373L, 235883L), BUNKER = c(350L, 405L, 276L, > 350L, 373L, 355L, 370L, 343L, 345L, 288L, 313L, 358L, 440L, 292L, > 318L, 360L, 318L, 288L, 350L, 349L, 350L, 318L, 345L, 313L, 313L, > 378L, 298L, 363L, 315L, 435L, 423L, 440L, 343L, 355L, 313L, 318L, > 435L, 313L, 345L, 318L, 349L, 353L, 368L, 362L, 348L, 345L, 296L, > 313L, 365L, 355L, 368L, 362L, 378L, 348L, 313L, 418L, 348L, 418L, > 345L, 362L, 318L, 350L, 300L, 343L, 348L, 349L, 298L, 313L, 303L, > 388L, 370L, 360L, 362L, 338L, 313L, 350L, 313L, 423L, 313L, 343L, > 353L, 313L, 318L, 360L, 292L, 423L, 298L, 343L, 313L, 367L, 368L, > 303L, 355L, 353L, 370L, 296L, 303L, 355L, 343L, 313L, 353L, 370L, > 313L, 303L, 418L, 373L, 353L, 349L, 349L, 363L, 367L, 355L, 365L, > 443L, 440L, 350L, 363L, 318L, 423L, 364L, 313L, 422L, 358L, 430L, > 358L, 343L, 370L, 298L, 362L, 378L, 419L, 445L, 362L, 313L, 432L, > 373L, 355L, 318L, 353L, 283L, 338L, 255L, 276L, 276L, 430L, 313L, > 367L, 276L, 300L, 313L, 283L, 350L, 313L, 313L, 362L, 288L, 425L, > 313L, 348L, 426L, 345L, 313L, 353L, 355L, 443L, 355L, 423L, 343L, > 355L, 348L, 303L, 298L, 318L, 367L, 313L, 435L, 313L, 425L, 355L, > 318L, 368L, 370L, 343L, 430L, 348L, 300L, 313L, 423L, 350L, 443L, > 338L, 276L, 292L, 358L, 378L, 313L, 443L, 313L, 348L, 338L, 370L, > 313L, 318L, 360L, 363L, 358L, 345L, 353L, 318L, 313L, 338L, 345L, > 345L, 313L, 355L, 348L, 313L, 422L, 363L, 313L, 276L, 318L, 350L, > 363L, 313L, 292L, 350L, 368L, 418L, 298L, 375L, 313L, 315L, 353L, > 313L, 288L, 348L, 360L, 413L, 318L, 345L, 365L, 292L, 348L, 318L, > 362L, 426L, 313L, 365L, 367L, 315L, 368L, 425L, 276L, 345L, 360L, > 350L, 405L, 362L, 313L, 350L, 343L, 360L, 313L, 355L, 303L, 358L, > 419L, 350L, 298L, 367L, 313L, 343L, 405L, 419L, 345L, 303L, 367L, > 265L, 378L, 345L, 318L, 432L, 350L, 445L, 303L, 364L, 296L, 418L, > 365L, 370L, 313L, 362L, 318L, 313L, 353L, 373L, 360L, 345L, 313L, > 353L, 422L, 365L, 315L, 365L, 313L, 313L, 360L, 413L, 345L, 318L, > 338L, 355L, 313L, 349L, 418L, 360L, 303L, 313L, 355L, 313L, 318L, > 367L, 425L, 270L, 318L, 349L, 353L, 318L, 349L, 345L, 368L, 318L, > 313L, 362L, 338L, 303L, 296L, 345L, 364L, 283L, 368L, 368L, 343L, > 423L, 367L, 368L, 313L, 298L, 355L, 405L, 292L, 368L, 355L, 440L, > 313L, 313L, 313L, 438L, 358L, 313L, 292L, 338L, 313L, 313L, 373L, > 360L, 345L, 423L, 348L, 370L, 292L, 303L, 345L, 265L, 364L, 315L, > 338L, 350L, 368L, 313L, 318L, 370L, 303L, 423L, 388L, 343L, 362L, > 355L, 426L, 350L, 365L, 345L, 355L, 343L, 443L, 313L, 270L, 360L, > 350L, 435L, 445L, 313L, 348L, 355L, 430L, 362L, 349L, 349L, 298L, > 313L, 292L, 375L, 367L, 318L, 315L, 368L, 296L, 300L, 318L, 296L, > 425L, 355L, 288L, 353L, 370L, 362L, 355L, 318L, 313L, 435L, 343L, > 435L, 292L, 355L, 440L, 338L, 313L, 355L, 288L, 440L, 435L, 303L, > 360L, 270L, 435L, 283L, 373L, 353L, 265L, 265L, 425L, 367L, 353L, > 367L, 448L, 368L, 283L, 350L, 343L, 353L, 303L, 355L, 368L, 373L, > 343L, 375L, 348L, 413L, 362L, 303L, 298L, 313L, 300L, 440L, 349L, > 355L, 318L, 355L, 388L, 363L, 440L, 292L, 373L, 349L, 300L, 315L, > 338L, 373L, 353L, 348L, 370L, 362L, 338L, 440L, 440L, 350L, 296L, > 343L, 368L, 349L, 423L, 364L, 348L, 349L, 423L, 353L, 345L, 370L, > 292L, 355L, 349L, 355L, 276L, 440L, 283L, 358L, 375L, 348L, 440L, > 355L, 423L, 445L, 368L, 348L, 355L, 367L), CHARTERVALUE = c(14000L, > 12825L, 10475L, 11850L, 13250L, 12100L, 11875L, 14500L, 12500L, > 10500L, 13375L, 14500L, 13400L, 11000L, 12750L, 11625L, 11875L, > 10500L, 11850L, 11900L, 11850L, 12750L, 12500L, 12000L, 12250L, > 12750L, 10450L, 12900L, 12425L, 13375L, 12075L, 13400L, 12625L, > 11125L, 12000L, 11875L, 13400L, 12000L, 12500L, 12750L, 11900L, > 13625L, 12750L, 11800L, 12500L, 12500L, 9850L, 12000L, 12350L, > 11125L, 12750L, 11800L, 12750L, 12500L, 12250L, 13125L, 13125L, > 13125L, 12500L, 11800L, 11875L, 11850L, 11500L, 12625L, 13125L, > 11900L, 10425L, 12250L, 12375L, 12400L, 11875L, 11625L, 11800L, > 12400L, 12000L, 14000L, 12000L, 13125L, 12250L, 12625L, 13875L, > 12400L, 11875L, 11625L, 11000L, 12075L, 10450L, 12625L, 13375L, > 12875L, 13125L, 12375L, 11125L, 13625L, 11875L, 9850L, 12375L, > 12100L, 14500L, 12000L, 13875L, 11875L, 12400L, 12375L, 13125L, > 13250L, 13875L, 11900L, 11900L, 12900L, 12875L, 12100L, 12350L, > 12375L, 13125L, 11850L, 12900L, 12750L, 13125L, 13875L, 13375L, > 13025L, 14500L, 13400L, 14500L, 12625L, 11875L, 10450L, 11800L, > 12750L, 12625L, 12250L, 11800L, 12250L, 13250L, 13250L, 12100L, > 12750L, 13625L, 11125L, 12400L, 10250L, 10475L, 10475L, 13400L, > 13375L, 12875L, 10475L, 11500L, 12400L, 11125L, 11850L, 13375L, > 12400L, 11800L, 10500L, 13375L, 13375L, 12500L, 12625L, 12500L, > 12000L, 13875L, 11125L, 12375L, 11125L, 13125L, 12625L, 12100L, > 12500L, 12375L, 10450L, 12750L, 12875L, 12250L, 13400L, 12250L, > 13375L, 11125L, 12750L, 12750L, 11875L, 14500L, 13400L, 12500L, > 11500L, 13375L, 13125L, 11850L, 12375L, 12400L, 10475L, 11000L, > 14500L, 12750L, 12400L, 12375L, 12250L, 12500L, 12400L, 11875L, > 12250L, 12750L, 11625L, 12900L, 14500L, 12500L, 13875L, 11875L, > 13375L, 12400L, 12500L, 12500L, 13375L, 12100L, 12500L, 12400L, > 13025L, 12900L, 12400L, 10475L, 12750L, 11850L, 12900L, 13375L, > 11000L, 11850L, 13125L, 13125L, 10450L, 12500L, 12250L, 12425L, > 13875L, 12000L, 10500L, 13125L, 11625L, 12975L, 12750L, 12500L, > 12350L, 11000L, 13125L, 12750L, 11800L, 12625L, 13375L, 12350L, > 12875L, 12425L, 12750L, 12675L, 10475L, 12500L, 11625L, 11850L, > 12825L, 11800L, 13375L, 14000L, 12625L, 11625L, 12400L, 11125L, > 12375L, 14500L, 12625L, 14000L, 10425L, 12875L, 13375L, 14500L, > 12825L, 12625L, 12500L, 12375L, 12875L, 9875L, 12750L, 12500L, > 12750L, 13250L, 11850L, 12250L, 12375L, 13875L, 9850L, 13125L, > 12350L, 11875L, 12000L, 11800L, 11875L, 12000L, 13875L, 13250L, > 11625L, 12500L, 12400L, 13625L, 13025L, 12350L, 12425L, 12350L, > 12400L, 12400L, 11625L, 12975L, 12500L, 11875L, 12400L, 11125L, > 12000L, 11900L, 13125L, 11625L, 12375L, 12250L, 12100L, 13375L, > 12750L, 12875L, 12675L, 10000L, 11875L, 11900L, 13625L, 12750L, > 11900L, 12500L, 12750L, 12750L, 12000L, 11800L, 12400L, 12375L, > 9850L, 12500L, 13875L, 11125L, 12750L, 12750L, 12625L, 12075L, > 12875L, 13125L, 12250L, 10450L, 11125L, 12825L, 11000L, 13125L, > 11125L, 13125L, 12000L, 12000L, 13375L, 13375L, 14500L, 12000L, > 11000L, 12400L, 12250L, 12000L, 13250L, 11625L, 12500L, 13125L, > 13125L, 11875L, 11000L, 12375L, 12500L, 9875L, 13875L, 12425L, > 12400L, 14000L, 12750L, 12400L, 11875L, 11875L, 12375L, 12075L, > 12400L, 14500L, 11800L, 12100L, 12625L, 14000L, 12350L, 12500L, > 12100L, 12625L, 12375L, 12250L, 10000L, 11625L, 14000L, 13375L, > 12250L, 13375L, 12500L, 11125L, 13400L, 11800L, 11900L, 11900L, > 10425L, 12400L, 11000L, 12500L, 12875L, 12750L, 12425L, 12750L, > 9850L, 11500L, 12750L, 9850L, 13375L, 11125L, 10500L, 13875L, > 11875L, 11800L, 11125L, 12750L, 12250L, 13375L, 12625L, 13375L, > 11000L, 11125L, 13125L, 12400L, 13375L, 12100L, 10500L, 13075L, > 13375L, 12375L, 11625L, 10000L, 13400L, 11125L, 13250L, 13875L, > 9875L, 9875L, 13375L, 12875L, 13875L, 12875L, 13500L, 12750L, > 11125L, 11850L, 12625L, 13875L, 12375L, 12100L, 13125L, 13250L, > 12625L, 12500L, 13125L, 12975L, 11800L, 12375L, 10425L, 13375L, > 11500L, 13075L, 11900L, 12100L, 11875L, 11125L, 12400L, 12900L, > 13400L, 11000L, 13250L, 11900L, 11500L, 12425L, 12400L, 13250L, > 13625L, 12500L, 11875L, 11800L, 12400L, 13125L, 13075L, 14000L, > 9850L, 14500L, 13125L, 11900L, 13125L, 13875L, 13125L, 11900L, > 13125L, 13875L, 12500L, 11875L, 11000L, 11125L, 11900L, 11125L, > 10475L, 13075L, 11125L, 14500L, 12500L, 13125L, 13125L, 12100L, > 13125L, 12250L, 13125L, 12500L, 11125L, 12875L)), class = "data.frame", > row.names = c(NA, > -527L)) > > > Best regards, > > Paul > > [[alternative HTML version deleted]] > > ______________________________________________ > R-help at r-project.org mailing list -- To UNSUBSCRIBE and more, see > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide > http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. >[[alternative HTML version deleted]]