Displaying 2 results from an estimated 2 matches for "counts2".
Did you mean:
counts
2010 Feb 05
1
Appending results via for loop
...to assign the size of the dataframe beforehand, or can it be dynamic? I've read the help files and the mailing lists, but I think I am not wording my question properly to find the right answers, so I'm opting for this instead : )
Here's my code:
window <- 50
size <- 0
if (length(counts2$Count)<length(counts1$Count)) { size <-length(counts2$Count) } else { size <- length(counts1$Count) }
a <- 0
b <- 0
alpha <- .05
n <- window
pval <- qchisq((1-(alpha / n)), 1)
results <- c(seq(1,size,window))
loop <- 1
for (i in seq(1,size,window)) {
a <- i
b <...
2010 Feb 09
2
Double Integral Minimization Problem
...869 63.883 44.524 31.918
23.500 17.877 13.915 11.032 8.881 7.245 5.978 4.972 4.175112)
countserr=c(80.085 12.181 4.247 3.148 1.963 1.279 0.448 0.242 0.097
0.055 0.034 0.022 0.015 0.011 0.0082 0.006 0.0043)
intensity=function(x,counts0,index,radius_eff){
thita=x[1]
radiusb=x[2]
counts2=(1-ellip)*(vita-1)/(pi*alpha^2)*counts0*radiusb*exp(-(2*index-0.324)*(radiusb/radius_eff)^(1/index))*(1+(((majaxis_pix^2)+(radiusb^2)-2*majaxis_pix*radiusb*cos(-thita)+(ellip^2-2*ellip)*((radiusb*sin(thita))^2))/(alpha^2)))^(-vita)
return(counts2)
}
sersic=function(p)
{
counts0=p[1]
index=p[2...