search for: res4

...is no result value for this. For example- source Basic block : continuation: ; preds = %else, %then %iftmp = phi i32 [ 5, %then ], [ 9, %else ] store i32 %iftmp, i32* %datasize ; 3 instructions below being copied store i8000000 0, i8000000* %res %res4 = load i8000000, i8000000* %res ret i8000000 %res4 After copying the 3 instructions above to destination basic block. destination basic block : else: ; preds = %entry store i8000000 0, i8000000* %res %1---->issue = load i8000000, i8000000* %...

...in (1:length(lalgo))) { current_algo<-lalgo[k] res2<-res[res$instance==current_instance & res$try==current_try & res$idalgo==current_algo,] # res2 contains for a given instance, a given try, and a given algo, all results res3<-res2[res2$best==min(res2$best),] res4<-res3[res3$time==min(res3$time),] if (nrow(res4)>1) { res4<-res4[1,] } if (nrow(res4)==1) { res4$best<-(res4$best*100/optimal_values[optimal_values $instance==linstance[i],]$optimum)-100 print(res4) bestalgo<-rbind(bestalgo,res4) } } } eps...

...method.hclust="complete", method.dist="euclidian", nboot=500) a.res2 <- pvclust(t(s), method.hclust="complete", method.dist="maximum", nboot=500) a.res3 <- pvclust(t(s), method.hclust="complete", method.dist="canberra", nboot=500) a.res4 <- pvclust(t(s), method.hclust="ward", method.dist="euclidian", nboot=500) a.res5 <- pvclust(t(s), method.hclust="ward", method.dist="maximum", nboot=500) a.res6 <- pvclust(t(s), method.hclust="ward", method.dist="canberra", nboot...

Just for the fun of it (actually in order to test some Fortran more extensively), I've compiled and installed g77, yesterday on our Solaris 2.5.1 machines. After doing (in R's toplevel directory) make distclean ./configure now chooses g77 (instead of Sun's f77 which was chosen before). I now do make make check cd tests/Examples and carefully compare base-Ex.Rout to an

...case? data = read.table("data.txt", header=T, nrows=90) attach(data) res1 = pairwise.t.test(all, location, p.adj="bonf") print(res1) res2 = pairwise.t.test(M, location, p.adj="bonf") print(res2) res3 = pairwise.t.test(OA, location, p.adj="bonf") print(res3) res4 = pairwise.t.test(UE, location, p.adj="bonf") print(res4) res1 = t.test(all~job) print(res1) res2 = t.test(M~job) print(res2) res3 = t.test(OA~job) print(res3) res4 = t.test(UE~job) print(res4) I'd also like to compare the four exposure metrics - how to do that? One potential proble...

...0,4, add = TRUE, col = 2) ## If we repeat the procedure, using the fix intial value 1, ## the situation is even worse res3 <- NULL for(i in 1:1000) res3[i] <- arms(1, logDichteGam, function(x) (x>0)&(x<100), 1) ## Generating a sample of thousand observations with 1 call of arms res4 <- arms(1, logDichteGam, function(x) (x>0)&(x<100), 1000) ## Plot of the samples par(mfrow = c(2,2)) plot(res3, log = "y") plot(res4, log = "y") hist(res3, freq = FALSE, xlim = c(0,4), breaks = seq(0,100,by = 0.1), ylim = c(0,1)) curve(DichteGam, 0,4, add = T...

...S.Y) if(toto$Rsquare > R2cum) { R2cum <- toto$Rsquare no.sup <- j } } mm <- 1 FP <- FPval(R2cum,R2prev,n,mm,p) if(FP$pval <= alpha) { adjRsq <- RsquareAdj(R2cum,n,mm) res1 <- var.names[no.sup] res2 <- no.sup res3 <- R2cum res4 <- R2cum res5 <- adjRsq res6 <- FP$Fstat res7 <- FP$pval X.out[no.sup] <- 0 delta <- R2cum } else { stop("Procedure stopped (alpha criterion): pvalue for variable ",no.sup," is ",FP$pval) } ## Add variables X to the model while(...

