search for: walmes

...tes #------------------------------------------------------------------------------------------ # draw an arrow with the coordinates extracted arrows(id$x[1], id$y[1], id$x[2], id$y[2]) #------------------------------------------------------------------------------------------ At your disposal. Walmes. ============================================================ Walmes Marques Zeviani LEG (Laboratório de Estatística e Geoinformação) Departamento de Estatística - Universidade Federal do Paraná fone: (+55) 41 3361 3573 VoIP: (3361 3600) 1053 1173 e-mail: walmes@ufpr.br / @walmeszevian...

...v(y~cult*fert, subset(da, cult!="A")) summary(m0) coef(m0) # adjust incomplete factorial m1 <- aov(y~cult*fert, da) summary(m1) coef(m1) require(multcomp) glht(m0, linfct=matrix(c(1,1,10,0), nrow=1)) # work glht(m1, linfct=matrix(c(1,1,0,10,0), nrow=1)) # don't work Thank you. Walmes. ========================================================================== Walmes Marques Zeviani LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) Departamento de Estatística - Universidade Federal do Paraná fone: (+55) 41 3361 3573 VoIP: (3361 3600) 1053 1173 e-mail: wa...

...plot result like this plot(1, main="aaaaaaaaaaaaaaaaaaaa bbbbbbbbbbbbbbbbbbbb") I remember once I saw a meta character like "\n" that avoid this breack line plot(1, main="aaaaaaaaaaaaaaaaaaaa\(?) bbbbbbbbbbbbbbbbbbbb") Does someone know that? Thanks. Walmes ========================================================================== Walmes Marques Zeviani LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) Departamento de Estatística - Universidade Federal do Paraná fone: (+55) 41 3361 3573 VoIP: (3361 3600) 1053 1173 e-mail: wal...

...omized plot, I want semitransparent colors # and use the colors of trellis.par.set("superpose.polygon") to fill densityplot(~y|B, groups=A, data=dt, plot.points="rug", col=2:3, panel=panel.superpose, panel.groups=my.panel.densityplot) Thanks! Walmes. ========================================================================== Walmes Marques Zeviani LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) Departamento de Estatística - Universidade Federal do Paraná fone: (+55) 41 3361 3573 VoIP: (3361 3600) 1053 1173 e-mail: wa...

...space) #----------------------------------------------------------------------------- 1 & a & \textcolor{red}{ 2.00 } \\ 2 & b & \textcolor{red}{ 4.00} \\ #----------------------------------------------------------------------------- Any contribution will be useful. Thanks. Walmes. ========================================================================== Walmes Marques Zeviani LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) Departamento de Estatística - Universidade Federal do Paraná fone: (+55) 41 3361 3573 VoIP: (3361 3600) 1053 1173 e-mail: wa...

...t;) # copy this result curve(3 * x^2) # paste inside # my failed attempts curve(as.expression(D(expression(x^3), "x"))) curve(as.character(as.expression(D(expression(x^3), "x")))) curve(noquote(as.character(as.expression(D(expression(x^3), "x"))))) Thanks in advance. Walmes Zeviani, Brasil. _________________________________________________________________ Quer deixar seus vídeos mais divertidos? Com o Movie Maker isso fica fácil. ndows Live:Dicas - Movie Maker:Hotmail:Tagline:1x1:Titulo Legendas Creditos [[alternative HTML version deleted]]

...ultC:fert NA NA NA NA ... > require(gmodels) > estimable(m0, cm=c(1,0,0,100,0,0)) Erro em estimable.default(m0, cm = c(1, 0, 0, 100, 0, 0)) : Dimension of structure(c(1, 0, 0, 100, 0, 0), .Dim = c(1L, 6L)): 1x6, not compatible with no of parameters in m0: 5 Thanks. Walmes. ========================================================================== Walmes Marques Zeviani LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) Departamento de Estatística - Universidade Federal do Paraná fone: (+55) 41 3361 3573 VoIP: (3361 3600) 1053 1173 e-mail: wa...

I am trying to calculate Poisson regression based adjusted rate ratios in R, but R's default in glm does not code the intercept as the global rate. In SAS I use "cell means" coding so that the intercept is the global rate, but I do not know how to do this in R. If anyone knows a way to make glm use "cell means", or an how to find adjusted rate ratios I would be grateful.

