search for: walmeszeviani

...====================================================== 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: walmes@ufpr.br twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ========================================================================== [[alternative HTML version deleted]]

...osal. 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 / @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes ============================================================ > Date: Mon, 14 Feb 2011 20:05:54 +0000 > From: Colin.bleay@bristol.ac.uk > To: walmeszeviani@hotmail.com > Subject: ternary contour plot > > Dear Walmes, > > Firs...

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

...====================================================== 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: walmes@ufpr.br twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ========================================================================== [[alternative HTML version deleted]]

...========================================================= 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: walmes@ufpr.br skype: walmeszeviani twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ========================================================================== [[alternative HTML version deleted]]

Hello, I look for a simple command computing multivariate group means and returning an object of class "matrix" rather than "list". Does any such function exist in standard packages? I'm beginning with R, so I'm sorry if the solution is trivial. Ondra Mikula

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

Hi, I wanted to use the predict.lm() function to compare the empirical data with the predicted values. The problem is that I have NAs in my data. I wanted to cbind my data.frame with the empirical values with the vector I get from predict.lm. But they don't have the same length because predict.lm just skip NA-predictions. Is there a way to get a vector with predicted values of the same

Hi, I'm trying to fit nonlinear mixed effects model using nlme function but getting an error message. Here is what I have: fitted_model = nlme(scores~spline(b1,b2,b3,kt,time), fixed = list(b1~1, b2~1, b3~1, kt~1), random = b1+b2+b3~1, groups= ~id, data = sdat, start = c(b1=3.5,b2=2,b3=.60,kt=3.5),verbose=T) Error: Error in

Hi, I am currently running the following negative binomial GLM: glm89.nb <- glm.nb(AvGUD ~ Year*Trt*Micro + (0 + Micro/Trt/Year)) where Year has 3 levels, Trt has 2 levels, and Micro has 3 levels. >From what I have read the above model has a 3 way interaction (Year*Trt*Micro), and Micro is nest within Trt and Trt is nested with Year (0 + Micro/Trt/Year). I was hoping someone could

Dear all, I am interested to calculate confidence interval for fitted values in general for non-linear regressions. Lets say we have y=f(x1,x2,..xN) where f() is a non-linear regression. I would like to calculate a confidence interval for new prediction f(a1,..,aN). I am aware of techniques for calculating confidence intervals for coeffiecients in specific non-linear regressions and with them

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.

...====================================================== 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: walmes@ufpr.br twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ========================================================================== [[alternative HTML version deleted]]

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

pairwise.t.test is returning NAs when one of the samples only has one entry, while TukeyHSD returns results (maybe not trustworthy or believable, but results). I stumbled on this because I did not realize one of my samples only had one entry while most of the others had several hundred, so I realize this is not a desirable situation. I'm really just curious about the difference between how

Do you know what steps I need to take to add a scale to a plot? I'm pulling my example out of "An Introduction to R: Software for StatisticalModelling & Computing" (see the R code around Figure 76). I would like to add a scale to the image produced by the following code.? I would like to the scale to list the distance?and units:? data(volcano) x <- 10*(1:nrow(volcano)) y