Skip to main content
added 59 characters in body
Source Link
Bob Hanlon
  • 167.4k
  • 7
  • 83
  • 213
confidInterval = Thread[{c1, c2} -> Around@*Interval /@ nlm["ParameterConfidenceIntervals"] nlm["ParameterConfidenceIntervals"]] (* {c1 -> Around[9.297969154154664, 1.0182007406816869`],   c2 -> Around[22.62894607083057, 3.89519370692895]} *)   Plot[nlm[log10q], {log10q, Min[data[[All, 1]]], Max[data[[All, 1]]]}, Epilog -> {Red, AbsolutePointSize[4], Point[data]}, Frame -> True, Axes -> False, GridLines -> Automatic, GridLinesStyle -> Lighter[Gray, .8], FrameLabel -> (Style[#, 14, Bold] & /@ {HoldForm@Log10["q (K/s)"], Row[{Subscript[T, f'], " (", Degree, " C)"}]})] 
Around@*Interval /@ nlm["ParameterConfidenceIntervals"] (* {Around[9.297969154154664, 1.0182007406816869`], Around[22.62894607083057, 3.89519370692895]} *) Plot[nlm[log10q], {log10q, Min[data[[All, 1]]], Max[data[[All, 1]]]}, Epilog -> {Red, AbsolutePointSize[4], Point[data]}, Frame -> True, Axes -> False, GridLines -> Automatic, GridLinesStyle -> Lighter[Gray, .8], FrameLabel -> (Style[#, 14, Bold] & /@ {HoldForm@Log10["q (K/s)"], Row[{Subscript[T, f'], " (", Degree, " C)"}]})] 
confidInterval = Thread[{c1, c2} -> Around@*Interval /@  nlm["ParameterConfidenceIntervals"]] (* {c1 -> Around[9.297969154154664, 1.0182007406816869`],   c2 -> Around[22.62894607083057, 3.89519370692895]} *)   Plot[nlm[log10q], {log10q, Min[data[[All, 1]]], Max[data[[All, 1]]]}, Epilog -> {Red, AbsolutePointSize[4], Point[data]}, Frame -> True, Axes -> False, GridLines -> Automatic, GridLinesStyle -> Lighter[Gray, .8], FrameLabel -> (Style[#, 14, Bold] & /@ {HoldForm@Log10["q (K/s)"], Row[{Subscript[T, f'], " (", Degree, " C)"}]})] 
added 187 characters in body
Source Link
Bob Hanlon
  • 167.4k
  • 7
  • 83
  • 213
Around@*Interval /@ nlm["ParameterConfidenceIntervals"] (* {Around[9.297969154154664, 1.0182007406816869`], Around[22.62894607083057, 3.89519370692895]} *) Plot[nlm[log10q], {log10q, Min[data[[All, 1]]], Max[data[[All, 1]]]}, Epilog -> {Red, AbsolutePointSize[4], Point[data]}, Frame -> True, Axes -> False, GridLines -> Automatic, GridLinesStyle -> Lighter[Gray, .8], FrameLabel -> (Style[#, 14, Bold] & /@ {HoldForm@Log10["q (K/s)"], Row[{Subscript[T, f'], " (", Degree, " C)"}]})] 
Plot[nlm[log10q], {log10q, Min[data[[All, 1]]], Max[data[[All, 1]]]}, Epilog -> {Red, AbsolutePointSize[4], Point[data]}, Frame -> True, Axes -> False, GridLines -> Automatic, GridLinesStyle -> Lighter[Gray, .8], FrameLabel -> (Style[#, 14, Bold] & /@ {HoldForm@Log10["q (K/s)"], Row[{Subscript[T, f'], " (", Degree, " C)"}]})] 
Around@*Interval /@ nlm["ParameterConfidenceIntervals"] (* {Around[9.297969154154664, 1.0182007406816869`], Around[22.62894607083057, 3.89519370692895]} *) Plot[nlm[log10q], {log10q, Min[data[[All, 1]]], Max[data[[All, 1]]]}, Epilog -> {Red, AbsolutePointSize[4], Point[data]}, Frame -> True, Axes -> False, GridLines -> Automatic, GridLinesStyle -> Lighter[Gray, .8], FrameLabel -> (Style[#, 14, Bold] & /@ {HoldForm@Log10["q (K/s)"], Row[{Subscript[T, f'], " (", Degree, " C)"}]})] 
Added ParameterTable
Source Link
Bob Hanlon
  • 167.4k
  • 7
  • 83
  • 213
nlm = NonlinearModelFit[data, {model2, c1 > 5, c2 > 5}, {c1, c2}, log10q]; param = nlm["BestFitParameters"] (* {c1 -> 9.29797, c2 -> 22.629} *) nlm["ParameterTable"] // Quiet 

enter image description here

Plot[nlm[log10q], {log10q, Min[data[[All, 1]]], Max[data[[All, 1]]]}, Epilog -> {Red, AbsolutePointSize[4], Point[data]}, Frame -> True, Axes -> False, GridLines -> Automatic, GridLinesStyle -> Lighter[Gray, .8], FrameLabel -> (Style[#, 14, Bold] & /@ {HoldForm@Log10["q (K/s)"], Row[{Subscript[T, f'], " (", Degree, " C)"}]})] 
nlm = NonlinearModelFit[data, {model2, c1 > 5, c2 > 5}, {c1, c2}, log10q]; param = nlm["BestFitParameters"] (* {c1 -> 9.29797, c2 -> 22.629} *) Plot[nlm[log10q], {log10q, Min[data[[All, 1]]], Max[data[[All, 1]]]}, Epilog -> {Red, AbsolutePointSize[4], Point[data]}, Frame -> True, Axes -> False, GridLines -> Automatic, GridLinesStyle -> Lighter[Gray, .8], FrameLabel -> (Style[#, 14, Bold] & /@ {HoldForm@Log10["q (K/s)"], Row[{Subscript[T, f'], " (", Degree, " C)"}]})] 
nlm = NonlinearModelFit[data, {model2, c1 > 5, c2 > 5}, {c1, c2}, log10q]; param = nlm["BestFitParameters"] (* {c1 -> 9.29797, c2 -> 22.629} *) nlm["ParameterTable"] // Quiet 

enter image description here

Plot[nlm[log10q], {log10q, Min[data[[All, 1]]], Max[data[[All, 1]]]}, Epilog -> {Red, AbsolutePointSize[4], Point[data]}, Frame -> True, Axes -> False, GridLines -> Automatic, GridLinesStyle -> Lighter[Gray, .8], FrameLabel -> (Style[#, 14, Bold] & /@ {HoldForm@Log10["q (K/s)"], Row[{Subscript[T, f'], " (", Degree, " C)"}]})] 
Changed approach to correct error
Source Link
Bob Hanlon
  • 167.4k
  • 7
  • 83
  • 213
Loading
Added "BestFitParameters"
Source Link
Bob Hanlon
  • 167.4k
  • 7
  • 83
  • 213
Loading
Source Link
Bob Hanlon
  • 167.4k
  • 7
  • 83
  • 213
Loading