Re: Statistical modeling of novel chemical systems( AKA Chemist with lots of stats questions)




Jeremy Boden wrote:
On Tue, 2006-08-29 at 13:06 -0700, ghostwriter wrote:
Let me summerize.

I managed to figure out a good statistical model to one of the major
reactions in my work. Everyone is extremely happy with that fact,
except for the PhD that I bumped from this years call for papers.
Unfortunatly said PhD is also the official go to guy for stats
problems, and while he hasnt snubbed me completely getting my questions
answered has been pulling teeth.

I am trying to iterate in order to maximize the fit of my curve to the
data, but what equation should I use to determine when I have the best
fit? Maximizing for R^2 gives systematic weaknesses in the model due
to the fact that the reaction is exponential, so I don't wish the
points weighted equally thoughout the entire reaction. Ideally I want
to maximize the confidence intervals, and while getting enough repeats
for each time isnt difficult, the curve gives only a single data point
for each time so a T-test isnt possible.

The equation is

f(x)=exp(a*x^b)

a is always negative and b is usally between -1 and -2

f(x) is therefor always between 0 and 1 and approaches 1 as time goes
to infinity.

Any help greatly appreciated.

How were these things done before the introduction of computers?

--
Jeremy Boden

Much slower.

Actually I figured out the model simply by repeatly graphing different
transformations of the data which is exactly how I imagine most
modeling was done prior to computers.

I am not a mathamatican, I am a chemist, if you want to discuss
reaction mechanisms, and energy diagrams of polymer crosslinking I can
reasonably claim some expertise.

Frankly in stats I beat the crap out of 90% of the population but you
guys are part of the 10% that beats the crap out of me. Thats the
simple truth.

My current assumption is that if I take multiple tests at each point in
time and then use that data to generate a hypothetical bell curve at
each of those times, I can then fit the curve to the different bell
curves not the individual points. Then R^2 will become the percent fit
to the therotical and be confidence interval of a sort.

That will likley work for me until I can one of the stats guys from
corporate to sit down with me. Or one of you guys to give me a straight
answer.

Ghostwriter

.



Relevant Pages

  • Re: www.letime.net
    ... >of appreciating and measuring time. ... >is why I have identified and plotted each parameter on a curve. ... >The law for a factor on a stable physical phenomenon is as follows: ... >The experimental results of the curve of reaction of the H.I.V. ...
    (sci.math.num-analysis)
  • Re: Pythonic indentation (or: beating a dead horse)
    ... first time, too. ... my reaction when I first learned of ... So contrary to how you feel about this, one size doesn't fit all - ... alternative was to use a preprocessor? ...
    (comp.lang.ruby)
  • Re: {kjd-sys} Races, classes, and advancement
    ... You were saying this in reaction to how easily a rich person could get +28 in all stats -- but it doesn't feel quite right to me for high plusses to be impossible. ... It'd fit with something like E6, which is meant for a heric-but-not-demigodish campaign. ...
    (rec.games.frp.dnd)
  • Re: Prevent controls being added to UserControl
    ... >I took some time to make a sample that would fit your previous problem, ... >however not any reaction on that. ... Sorry Cor, your example came after i had resolved it with the help of ... Sorry to offend, you help and examples have always proven ...
    (microsoft.public.dotnet.languages.vb)
  • Re: Humbled again
    ... sportbikes came blasting around me like I was standing still. ... reaction was how dangerous and irresponsible it was to pass on a blind ... curve like they did. ...
    (rec.motorcycles)