Hello all,
I'm trying to use confint from the MASS package to compute confidence
intervals for an nls object. When I plot the results, however, they don't
make sense - lines cross over the fitted model or just don't match the
data. Code is :
Thanks for help
dat<-data.frame(a,b)
with(dat, plot(a,b))
model<-(nls(b~(1/exp(a*x))*n, data=dat, start=list(x=.001,n=20),
na.action="na.omit"))
with(dat, plot(a,b))
x1<-summary(model)$coefficients[1,1]
n1<-summary(model)$coefficients[2,1]
lines((1/exp(c(0:30)*x1)*n1))
int<-MASS:::confint.nls(model)
lines((1/exp(c(0:30)*int[1,1])*int[2,1]),lty=2)
lines((1/exp(c(0:30)*int[2,1])*int[2,2]),lty=2)
"","a"
"1",NA
"2",5.8
"3",5.15
"4",NA
"5",5.7
"6",12.35
"7",10.9
"8",10.65
"9",18.7
"10",5.35
"11",12
"12",NA
"13",14.8
"14",13.3
"15",15.6
"16",19.8
"17",NA
"18",NA
"19",NA
"20",NA
"21",NA
"22",NA
"23",NA
"24",NA
"25",NA
"26",NA
"27",NA
"28",NA
"29",7.1
"30",NA
"31",NA
"32",15.15
"33",NA
"34",11.6
"35",NA
"36",NA
"37",6.35
"38",18.45
"39",16.65
"40",NA
"41",9.8
"42",7.35
"43",7.35
"44",11.55
"45",18.1
"46",5.65
"47",4.25
"48",NA
"49",6.45
"50",5.2
"51",5.45
"52",NA
"53",NA
"54",9.8
"55",6.35
"56",21.95
"57",NA
"58",NA
"59",11.4
"60",NA
"61",NA
"62",10.65
"63",NA
"64",NA
"65",4.8
"66",7.35
"67",17.55
"68",5.55
"69",8.8
"70",8.3
"71",10.8
"72",11.1
"73",5.3
"74",NA
"75",NA
"76",7.55
"77",8.5
"78",NA
"79",21.55
"80",NA
"81",7
"82",5.4
"83",NA
"84",NA
"85",7.45
"86",5.5
"87",8.05
"88",16.35
"89",5.25
"90",6.85
"91",4.7
"92",3.25
"93",5.45
"94",8.15
"95",2.9
"96",3.05
"97",7.9
"98",6.2
"99",5.85
"100",6.7
"101",17.45
"102",4.3
"103",9.85
"104",12.6
"105",12.75
"106",3.95
"107",NA
"108",14.4
"109",7.4
"110",NA
"111",4.95
"112",4.75
"113",4.2
"114",5.5
"115",6.2
"116",9.6
"117",NA
"118",12.95
"119",9.1
"120",25.9
"121",7.55
"122",12.3
"123",11.2
"124",8.35
"125",NA
"126",6.95
"127",20.9
"128",NA
"129",NA
"130",6.75
"131",6.55
"132",NA
"133",10.25
"134",8.3
"135",NA
"136",10.35
"137",4.8
"138",12.25
"139",9.45
"140",7.55
"141",5.15
"142",9
"143",4.95
"144",NA
"145",8.4
"146",6.65
"147",11.15
"148",7.65
"149",NA
"150",NA
"151",NA
"152",NA
"153",12.95
"154",NA
"155",4.6
"156",4.8
"157",10.95
"158",4.25
"159",4.1
"160",5.3
"161",3.8
"162",12.75
"163",NA
"164",10.8
"165",NA
"166",8.65
"167",25.45
"168",NA
"169",7.65
"170",18.4
"171",9.1
"172",5.7
"173",6.15
"174",NA
"175",NA
"176",4.5
"177",8.45
"178",NA
"179",NA
"180",NA
"181",8.55
"182",9.15
"183",12.9
"184",27.65
"185",NA
"186",8.3
"187",12.8
"188",NA
"189",5.85
"190",NA
"191",7
"192",11.65
"193",8.15
"194",24.45
"195",15.25
"196",10.5
"197",NA
"198",3.6
"199",NA
"200",5.55
"201",NA
"202",10.3
"203",NA
"204",4.6
"205",5.95
"206",10.15
"207",5.7
"208",23.2
"209",11.95
"210",NA
"211",6.05
"212",5.55
"213",NA
"214",NA
"215",NA
"216",NA
"217",16.4
"218",4.65
"219",16.35
"220",7.95
"221",7.45
"222",NA
"223",7.9
"224",21.15
"225",NA
"226",NA
"227",NA
"228",11.15
"229",9.15
"230",11.35
"231",4.75
"232",3.85
"233",NA
"234",5.05
"235",NA
"236",12.7
"237",6.35
"238",NA
"239",NA
"240",NA
"241",10.15
"242",5.15
"243",NA
"244",NA
"245",6.4
"246",7.15
"247",NA
"248",16.5
"249",NA
"250",NA
"251",NA
"252",7.15
"253",6.7
"254",NA
"255",NA
"256",NA
"257",NA
"258",NA
"259",NA
"260",NA
"261",4.1
"262",4.55
"263",4.45
"264",NA
"265",7.2
"266",6.1
"267",6
"268",NA
"269",8.75
"270",NA
"271",5.6
"272",NA
"273",NA
"274",NA
"275",NA
"276",NA
"277",NA
"278",NA
"279",NA
"280",NA
"281",NA
"282",NA
"283",NA
"284",NA
"285",NA
"286",NA
"287",NA
"288",8.1
"289",NA
"290",NA
"291",NA
"292",NA
"293",NA
"294",NA
"295",NA
"296",NA
"297",8.15
"298",7.7
"299",NA
"300",NA
"301",26.75
"302",21
"303",10.45
"304",NA
"305",6.7
"306",20.5
"307",NA
"308",10.5
"309",NA
"310",NA
"311",NA
"312",NA
"313",NA
"314",5.65
"315",12.75
"316",5
"317",5.85
"318",3.55
"319",4.05
"320",6.4
"321",4.75
"322",NA
"323",NA
"324",NA
"325",2.75
"326",NA
"327",4.9
"328",5.05
"329",5.3
"330",NA
"331",NA
"332",NA
"333",NA
"334",NA
"335",6.95
"336",6.4
"337",6.1
"338",NA
"339",11.75
"340",3.35
"341",7.45
"342",6.9
"343",5.75
"344",6.1
"345",NA
"346",NA
"347",5
"348",5.85
"349",NA
"350",5.5
"351",6.65
"352",11.1
"353",10.95
"354",16.6
"355",7.2
"356",6.5
"357",NA
"358",8.45
"359",NA
"360",8.3
"361",17.05
"362",NA
"363",NA
"364",NA
"365",NA
"366",12.7
"367",NA
"368",8.05
"369",9.7
"370",NA
"371",12.35
"372",4.85
"373",5.5
"374",NA
"375",4.9
"376",NA
"377",NA
"378",15.05
"379",NA
"380",NA
"381",8.4
"382",NA
"383",NA
"384",5.15
"385",11.4
"386",NA
"387",8
"388",7.45
"389",8.7
"390",9.05
"391",15.1
"392",NA
"393",NA
"394",NA
"395",NA
"396",7.4
"397",10.6
"398",20.2
"399",10.95
"400",7.6
"401",NA
"402",NA
"403",NA
"404",8.3
"405",9.65
"406",14.75
"407",NA
"408",2.9
"409",NA
"410",NA
"411",NA
"412",NA
"413",NA
"414",NA
"415",NA
"416",2.55
"417",3.2
"418",4.05
"419",9.05
"420",3.25
"421",6.7
"422",NA
"423",NA
"424",NA
"425",7.25
"426",NA
"427",NA
"428",NA
"429",8.8
"430",5.05
"431",3.6
"432",7.7
"433",14.45
"434",10.2
"435",NA
"436",NA
"437",12.95
"438",8.8
"439",4.5
"440",5.3
"441",17.7
"442",10.2
"443",5.55
"444",3.9
"445",20
"446",12.2
"447",10.8
"448",NA
"449",NA
"450",NA
"451",NA
"452",NA
"453",NA
"454",NA
"455",NA
"456",NA
"457",NA
"458",NA
"459",NA
"460",NA
"461",NA
"462",NA
"463",NA
"464",NA
"","b"
"1",NA
"2",243.085758119756
"3",274.339641306582
"4",NA
"5",418.039560741661
"6",46.5033624229099
"7",79.9087810465099
"8",215.216454070312
"9",107.035742459874
"10",278.567219708664
"11",283.033857603722
"12",NA
"13",98.2911978484231
"14",173.670430606987
"15",174.110803357204
"16",95.6489613471215
"17",NA
"18",NA
"19",NA
"20",NA
"21",NA
"22",NA
"23",NA
"24",NA
"25",NA
"26",NA
"27",NA
"28",NA
"29",222.264934489499
"30",NA
"31",NA
"32",181.181931517831
"33",NA
"34",189.486103379064
"35",NA
"36",NA
"37",389.654391470535
"38",41.7095904848339
"39",199.224632198148
"40",NA
"41",144.781978192755
"42",118.095389529608
"43",97.787914705318
"44",83.4191809696678
"45",22.9497113255919
"46",164.196125438034
"47",201.514570499276
"48",NA
"49",273.786029849166
"50",350.184410972518
"51",271.772897276746
"52",NA
"53",NA
"54",57.3742783139797
"55",135.005703137939
"56",73.4793388933424
"57",NA
"58",NA
"59",90.5909657589153
"60",NA
"61",NA
"62",204.685254300838
"63",NA
"64",NA
"65",296.182129717343
"66",237.901941745774
"67",26.2713800700854
"68",290.94798502905
"69",196.657888168312
"70",117.403375207839
"71",62.6210050808502
"72",238.367478653146
"73",423.550511158662
"74",NA
"75",NA
"76",86.4011335925654
"77",189.901311972126
"78",NA
"79",42.2757840208271
"80",NA
"81",269.961077961568
"82",330.355055134178
"83",NA
"84",NA
"85",109.489247782511
"86",398.85189091078
"87",83.9098820341953
"88",19.099595280838
"89",335.299812015185
"90",100.895688113992
"91",1118.39580060812
"92",849.554527639962
"93",719.97170036898
"94",194.22954700283
"95",769.06697097888
"96",616.307954967909
"97",87.5712669002847
"98",209.491608317492
"99",346.485279870696
"100",209.038653488697
"101",38.4885783689614
"102",482.396892666224
"103",46.6543473658414
"104",86.0614174709695
"105",94.6675592180665
"106",364.691547572522
"107",NA
"108",22.1444582966237
"109",65.2632415821519
"110",NA
"111",256.649238826438
"112",358.488582833752
"113",167.480047946795
"114",133.395197080003
"115",206.497073616016
"116",309.82110289549
"117",NA
"118",58.8841277432949
"119",135.420911731001
"120",210.850472803875
"121",221.130031001797
"122",124.71356286144
"123",101.575120357184
"124",122.952071860572
"125",NA
"126",168.29788305434
"127",58.3305162858794
"128",NA
"129",NA
"130",181.559393875159
"131",331.311293106077
"132",NA
"133",109.86671013984
"134",121.039595916773
"135",NA
"136",85.2561644420014
"137",312.664652654034
"138",65.5274652322821
"139",101.91483647878
"140",167.593286653993
"141",280.265800316644
"142",216.462079849497
"143",548.087924920015
"144",NA
"145",131.998586357886
"146",286.619749998346
"147",102.141313893177
"148",117.264972343485
"149",NA
"150",NA
"151",NA
"152",NA
"153",113.125468491445
"154",NA
"155",352.411438880758
"156",211.378920104136
"157",79.3174233533614
"158",305.392211236165
"159",379.777459787097
"160",306.763657801127
"161",138.264461489544
"162",57.6510840426875
"163",NA
"164",114.05654230619
"165",NA
"166",73.3535181075661
"167",30.1969885863051
"168",NA
"169",74.727481088243
"170",27.4490626249513
"171",60.5147651269554
"172",124.529864514207
"173",324.04388451964
"174",NA
"175",NA
"176",76.7305479978012
"177",141.76479574984
"178",NA
"179",NA
"180",NA
"181",249.359182498561
"182",476.357494948963
"183",227.926869849431
"184",17.5796801886606
"185",NA
"186",170.49723038971
"187",282.643813167816
"188",NA
"189",641.691040290414
"190",NA
"191",63.4136760312407
"192",275.109664515532
"193",75.618292251539
"194",14.305823342762
"195",38.7276378619363
"196",25.9694101842224
"197",NA
"198",478.408373757116
"199",NA
"200",869.069331513861
"201",NA
"202",182.314318589817
"203",NA
"204",892.29584856816
"205",305.744509436339
"206",672.335950874083
"207",445.707551533863
"208",12.5317502633166
"209",38.4759962903837
"210",NA
"211",212.57421756901
"212",331.56293467763
"213",NA
"214",NA
"215",NA
"216",NA
"217",34.9152680529153
"218",328.417415033223
"219",34.1980895739905
"220",119.856880530476
"221",396.637445081118
"222",NA
"223",252.358750031467
"224",9.66303634761762
"225",NA
"226",NA
"227",NA
"228",168.297883054341
"229",48.8184648811932
"230",78.6128269530143
"231",86.3382231996773
"232",1503.6967928908
"233",NA
"234",723.909890963777
"235",NA
"236",15.6017774362576
"237",166.007944753212
"238",NA
"239",NA
"240",NA
"241",150.632644731352
"242",733.182882875488
"243",NA
"244",NA
"245",330.719935412929
"246",417.133651084072
"247",NA
"248",61.7150954232608
"249",NA
"250",NA
"251",NA
"252",141.800025569858
"253",183.195064090251
"254",NA
"255",NA
"256",NA
"257",NA
"258",NA
"259",NA
"260",NA
"261",864.388798282983
"262",285.05957225472
"263",246.118039056964
"264",NA
"265",52.3414468829288
"266",182.767273418611
"267",432.571861498821
"268",NA
"269",169.883224955121
"270",NA
"271",69.8934464987186
"272",NA
"273",NA
"274",NA
"275",NA
"276",NA
"277",NA
"278",NA
"279",NA
"280",NA
"281",NA
"282",NA
"283",NA
"284",NA
"285",NA
"286",NA
"287",NA
"288",246.357098549939
"289",NA
"290",NA
"291",NA
"292",NA
"293",NA
"294",NA
"295",NA
"296",NA
"297",162.258485337079
"298",111.980499340881
"299",NA
"300",NA
"301",413.635833239492
"302",309.745610424024
"303",337.97979475222
"304",NA
"305",567.489490086716
"306",105.94110162362
"307",NA
"308",801.46582331627
"309",NA
"310",NA
"311",NA
"312",NA
"313",NA
"314",91.949830245299
"315",168.021077325632
"316",671.958488516754
"317",119.227776601595
"318",1396.20809560213
"319",671.279056273563
"320",44.8425280506631
"321",443.669254804288
"322",NA
"323",NA
"324",NA
"325",504.138724448364
"326",NA
"327",183.597690604735
"328",536.487248471443
"329",481.30225182997
"330",NA
"331",NA
"332",NA
"333",NA
"334",NA
"335",91.2200696877966
"336",130.224513278441
"337",66.1062408468529
"338",NA
"339",23.2516812114549
"340",187.774940692507
"341",119.215194523017
"342",246.105456978387
"343",352.235289780671
"344",139.77431091886
"345",NA
"346",NA
"347",377.877565921875
"348",203.011837850014
"349",NA
"350",128.941141263523
"351",347.756069807037
"352",164.825229366915
"353",178.791336588081
"354",123.807653203851
"355",552.302921243519
"356",542.539228267281
"357",NA
"358",239.059492974915
"359",NA
"360",127.267724812698
"361",245.841233328256
"362",NA
"363",NA
"364",NA
"365",NA
"366",188.680850350096
"367",NA
"368",271.772897276746
"369",44.6915431077315
"370",NA
"371",263.443561258357
"372",493.418793500225
"373",765.493660662834
"374",NA
"375",357.960135533491
"376",NA
"377",NA
"378",28.0957814638413
"379",NA
"380",NA
"381",214.436365198499
"382",NA
"383",NA
"384",361.898326128289
"385",68.8868802125085
"386",NA
"387",273.282746706061
"388",246.860381693044
"389",468.619516623722
"390",237.008614166762
"391",75.8699338230915
"392",NA
"393",NA
"394",NA
"395",NA
"396",268.350571903631
"397",53.9016246265545
"398",243.048011884023
"399",94.6801412966441
"400",485.882128432227
"401",NA
"402",NA
"403",NA
"404",651.638431613887
"405",521.099366371005
"406",149.651242602297
"407",NA
"408",524.924318258604
"409",NA
"410",NA
"411",NA
"412",NA
"413",NA
"414",NA
"415",NA
"416",549.534863956443
"417",269.394884425574
"418",373.889047012767
"419",93.6861570890115
"420",257.11477573381
"421",167.492630025372
"422",NA
"423",NA
"424",NA
"425",301.34078193417
"426",NA
"427",NA
"428",NA
"429",143.850904378011
"430",239.285970389312
"431",380.104593830115
"432",181.181931517831
"433",31.4426143654902
"434",52.6685809259471
"435",NA
"436",NA
"437",17.6652383229885
"438",97.1588107764366
"439",305.744509436339
"440",332.166874449356
"441",34.8775218171824
"442",65.8294351181451
"443",153.249717075498
"444",263.279994236848
"445",15.287225471817
"446",44.880274286396
"447",71.6675195781641
"448",NA
"449",NA
"450",NA
"451",NA
"452",NA
"453",NA
"454",NA
"455",NA
"456",NA
"457",NA
"458",NA
"459",NA
"460",NA
"461",NA
"462",NA
"463",NA
"464",NA
[[alternative HTML version deleted]]
______________________________________________
R-help@r-project.org mailing list
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.
