Errors Generated by Function gof from LogisticDx Package

classic Classic list List threaded Threaded
2 messages Options
Reply | Threaded
Open this post in threaded view
|

Errors Generated by Function gof from LogisticDx Package

PaulJr
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]]

______________________________________________
[hidden email] 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.
Reply | Threaded
Open this post in threaded view
|

Re: Errors Generated by Function gof from LogisticDx Package

Bert Gunter-2
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 <[hidden email]> 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]]
>
> ______________________________________________
> [hidden email] 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]]

______________________________________________
[hidden email] 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.