search for: ssum

...new<-ttab+pseudo/(ncol(ttab)*nrow(ttab)) ttabnew<-prop.table(ttabnew) pah<-(margin.table(ttabnew,1)) pbh<-(margin.table(ttabnew,2)) s<- foreach (ii =1:4) %dopar% { for (i in 1:ncol(ttabnew)){ for (j in 1:ncol(ttabnew)){ s<-rep(NA,2) s[ii]<-d[i,j]*pah[i]*pbh[j] } ssum<-sum(as.vector(unlist(s))) pbhh<-NA pahh<-NA pbhh[j]<- d[i,j]*pbh[j] pahh[i]<-d[i,j]*pah[i] sumpbh<-sum(pbhh) sumpah<-sum(pahh) pah[i]<-pa[i]/(1-alphanew+(alphanew*sumpbh/ssum) ) pbh[j]<-pb[j]/(1-alphanew+(alphanew*sumpah/ssum) ) # s<-sum(d[i,j]*pah[i]*pb...

Hi there, Is there an implementation of the Gaveh Stehfest algorithm in R somewhere ? Or some other inversion ? Thanks, Tolga

See attached file. This was on an old DELL with removable drives. -------------- next part -------------- An embedded and charset-unspecified text was scrubbed... Name: anacdump.txt URL: <http://lists.centos.org/pipermail/centos/attachments/20050418/1d7e0afc/attachment-0001.txt>

...2 line in - let cmd = sprintf "sha1sum %s" (quote (mf_folder // disk)) in + let sha_num = Str.matched_group 1 line in + let disk = Str.matched_group 2 line in + let expected = Str.matched_group 3 line in + let sha = sprintf "sha%ssum" sha_num in + let cmd = sprintf "%s %s" sha (quote (mf_folder // disk)) in let out = external_command cmd in match out with | [] -> - error (f_"no output from sha1sum command, see previous errors") +...