...s/gilp/R/doc/html/packages.html These are good. DSE documentation is added and R documentation is not clobbered. After help.start() ~/.R has doc/ library/ /home/res/gilp : ls -l .R/doc/html/pac* lrwxrwxrwx 1 gilp res 61 Nov 18 14:41 .R/doc/html/packages-foot.html -> /home/res4/gilp/R-versions/R-0.63/doc/html/packages-foot.html lrwxrwxrwx 1 gilp res 61 Nov 18 14:41 .R/doc/html/packages-head.html -> /home/res4/gilp/R-versions/R-0.63/doc/html/packages-head.html -rw-r--r-- 1 gilp res 1997 Nov 18 14:41 .R/doc/html/packages.html ls -l .R/doc/ht...

...#Second function f<-function(listRes,Toselect){ res2<-split(listRes,gsub("[0-9]","",names(listRes))) res3<-lapply(seq_along(res2),function(i) lapply(res2[[i]],function(x) x[x[["FDR"]]<0.01,c("Seq","Mod","z","spec")])) res4<-lapply(res3,function(x) x[names(x)[names(x)%in%Toselect]]) res4New<- lapply(res4,function(x) lapply(names(x), function(i) do.call(rbind,lapply(x[i],function(x) cbind(folder_name=i,x))) )) library(plyr) library(data.table) res5<-lapply(res4New,function(x) lapply(x,function(x1){ x1<- d...

Thanks. Yes. Your approach can identify: Glaxy ace S 5830 and S 5830 Glaxy ace But you can not identify using same program: Iphone 4S 16 G Iphone 4S 16G How should I solve both in same time. Kind regards,Tammy [[alternative HTML version deleted]]

...ices (not only in terms of sign and row order). I'm not a specialist, so would appreciate any advice on whether this should really be the case: > res3 = > fastICA(af[,2:20],4,alg.typ="parallel",fun="logcosh",alpha=1,method="C",row.norm=TRUE) colstandard > res4 = > fastICA(af[,2:20],4,alg.typ="parallel",fun="logcosh",alpha=1,method="C",row.norm=TRUE) colstandard > res3$A [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [1,] 0.34926...

