search for: propor

...The purpose of my analysis is to compare the pattern of the whole genome 'gene expression' between the two tissue types. I have 'age' and 'sex' as covariates. Since 'age' and 'sex' definitely have influence on gene expression, I need to first filter out the proporation of 'gene expression' attributable to 'age' and 'sex', then compare the 'remaining gene expression' value between 'liver' and 'brain'. How to do the covariate adjustment to remove the effects of these two covariates? Should I do a 'step-wise...

Hi, I have a large matrix (SNPs) that I want to cycle over with logistic regression with interaction terms. I have made a loop but I am struggling to identify to the formula the name of the column in a way which is meaningful to the formula. It errors becasue it is not evaluated proporly. (below is a pilot with only 7 to 33 columns, my actual has 200,000 columns) My attempts: for (i in 7:33) { label <- colnames(n)[i] model1 <- glm(AS~label*interaction,family=binomial("logit"),data=n) X <- summary(model1)$coefficients[2,1] Y <- c(label,X) vector &lt...

...ite bitmap.c/bitmap.h, obstack.c and obstack.h, into >C++. They > are currently just modified versions of what gcc is >using. > You can get rid of the obstacks, of course. > > Using set<u32> or BitVector or something like it will >result in a > slowdown of mammoth proporations, and about a 10x memory >increase. > Trust me here, you don't want to use a non-sparse >bitmap. > BDD's are acceptable, but are about a 2x slowdown. With >various > optimizations that are going to be published in the next >few months, > you can get bitmaps...

...e an error. The subset only work with length < 4, it is correct? > plot((ocorrencia/isca)~frag,subset=especieama==c("grupo1","grupo2","grupo3","AnoplotermesSp1","NeocapritermesOpacus"),pch=c(1,2,3,4,5),xlab="?rea (ha)",ylab="Propor??o de iscas ocupadas por t?rmitas") Warning messages: 1: longer object length is not a multiple of shorter object length in: is.na(e1) | is.na(e2) 2: longer object length is not a multiple of shorter object length in: "==.default"(especieama, c("grupo1", "grupo2&q...

I am having an issue with Builder.node where it does not produce an error in firefox 1.5 however it just hangs and stops executing. For example if i run the following code, alert start pops up however Builder.node never pops up. alert(''start''); alert(Builder.node(''b'', ''test'')); Anyone have any ideas

...;mmout","rb") >readBin(mm5file,integer(),n=1,size=1) This gives output of "0", which is expect for this file. However, when I started to read the rest of the file, I cannot get the result as expected. The next bit of data is an 50 by 20 integer array, which I cannot read proporly. Here is the command I tried. >readBin( mm5file,integer(),n=1) Another way I have tried uses .Fortran (). The fortran code were. 1.f subroutine readmm5(mm5file,IFLAG,MIF) integer iflag,MIF(50,20),MRF(20,20) CHARACTER*80 mm5file,MIFC(50,20),MRFC(20,20) OPEN (11,FILE=mm5f...

...;)? Please, feel free to comment/suggest. Regards, Eduardo Esteves # Dados CO2<-c(141,172,181,227,309,414,641,936) Prop<-c(0.25,0.34,0.34,0.68,0.85,0.99,0.98,0.99) # Diagrama dispersão plot(Prop~CO2, las=1, xlim=c(0,1000),ylim=c(0,1),pch=16,cex=1.5, xlab="CO2 (ppm)", ylab="Proporção de respostas correctas") # Estimaçao (Método Mínimos Quadrados) ajuste<-nls(Prop~(1/3+exp(B*(T-log(CO2))))/(1+exp(B*(T-log(CO2)))), data=data.frame(CO2=CO2,Prop=Prop),start=list(B=-10,T=5)) summary(ajuste) # Ilustracao do ajuste PropEsp<-predict(ajuste,newdata=list(CO2=seq(0,1000...

..., however, is that someone needs to rewrite bitmap.c/bitmap.h, obstack.c and obstack.h, into C++. They are currently just modified versions of what gcc is using. You can get rid of the obstacks, of course. Using set<u32> or BitVector or something like it will result in a slowdown of mammoth proporations, and about a 10x memory increase. Trust me here, you don't want to use a non-sparse bitmap. BDD's are acceptable, but are about a 2x slowdown. With various optimizations that are going to be published in the next few months, you can get bitmaps to use just as little memory as BDDs wo...