similar to: big model on modest computer.

Hi, I make an analysis and depending of the order of the variables, the significance change, look. m1 <- glm((infec/ntot)~idade+sexo+peso,family=binomial,weights=ntot) > anova(m1,test="F") Analysis of Deviance Table Model: binomial, link: logit Response: (infec/ntot) Terms added sequentially (first to last) Df Deviance Resid. Df Resid. Dev F Pr(>F)

Hello R users, I'm analyzing an experiment in a balanced incomplet block design (BIB). The effect of blocks are assumed to be random, so I'm using nlme::lme for this. I'm analysing another more complex experiments and I notice some diferences from doBy::popMeans() compared multcomp::glht() and contrast::contrast(). In my example, glht() and contrast() were equal I suspect popMeans()

Hi, in Latex I get the table using: \begin{table}[H] \centering \renewcommand{\arraystretch}{1.3} \setlength{\tabcolsep}{18pt} \begin{tabular}{cc} \hline Idade & Frequ?ncia \\ \hline $18 \vdash 26$ & 11 \\ $26 \vdash 34$ & 8 \\ $34 \vdash 42$ & 26 \\ $42 \vdash 50$ & 20 \\ $50 \vdash 58$ & 23 \\ $58 \vdash 66$ & 30 \\ $66

Hi, I'm trying to build a thesaurus that will sensible values for rare words. I suspect the best algorithm to use is k-means although I'm not sure about that -- I would have preferred a k dimensional space with a binary cluster in each dimension so a word can belong to 0..k clusters, but I digress... I can measure the strength of correlation between words fairly easily by counting

Hi, I have some simple statistics to calculate for a large number of variables. I created a simple function to apply to variables. I would like the variable name to be placed automatically. I tried the following function but is not working. desc = function(x){ media = mean(x, na.rm=T) desvio = sd(x, na.rm=T) cv = desvio/media*100

Dear R list, I'm having some problems with font encodings when using R+Sweave+Latex in my native language: Portuguese. My environment: Kubuntu 5.10 Linux $> uname -a Linux nassa 2.6.12-9-686 #1 Mon Oct 10 13:25:32 BST 2005 i686 GNU/Linux R> R.version _ platform i486-pc-linux-gnu arch i486 os linux-gnu system i486, linux-gnu

Hello I wonder how to change the font of chart to Times New Roman and size 9. plot(c(0,100,20),c(0,600,50), xlab= 'Idade(meses)', ylab="Peso(kg)", type = "n", axes=F) axis(1, pos=0, at=seq(0,100,20)) axis(2, pos=0, at=seq(0,600,100)) t<- seq(0,100,1) TA=543.56*(1-0.8976*exp(-0.0522*t)) NI=498.97*(1-0.9259*exp(-0.0494*t)) RC=514.57*(1-0.9112*exp(-0.0499*t))

Looking for cheap high-quality software? that will be accessing this feed. Poets are the unacknowledged legislators of the world. I'd give my right arm to be ambidextrous. By all means, let's not confuse ourselves with the facts! Life is anything that dies when you stomp on it. Big thinking precedes great achievement. To be idle requires a strong sense of

Hi I have tryed to do the meta-regression in metafor package, but I would like to get the standardized coefficients for each variable, however in command:   Ø  res<-rma.uni (yi, vi, method="REML", mods=~cota+DL+uso+gadiente+idade, data= turbidez)   I just have the coefficients no standardized (estimate) of the multiple regression. What I need to do? Thanks Fernanda Melo

-----BEGIN PGP SIGNED MESSAGE----- Hash: SHA1 Hi all, this is probably simple and I'm just doing something stupid, sorry about that :-) I'm trying to convert words (strings of letters) into a fairly small dimensional space (say 10, but anything between about 5 and 50 would be ok), which I will call a feature vector. The the distance between two words represents the similarity of the

Hi, I'm using kinship2 to calculate heritabilty, but I would like calculate in GenABEL too. I trying the code: > require(kinship2) > require(GenABEL) > pedig = with(Dados, pedigree(id=IID, dadid=PAT, momid=MAT, sex=SEX, famid=FID, missid=0)) > kmat = kinship(pedig) > (mod1 = polygenic(altura ~ SEX + idade, data=Dados, kin=kmat)) Erro em intI(i, n = d[1],

Hi, I have a experiment with block, plots, sub-plots, and sub-sub-plots with repeated measures and 3 factors (factorial design) when we have been observed diameter (mm), high (cm) and leaves number (count). However, we don't have one treatment in one factor, so, my design is unbalanced. On a previous message here, a friend tell me that "It appears to me that your design is a split-split

Hello R helpers, I would like to automate this code for many files of the same type. But I don´t know how to make it. In particular, i don´t know how to read many files each one as an r object with the name of the file. Then a for loop would be sufficient, right? Many thanks and a happy new year. Dominic datos <- read.table('global2001.csv',head=T,sep=';',stringsAsFactors=F)

Dear R users, I am running R2.1 on windows XP SP2. I have a problem that I fear I am providing insufficient information to resolve. But if anyone can suggest directions to understand the problem I would be most grateful. I import modest datasets which create a list of data frames (primary data object) amounting to 15.9Mb when measured internally by object.size() . The sizes of the list

Taking David's example further, here is a first try: id | name | description | more info -----+---------+-------------------+----------------- 6 | Inset | An inset element | just one element 8 | Stories | Another element | another element : with 2 lines, without : colons on the left. 9 | Other

Dear list member, My question is related to input file format to an Anova from car package. Here is an example of what I did: My file format is like this (and I dislike the idea that I will need to recode it): Hormone day Block Treatment Plant Diameter High N.Leaves SH 23 1 1 1 3.19 25.3 2 SH 23 1 1 2 3.42 5.5 1 SH 23 1 2 1 2.19 5.2 2 SH 23 1 2 2 2.17 7.6 2 CH 23 1 1 1 3.64 6.5 2 CH 23 1 1 2

https://bugzilla.mindrot.org/show_bug.cgi?id=1585 emlyn.corrin at gmail.com changed: What |Removed |Added ---------------------------------------------------------------------------- CC| |emlyn.corrin at gmail.com --- Comment #13 from emlyn.corrin at gmail.com --- +1 Any progress on this? -- You are receiving this

Howdy, I've been thinking a lot about definition list syntax, and what I did and didn't like about the PHP Markdown Extra syntax (which seems to be a de-facto standard, discount notwithstanding, yes?). What I came up with is a single character change to the PHP Markdown Extra syntax. I just published a detailed explanation of my thoughts and reasoning for this on [my

Buenas tardes. Tengo un ensayo de dos factores (Variedad y dosis) Quisiera hacer anovas para cada una de las dos variedades que estoy evaluando. Lo puedo hacer cogiendo los datos primero de una variedad y luego de la otra. Pero creo que hay un comando "by variedad" que lo hace si necesidad de andar modificando el archivo de datos. Pero no sé dónde hay que colocarlo. Si alguno me puede

msz De: Fern?ndez-L?pez, Ambito Iberoamericano, Paseo de la Castellana 223, Madrid. [1]InEnglish - [2]EnEspanol Caros amigos luso-brasileiros, ? de se perguntar se as "patrulhas ideol?gicas" esquerdistas estar?o impedindo que os ?ltimos artigos do ex preso pol?tico e escritor cubano Armando Valladares - que abordam delicados aspetos das rela??es entre o regime