Displaying 2 results from an estimated 2 matches for "temppow".
Did you mean:
tempo
2012 May 31
1
fitting allometric equation using a for a power model
...lysis of logâtransformed data: Statistical
bias and its correction. Environmental Toxicology and Chemistry. Retrieved
from http://onlinelibrary.wiley.com/doi/10.1002/etc.5620120618/abstract(page
1130)
but that doesn't seem to help much - if anything the fit looks worse. Data
and code below:
temppow<-lm(log(y)~log(x))
plot(log(y)~log(x))
plot(residuals(temppow), main="pow")
abline(temppow)
plot(y~x, main="pow")
tempsum<-summary(temppow)$adj.r.squared
tempint<-summary(temppow)$coefficients[1,1] #intercept of power function
tempslope<-summary(temppow)$coefficients...
2012 Jun 01
1
trouble with append() in a for loop
...7578,502
37121.23438,461
29711.42749,565
33400.49145,510
and the code is:
xy<-read.csv(xy.csv)
x<-xy$x
y<-xy$y
cross.val.error.temp<-numeric(0)
b<-c(1:length(x))
for(i in b){
frame.val.temp<-data.frame(x,y)
frame.val<-frame.val.temp[-i,]
yt<-frame.val$y
xt<-frame.val$x
temppow.val<-lm(log(yt)~log(xt))
tempint.val<-summary(temppow.val)$coefficients[1,1] #intercept of power
function
tempslope.val<-summary(temppow.val)$coefficients[2,1]
nls.model.val<-nls(yt~exp(intercept)*((xt)^slope),
start=list(intercept=tempint.val, slope=tempslope.val))
prediction.cross.val...