...da, start = list(A = c(10, : fator de passos 0.000488281 reduzido abaixo de 'minFactor' de 0.000976562 #------------------------------------------------------------ That indicates the categorical variable was accounted. But values of B didn't change during iteration. Why? Thanks, Walmes. ========================================================================== Walmes Marques Zeviani LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) Departamento de Estatística - Universidade Federal do Paraná fone: (+55) 41 3361 3573 VoIP: (3361 3600) 1053 1173 e-mail: wa...

Dear all, I have been working in a plot based on figure 5.6 of the Lattice book (http://lmdvr.r-forge.r-project.org/figures/figures.html). I have already modified it to include the size of the circles as another variable, but I would like to modify the legend to show it (like they do it in http://www.jstatsoft.org/v15/i05/paper). I have divided my variable in intervals: DATA$s_Shape_2 <-

Hi I ran the following sweave file in R2.14.1 and upgraded to R2.15 yesterday with the same setup I got the following error message when I rand the following Sweave file > setwd("D:/Cic/Sweave/Parasite/Comb/12") > Sweave("D:/Cic/Sweave/Parasite/Comb/12/ParasiteComb12.Rnw") Error: c("'ParasiteComb12DS.Rnw'", "'ParasiteComb12.Rnw'")

...ing something diferent. This a proprosital difference? Does popMeans account the variance of random effects or something similar? The code below is reproducible, for easy evaluation the principal results are appended. Thank you. # reading BIB data da <- read.table("http://www.leg.ufpr.br/~walmes/data/pimentel_bib3.txt", + header=TRUE, sep="\t", colClasses=c("factor","factor",NA)) str(da) ## 'data.frame': 30 obs. of 3 variables: ## $ bloco : Factor w/ 10 levels "1","10","2","3",..: 1...

When I use barchart (with default formatting options), I get bars whose lengths/heights are not proportional to their value. For example: http://drop.io/wbagm6s/asset/capture-png Many of the values in this chart are 1; however, because the blue bars extend to the left of the "0" tick mark, those bars appear to represent higher numeric values. Is there a way to make the length of the

Hello, I want to do with Lattice functions (qqmath, histogram) a figure like this below. n <- 1000 x <- rnorm(n) qqnorm(x); qqline(x) op <- par(fig=c(.02,.5,.5,.98), new=TRUE) hist(x, xlab="", ylab="", main="", axes=FALSE) box() par(op) Is possible? Thanks. Walmes Zeviani. _________________________________________________________________ CANSADO DE ENTRAR EM TODAS AS SUAS DIFERENTES CONTAS DE EMAIL? JUNTE TODAS AGORA. [[elided Hotmail spam]] [[alternative HTML version deleted]]

...LSE, panel=function(x, y, subscripts, groups, ...){ panel.barchart(x, y, subscripts=subscripts, groups=groups, ...) d <- 0.22 # <------ how obtain "d" or coordinates? panel.text(x+c(-d,0,d), y, label=y, pos=3) }) Thanks in advance. Walmes Zeviani, Lavras - MG, Brasil. _________________________________________________________________ os. dium=Tagline&utm_campaign=InfuseSocial [[alternative HTML version deleted]]

----- Original Message ----- From: "Silvano" <silvano at uel.br> To: <r-help at r-project.org> Sent: Thursday, June 30, 2011 9:07 AM Subject: Tables and merge > Hi, > > I have 21 files which is common variable CODE. > Each file refers to a question. > > I would like to join the 21 files into one, to construct > tables for each question by CODE. >

Hi, I was wondering: I've got a dataset where I've got student 'project's nested within 'school's, and 'division' (elementary, junior, or senior) at the student project level. (Division is at the student level and not nested within schools because some students are registered as juniors & others as seniors within the same school.) So schools are random,

Dear useRs, How could I obtain the confidence intervals for the means of my treatments, when my data was fitted to a GLM? I need the CI's for the Poisson and Negative Binomial distributions. Here's what I have: mydata1 <- data.frame('treatments'=gl(4,20), 'value'=rpois(80, 1)) model1 <- glm(value ~ treatments, data=mydata1, family=poisson) means1 <-

Hi, How can I find and download a function in R to do the LU decompostion for finding the upper and lower triangular matrix.  Thank you so much. Joe ___________________________________________________ 您的生活即時通 － 溝通、娛樂、生活、工作一次搞定！ [[alternative HTML version deleted]]

I am using levelplot() function from pkg:lattice and I want to add some points to the plot generated by levelplot() similarly as in basic R graphic sequence: plot(...); points(...), but it does not work. I was reading documentation for lattice, but so far without much success. I would very much appreciate some tips. Best regards, Ryszard