similar to: Repeated Measures/Linear Mixed Effects function

My data is the following: Time Resistance Temperature 5 2000 4 10 2200 8 15 2500 14 20 2900 20 25 3000 29 30 3100 38 35 3500 46 40 3800 47 45 3900 50 50 4000 51 I would like to create a scatter plot with Time on the x axis, Resistance on the y axis and Temperature on a second y axis. I have tried Plot(Time,Resistance,ylim=range(0,4500)

Quick question about the usage of glht. I'm working with a data set from an experiment where the response is bounded at 0 whose variance increases with the mean, and is continuous. A Gamma error distribution with a log link seemed like the logical choice, and so I've modeled it as such. However, when I use glht to look for differences between groups, I get significant

Dear R People: I have an AOV model that I get confidence intervals from via > confint(chick1.aov1) 2.5 % 97.5 % trtA 1.472085 1.607915 trtB 1.512085 1.647915 trtC 1.328751 1.464582 > I am using R2HTML to produce HTML output. However, the HTML code itself just has rounded values, i.e., 1.5 and 1.6. Has anyone run across this, please? Any suggestions would be much appreciated.

Dear R-help-list, I have a problem in which the explanatory variables are categorical, the response variable is a proportion, and experiment contains technical replicates (pseudoreplicates) as well as biological replicated. I am new to both generalized linear models and mixed- effects models and would greatly appreciate the advice of experienced analysts in this matter. I analyzed the

The Department of Statistics at the Faculty of Mathematics, Informatics and Statistics of the Ludwig-Maximilians Universit?t in Munich invites applications for a Professorship (W2) (6 years/tenure track) for Statistics with Special Emphasis on Computational Statistics as soon as possible. The professorship represents the entire field of computational statistics in research and teaching.

A non-text attachment was scrubbed... Name: not available Type: text Size: 3126 bytes Desc: not available Url : https://stat.ethz.ch/pipermail/r-help/attachments/19991130/5cb00c0f/attachment.pl

Hello, I have the following data.frame structure(list(Study = structure(c(1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 13L, 14L, 15L, 16L, 17L, 18L, 19L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 13L, 14L, 15L, 16L, 17L, 18L, 19L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 13L, 14L, 15L, 16L, 17L, 18L, 19L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 13L, 14L,

Hello, R-i-zens. I'm working on an data set with a factorial ANOVA that has a significant interaction. I'm interested in seeing whether the simple effects are different from 0, and I'm pondering how to do this. So, I have my.anova<-lm(response ~ trtA*trtB) The output for which gives me a number of coefficients and whether they are different from 0. However, I want the

Dear everyone, I would like to create a kite chart in which I plot densities (width of the vertical kites) in relation to sediment depth (on reversed y-axis) for 6 different locations (Distances from seep site, on x-axis on top of the plot). The dataset I would like to use is: Distance_from_seep_site Sedimentdepth Density 1100 0 107.8 1100 1 264.6 1100 2 284.2

Hi How would I analyse time series with - different lengths (i.e. one has 9 entries and the other has 14 entries) - different frequency (i.e. dates are random - no repeated length) - multiple values for the same time entry (e.g. 2009-10-23 below) i.e. my data takes the form: 1st time series 2009-10-07 0.009378 2009-10-19 0.014790 2009-10-23 -0.005946 2009-10-23 0.009096 2009-11-08 0.004189

Dear All, I get the following error when using either SASxport or Hmisc to read an .xpt file: > w <- read.xport("D:/consult/Trophos/dnp/base/TRO_ds_20081006.xpt") Erreur dans factor(x, f$value, f$label) : invalid labels; length 15 should be 1 or 14 > z<- sasxport.get("D:/consult/Trophos/dnp/base/TRO_ds_20081006.xpt") Erreur dans factor(x, f$value, f$label) :

Dear Users! I encountered with some problem in data reading while I challenged R (and me too) in a validation point of view. In this issue, I tried to utilize some reference datasets ( http://www.itl.nist.gov/div898/strd/index.html). And the result departed a bit from my expectations. This dataset dedicated to challenge cancellation and accumulation errors (case SmLs07), that's why this

This may be of interest to some of you. Janusz Kawczak. Associate Professor in Biostatistics or Statistics The Department of Mathematics is soliciting applications for the tenure track position of Associate Professorship in BioStatistics or Applied Statistics, beginning August 2003. Demonstrated prior experience in multidisciplinary medical or health related funded research is required. Expected

Hello, We have running lustre 1.8.1 and have met "No space lest on device" error when uploading 500 Gb small files (less then 100 Kb each). The problem seems to depends on the number of files. If we remove one file, we can create one new file, even with Gb size; but if we haven''t remove something we can''t create even very little file, as an example using touch

Hello all, I have a question concerning how to get the P-value for a explanatory variables based on GLM. I'll run multiple regressions with GLM, and I'll need the P-value for the same explanatory variable from these multiple GLM results. I check the help and there are quite a few Value options but I just can not find anyone about the p-value. Could anyone help me with that?

> Dear useRs, > > by a circular diagram representation I mean what you will get by entering > this at your R promt: > > pie(1:5) > > Nice to have R as a lingua franca :-) > > The folowing quote is from page 360 in this very interesting paper: > > @article{SpenceI2005, > title = {No Humble Pie: The Origins and Usage of a Statistical Chart}, >

Dear R list members, I'm encountering a problem with expand.model.frame(): Suppose that I define the following simple function (meant just to illustrate the problem): > fun <- function(model){ + expand.model.frame(model, all.vars(formula(model))) + } > and I have the following model, created with an explicit data argument: > mod Call:

Dear R-Users, Have a question about reading in some data and manipulating dates. I have a data set in excel which looks like this: Date PEG ETN HSP PTC 13/10/2004 41.92 64.75 29.86 9.27 14/10/2004 41.93 61.79 29.98 9.14 15/10/2004 41.69 62.7 30.09 9.04 18/10/2004

Hi, I am trying to plot the Case-Shiller index found at: http://www2.standardandpoors.com/spf/pdf/index/CSHomePrice_History_022445.xls The way I'm importing it into R is as follows: library(gdata) W <- read.xls("http://www2.standardandpoors.com/spf/pdf/index/CSHomePrice_History_022445.xls ", header=TRUE) attach(W) To give you and idea of what the data looks like: >

Dear R-users, I am in the process of creating new custom functions and am quite puzzled by some discrepancies in execution time when I run some R scripts that call those new functions. So here is the situation: - let's assume I have created two custom functions, called myg and myf; - myg is mostly a plotting function, which makes a heavy use of grid and lattice functions; - myf is a function