similar to: Table and Sweave

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)

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))

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

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)

Hi, I have a pentium III 600 with 128 ram running Linux. I have a dataset with this variables in 12000 lines. coleta = response (count data) data = continuous data latitude = continuous data longitude = continuous data variedade = categorical with 14 levels solo = categorical with 22 levels idade = continuous data app = categorical with 4 levels I want to test several models beginning by the

Hi all, I'm from Brazil. I fit a Tobit model to FLUID MILK CONSUMPTION (DEPENDENT VARIABLE) data using survreg (attached). I am confused about the output interpretation and I would like yours explanations. Thanks, Marcio Roberto Silva -------------- next part -------------- A non-text attachment was scrubbed... Name: Tobit model.pdf Type: application/pdf Size: 9018 bytes Desc: not

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

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

I am checking a package with R CMD check and everything is fine until the pdf portion of the checking process and the errors documented at the end of this message happen. I have looked at the pdf that is created and the typesetting is not perfect, but it is acceptable (readable). I am no good at latex (something I plan to remedy in the next large paper I have to write), so is this acceptable for

Hi, I have a problem with Kolmogorov-Smirnov test fit. I try fit distribution to my data. Actualy I create two test: - # First Kolmogorov-Smirnov Tests fit - # Second Kolmogorov-Smirnov Tests fit see below. This two test return difrent result and i don't know which is properly. Which result is properly? The first test return lower D = 0.0234 and lower p-value = 0.00304. The lower 'D'

Hi, I'm building a graph (barplot) in which the X axis label disappears. I tried to use the option mgp of par() and I could not get the desired result. Note that want the axis labels horizontally. caes = c(37,20,19,16,75,103) names(caes) = c("Pinscher", "Pastor \n Alem?o", "Poodle", "Rottweiller", "SRD", "Outros") caess =

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

Hi all I am preparing a document using Sweave; a really useful tool. But I am having a problem. Consider this toy example Sweave file: \documentclass{article} \begin{document} <<echo=TRUE,results=verbatim>>= options(width=40) # Set width to 40 characters hide <- capture.output(example(glm)) # Create an example of the problem, but hide the output summary(glm.D93) #

> In a lattice plot like this: > > win.graph() > xyplot(tmx~frequ|as.factor(as.numeric(spf)),groups=as.factor(blm), > data=tmx,type="l",pch=16,xlab="frequency (N)",ylab="Area held (ha)", > auto.key=list(blm,points=F,lines=T,title="Blm factor",cex.title=0.7, > cex=0.7,corner=c(1,1)),main="Mangroves target=1573ha",layout=c(2,3))

does anyone know to rezize a table produzed by xtable? The size of my table is too big and I would like to resize it like using resizebox but it gives an erros when I try using it using it its fine \SweaveOpts{echo=false} <<results=tex>>= xtable(stats0,caption='Número de transacções dos artigos frequentes e infrequentes',label='tab:INEStats') @ but the size is too

... Thanks, anyway for telling as at least, it reflects your needs. > > > You want NT PtMP and i second that, > not being limited on the asterisk > side is a must in the > telephony ecosystem, since the legacy PABX aren't alwsys easy to > reconfigure. > > _______________________________________________ > -- Bandwidth and Colocation Provided by

Hi, My error isErro em hist.default(dados[[1]], freq = TRUE, seq(0, 30, 0.5), prob = FALSE, : some 'x' not counted; maybe 'breaks' do not span range of 'x' hist(dados[[1]],seq(0, 30, 0.5), prob=TRUE, xlab="chuva (mm/dia)",ylab="frequência", main="", cex.lab=1.6, cex.sub=3,cex.axis=3,cex.main=6) Someone help me? [[alternative