Displaying 2 results from an estimated 2 matches for "dfit".
Did you mean:
dfi
2012 Feb 10
0
a) t-tests on loess splines; b) linear models, type II SS for unbalanced ANOVA
...test the whether there is a significant difference in
glycerol concentration at any given od as follows:
#test for a difference at od=7
v.pred<-7
pred1<-predict(m1,newdata=data.frame(od=v.pred),se=T)
pred2<-predict(m2,newdata=data.frame(od=v.pred),se=T)
dfit<-data.frame(g=c("Ma","RILI2"),fit=c(pred1$fit,pred2$fit),se=c(pred1$fit,pred2$fit))
alpha<-0.05
dfit$lo<-with(dfit,fit-se*abs(qnorm(alpha/2)))
dfit$hi<-with(dfit,fit+se*abs(qnorm(alpha/2)))
Is this legitimate?
I might also try to tes...
2011 Mar 14
0
Fitting 4 moments distribution w/ Mixture Gaussian
...)
Taking into account that I fix the weight between the two gaussians at
(0.2;0.8) I implemented the below code in R:
distance <-function(parameter,x) {
u=mean(x)
s=sd(x)
sk=skewness(x)
kurt=kurtosis(x)
d1=dnorm(x,parameter[1],parameter[2])
d2=dnorm(x,parameter[3],parameter[4])
dfit=0.2d1+0.8d2
ufit=mean(dfit)
sdfit=sd(dfit)
skfit=skewness(fit)
kurtfit=kurtosis(fit)
abs((u-ufit)/ufit)+abs(s-sdfit)/sdfit)+abs((sk-skfit)/skfit)+abs((kurt-kurtfit)/kurtfit))
}
Parameter<-c(0,0.01,0,0.01) ' starting point of the optimization
opp<-optim(parameter,distance,x...