Sorry! I had omitted to copy in an essential line in the code below: On 11-Nov-04 Ted Harding wrote: > On 11-Nov-04 Wei Yang wrote: >> Hi, >> >> I have a list of numbers. For each of the numbers, I take >> sum of squares of the numbers centered on the number chosen. >> If it is less than a certain constant, I will take the >> average of the numbers chosen. >> >> Anyone can give me a sample code. You help will be greatly >> appreciated.
Y<-rep(NA,100) ##### This line is needed! x<-rnorm(100) const <-150 for( i in (1:100) ) {if(sum(x-x[i])^2<const) Y[i]<-x[i]} Y[!is.na(Y)] [1] 0.17096364 -0.32720155 0.19542299 0.13363724 [5] -0.19961480 -0.24486536 -0.31485802 -0.33369635 [9] 0.09981291 0.04263151 0.11127977 0.12144595 [13] -0.27767009 -0.01242218 0.06244776 0.11646301 mean(Y[!is.na(Y)]) [1] -0.04101397 > but I'm sure somebody out there will come up with a much > more elegant solution! > > Best wishes, > Ted. > > > -------------------------------------------------------------------- > E-Mail: (Ted Harding) <[EMAIL PROTECTED]> > Fax-to-email: +44 (0)870 094 0861 [NB: New number!] > Date: 11-Nov-04 Time: 20:56:08 > ------------------------------ XFMail ------------------------------ > > ______________________________________________ > [EMAIL PROTECTED] mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide! > http://www.R-project.org/posting-guide.html -------------------------------------------------------------------- E-Mail: (Ted Harding) <[EMAIL PROTECTED]> Fax-to-email: +44 (0)870 094 0861 [NB: New number!] Date: 11-Nov-04 Time: 21:23:33 ------------------------------ XFMail ------------------------------ ______________________________________________ [EMAIL PROTECTED] mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html