Skip to main content
edited tags
Link
smillig
  • 2.4k
  • 30
  • 31
Source Link

How do I fit a non-linear regression to only one factor level on my scatter plot?

My data is structured like so:

'data.frame': 50 obs. of 3 variables:

  • Project : Factor w/ 2 levels "A","B": 1 1 1 1 1
  • x: int 2 2 2 2 6 4 4 4 6 4 ...
  • y: num 0.622 0.425 0.363 0.344 0.346 ...

I have 'attached' the data and plotted both levels in my scatter plot using:

plot(x,y,pch=as.numeric(Project))
TSF<-x[order(x)]
SD<-y[order(x)]

And fitted a non-linear regression to the data

nls_fit <- nls(SD ~ a - (b*TSF)+ ((c*TSF)^(2)), start = list(a = 0.34, b = 0.017, 
+     c = 0.0003))
lines(TSF, predict(nls_fit), col = "red")

This works well...

But how can I fit this equation to only the factor "A" of the data? I'm very new at this so if you have the time an inclination it would be great if you could describe what any code you write does.

Thanks Kris