...mva.so dblcen.o hclust.o kmns.o mkdir ../../../../library/mva/libs mkdir ../../../library/stepfun mkdir ../../../library/stepfun/R Building system startup profile You should `make docs' now ... make: Fatal error: Don't know how to make target `*/man/*.Rd' Current working directory /home/res4/gilp/R-versions/R-0.63.2/src/library [17] /home/res/gilp/R-versions/R-0.63.2 : which perl /usr/bin/perl [18] /home/res/gilp/R-versions/R-0.63.2 : perl --version This is perl, version 5.004_04 built for sun4-solaris Copyright 1987-1997, Larry Wall Perl may be copied only under the terms of either...

...= NC60.DATA$Qe+NC60.DATA$sdQe ,ymin=NC60.DATA$Qe-NC60.DATA$sdQe) ,colour=alpha("black",0.4) ,width=1) + geom_errorbarh( aes(xmax = NC60.DATA$Ce+NC60.DATA$sdCe ,xmin=NC60.DATA$Ce-NC60.DATA$sdCe) ,colour=alpha("black",0.4) ,width=1) + geom_line(data=NC60.Res4.curve ,aes(x,y) ,size=1 ,colour=alpha("black",0.5)) + xlab(C[e]~(mmol/m^3)) + ylab(q[e]~(mmol/m^3)) print(NC60.iso.graph) dev.off() Regards/Cordialement ------------- Benoit Boulinguiez Ph.D student Ecole de Chimie de Rennes (ENSCR) Bureau 1.20 Equipe CIP UMR CNRS 6226 &quot...

...(y-y1,n-n1,p1H, log=FALSE)})? >>>> >>>>#term2_p0<- term1_p0*dbinom(x-x1,m-m1, p0L, log=FALSE)* dbinom(y-y1,n-n1,p0H, log=FALSE); >>>>#term2_p1<- term1_p1*dbinom(x-x1,m-m1, p1L, log=FALSE)* dbinom(y-y1,n-n1,p1H, log=FALSE)})? >>>> >>>>res4<-do.call(rbind,lapply(seq_len(nrow(res3)),function(i){ >>>>Pm2<-rbeta(1000,0.2+res3[i,"x"],0.8+res3[i,"m"]-res3[i,"x"]); >>>>Pn2<- rbeta(1000,0.2+res3[i,"y"],0.8+res3[i,"n"]-res3[i,"y"]); >>>>Fm...

Hi, #working directory data1 #changed name data to data1.? Added some files in each of sub directories a1, a2, etc. ?indx1<- indx[indx!=""] lapply(indx1,function(x) list.files(x)) #[[1]] #[1] "a1.txt"??????? "mmmmm11kk.txt" #[[2]] #[1] "a2.txt"??????? "mmmmm11kk.txt" #[[3]] #[1] "a3.txt"??????? "mmmmm11kk.txt" #[[4]] #[1]

Dear all, Without a loop, I would like transform 3 numeric vectors empty of 0/1 of same length Vec1 : transform 1 to A and 0 to "" Vec2 : transform 1 to B and 0 to "" Vec3 : transform 1 to C and 0 to "" to obtain only 1 vector Vec who is the paste of the 3 vectors (Ex : ABC, BC, AC, AB,...) Any idea ? Thank you for your help -- Michel ARNAUD

...tml2dos creating tests/Makefile creating tests/Examples/Makefile creating src/include/Platform.h R is now configured for sparc-sun-solaris2.6 Source directory: . Installation directory: /usr/local C compiler: gcc -g -O2 FORTRAN compiler: g77 -O2 [15] /home/res4/gilp/R : make make: Fatal error: Don't know how to make target `src/*/*.[chfy]' [16] /home/res4/gilp/R : make install make: Fatal error: Don't know how to make target `src/*/*.[chfy]' ... [31] /home/res4/gilp/R-versions/R-0.63 : make -d MAKEFLAGS value: Building all because it is...

...") res3 <- nlme(logconc~p2+p3+log(dose/(exp(p1)-exp(p2))* (exp(-exp(p2)*tm)-exp(-exp(p1)*tm))),start=list(fixed=c(5,-2,-0.1,0)), fixed=list(p1+p2~1,p3~bilirubin),random=pdDiag(p1+p2+p3~1), groups=~subj,data=datam,verbose=T,method="ML") However, when I add a second covariate, res4 <- nlme(logconc~p2+p3+log(dose/(exp(p1)-exp(p2))* (exp(-exp(p2)*tm)-exp(-exp(p1)*tm))), start=list(fixed=c(5,-2,-0.1,0,0)), fixed=list(p1+p2~1,p3~bilirubin+age),random=pdDiag(p1+p2+p3~1), groups=~subj,data=datam,verbose=T,method="ML") I get the error: Error in fixed[[nm]][[3]] !=...

Thank you very much for your reply. Your code work well with this example. I modified a little to fit my real data, I got an error massage. Error in split.default(x = seq_len(nrow(x)), f = f, drop = drop, ...) : Group length is 0 but data length > 0 On Thu, Jan 31, 2013 at 12:21 PM, arun kirshna [via R] < ml-node+s789695n4657196h87@n4.nabble.com> wrote: > Hi, > Try this: >

I have a large dataset (~1 million rows) of three variables: ID (patient's name), DATE (of appointment) and DIAGNOSIS (given on that date). Patients may have been assigned more than one diagnosis at any one appointment - leading to two rows, same ID and DATE but different DIAGNOSIS. The diagnoses may change between appointments. I want to subset the data in two ways: - define groups