Re: [R] Comparison of demographics between 2 study samples

ur goal is to demonstrate that the samples come from the same population then you probably should take a look at equivalence testing. Andrew -- Andrew Robinson Department of Mathematics and StatisticsTel: +61-3-8344-6410 University of Melbourne, VIC 3010 Austr

Re: [R] instal tar.gz package on windows

For cross-building a mature package, I stumbled across a convenient utility provided by the R Windows maintainer: _http://win-builder.r-project.org/ _This takes a little while, so for iteration it is much better to install the Rtools, but for a package or two it is great. Best, Andrew Harris

[R] New FITS file utility package

apart the headers and data in the file. Many thanks to the R core development team for developing and supporting such an excellent program for data exploration and analysis. Andrew Harris Universisty of Maryland __ R-help@r-project.org mailing list

Re: [R] Decomposing tests of interaction terms in mixed-effects models

On Mon, Aug 04, 2008 at 02:51:48PM +0200, Peter Dalgaard wrote: > Andrew Robinson wrote: > >On Mon, Aug 04, 2008 at 10:17:38AM +0200, Peter Dalgaard wrote: > > > >>Andrew Robinson wrote: > >> > > > >That is a neat idea, thanks, Peter, but

Re: [R] Decomposing tests of interaction terms in mixed-effects models

On Mon, Aug 04, 2008 at 10:17:38AM +0200, Peter Dalgaard wrote: > Andrew Robinson wrote: > >Dear R colleagues, > > > >a friend and I are trying to develop a modest workflow for the problem > >of decomposing tests of higher-order terms into interpretable sets of > &g

[R] Decomposing tests of interaction terms in mixed-effects models

ample) lme(Y ~ A + AC, random = ~ 1 | Block, data = example) ########## Are we doing anything obviously wrong? Is there another approach to the goal that we are trying to achieve? Many thanks, Andrew -- Andrew Robinson

[R] Changing values

Hello-- I am a relatively new user to R and I cannot find the information I need. Please help. I have a very large data set with values including letters, numbers, and symbols (sometimes within the same vector value [ie X9-]. I've imported the data using read.fwp and it arrives in list fo

Re: [R] Rf_error crashes entire program.

Switching all of my printfs to Rprintf fixed the problem the errors now proceed correctly. I was unable to reproduce the error, but now that is irrelevant, at least to me. Thanks for the help. -Andrew [[alternative HTML version deleted]] __ R

[R] Rf_error crashes entire program.

break; } and the output from running in batch mode that I get is this: Da: 0.18034.988e-017 PING and here the program crashes. I've tried this in multiple places and sometimes the error is thrown sometimes not. Does anyone have an idea of what is going on. I could not find any di

Re: [R] with lapply() how can you retrieve the name of the object

Thanks everyone for your suggestions (and sorry for the delay in the acknowledgement). Jorge and Jim, thanks for pointing out your approach. Andrew On Fri, Jul 18, 2008 at 7:02 PM, Jorge Ivan Velez <[EMAIL PROTECTED]>wrote: > > Dear Andrew, > > Following Jim Holtman

[R] help with integrate function

inite ones? Thanks very much Andrew Brown __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-cont

[R] with lapply() how can you retrieve the name of the object

out the names of the thing that is being lapplied? Thanks, Andrew [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org

[R] Help with error in "if then" statement

r than TRUE/FALSE? Am I using "if" incorrectly? In which case would it be more appropriate to perhaps create subsets of the data points based on < or > ycrit? Thanks in advance for any guidance-- Andy -- Andrew J. Rominger Department of Biological Sciences Stanford University [

[R] Non-normal data issues in PhD software engineering experiment

* Must have one observation per cell. * ERROR * Completion of computation impossible. 5) I have looked briefly at the non-parametric approaches to regression - there seems to be many (http://socserv.mcmaster.ca/jfox/Courses/Oxford-2005/R-nonparametric-regression.html) paths that can be taken. I need

Re: [R] problem with NA and if

Hi Keld you should read ?sum. sum(c(1,2,NA), na.rm=TRUE) Cheers Andrew On Fri, Jul 04, 2008 at 08:29:34AM +0200, Keld J?rn Simonsen wrote: > Hi > > I would like to sum a number of time series, some of them having NA's > > Standard action is here that if I sum a value

Re: [R] Problems with lm()

12 1682 23 121 11 1542 24 122 21 1247 25 131 22 1235 26 132 12 1605 27 141 11 1598 28 142 21 1718 Andrew On Thu, Jun 19, 2008 at 04:29:16PM +0800, leeznar wrote:

Re: [R] keeping original order in factor()

Thanks for everyone's suggestions. I think factor(foo, levels = unique(foo)) works best for my needs. (By the way, I'm still trying to figure out how to use the ordered option in factor().) Andrew On Wed, Jun 18, 2008 at 12:32 AM, Peter Alspach <[EMAIL PROTECTED]> wrote: >

[R] keeping original order in factor()

; "a" #not "a" "b" I thought having ordered=T would do the trick. Thanks, Andrew __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the p

Re: [R] in axis() suppressing axis line but keeping tick marks

Thanks for tracking this down. Andrew On Mon, Jun 16, 2008 at 4:49 PM, Peter Dalgaard <[EMAIL PROTECTED]> wrote: > Andrew Yee wrote: > >> I've been trying to figure out a parameter that will let you separately >> adjust the parameters for the axis line from the tick

[R] in axis() suppressing axis line but keeping tick marks

I've been trying to figure out a parameter that will let you separately adjust the parameters for the axis line from the tick mark. In the following example, I would like to suppress the axis line, but keep the tick marks. Thanks, Andrew foo <- data.frame(x=1:3, y=4:6) plot(foo$x, foo

Re: [R] "False convergence" in LME

These are great tips from Spencer. The other thing that I have found useful is to use a different optimizing algorithm. You can do this by: control=lmeControl(opt = "optim") Good luck! Andrew On Sat, Jun 14, 2008 at 09:45:22AM -0700, Spencer Graves wrote: > This is a c

Re: [R] controlling location of labels in axis()

Thanks for the reply. I think I've figured it out, you can set this with the mgp parameter. So I'd use the following statement instead: axis(1, at=foo$plot.x, labels=foo$plot.x, mgp=c(3,0.5,1)) #this brings the axis labels closer to the axis line Andrew On Thu, Jun 12, 2008 at 9:5

[R] controlling location of labels in axis()

s(1, at=foo$plot.x, labels=foo$plot.x) I'd like to be able the control the y location of the labels (i.e. move up or down in the graph). Thanks, Andrew [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.eth

[R] R loop

set<-as.Date(c("2008-01-01"), format = "%Y-%m-%d") set for(i in a$date){ if (a$date[i]>set & a$date[i]set +7){a$wk<-2}else {a$wk<-3} } a Kind regards andy Andrew McFadden MVS BVSc Incursion Investigator Investigation & Diagnostic Centres - Wallaceville B

Re: [R] Bimodal Distribution

nge the cutoff and the weights to suit your need. It's best to double-check by plotting an empirical density of the numbers generated. I hope that this helps, Andrew On Thu, May 29, 2008 at 04:05:29PM -0600, Mike Williams wrote: > Hello R Users, > > I am doing a Latin Hypercube t

Re: [R] can I do this with R?

ture of the dataset. As I recall, subset selection was preferred for finding a small number of large effects, lasso (L1) for finding a small to moderate number of moderate-sized effects, and ridge (L2) for many small effects. Can you provide any references to more up-to-date simulations that you

Re: [R] How to remove NAs and lme function

Jen, try na.action = na.exclude Andrew On Wed, May 28, 2008 9:26 pm, Jen_mp3 wrote: > > I am working on a project to find a model for the concentration of > dissolved > oxygen in the river clyde. Ive fitted a linear mixed model as > lme(DOW~Temperature+Salinity+Year+factor(

Re: [R] Need help for nlme

extraneous. Andrew. On Sun, May 25, 2008 at 01:36:49PM -0500, [EMAIL PROTECTED] wrote: > Hi everyone, > I try to write a module based on nlme however R always shows me the > error message > Error in eval(expr, envir, enclos) : object "y" not found. Does anyone &g

Re: [R] Some problems with Sweave

Martin, try omitting the results=tex argument. Andrew On Fri, May 23, 2008 at 10:16:33AM +0200, [EMAIL PROTECTED] wrote: > Dear R users, > I'm working in a brief R-tutorial to a group of students. To make that I'm > using Sweave but I've got two problems: > > Fir

Re: [R] Cumulative average

Jacques, you should be able to construct a solution from cumsum(). Cheers Andrew On Wed, May 21, 2008 at 08:48:29PM -0500, Jacques Wagnor wrote: > Dear List, > > Does there exist a function that calculates a cumulative average? > Neither running() from library(gregmisc) nor

[R] Legend in scatterplot and adding regression lines

panel=panel.superpose ,xlab="Time (weeks)",ylab="Farrowing rate", key = list(points = Rows(trellis.par.get("superpose.symbol"), 1:4), text = list(levels(Dataset$Operator)), columns = 4)) Regards Andy Andrew McFadden MVS BVSc Incursio

Re: [R] lme nesting/interaction advice

On Mon, May 12, 2008 at 10:50:03AM +0100, Federico Calboli wrote: > > On 12 May 2008, at 01:05, Andrew Robinson wrote: > > >On Mon, May 12, 2008 at 10:34:40AM +1200, Rolf Turner wrote: > >> > >>On 12/05/2008, at 9:45 AM, Andrew Robinson wrote: > >> >

Re: [R] lme nesting/interaction advice

On Mon, May 12, 2008 at 10:34:40AM +1200, Rolf Turner wrote: > > On 12/05/2008, at 9:45 AM, Andrew Robinson wrote: > > >On Sun, May 11, 2008 at 07:52:50PM +0100, Federico Calboli wrote: > >> > >>The main point of my question is, having a 3 way anova (or anc

Re: [R] lme nesting/interaction advice

y with the request at the bottom of each R-help email, and provide commented, minimal, self-contained, reproducible code, that actually ran, ideally with fewer value judgements, you might get more attention from the people who are smarter than you and me, but have less time than either o

[R] Finding unmatched data between two dataframe using several factors

, "zane"),z=c(10,20,30,40,50)) dat2 I would like to know data that from dat2 that doesn't appear in dat1 based on both the x and y factors ie (A5, dan, 40) and (A6,zane, 50). I have tried two approaches but have not been successful ? nomatch <- subset(dat1, is.element(?

Re: [R] lme() with two random effects

helpful. Cheers Andrew On Fri, May 09, 2008 at 04:09:07PM -0300, Mike Lawrence wrote: > Hi all, > > I have collected response time data from 178 participants ('sub') for > each combination of 4 within-Ss factors ('con','int','tone','cu

[R] p values for polychor

ration of the two tails of the distribution? > x <- as.integer(runif(20) > 0.5) > y <- as.integer(runif(20) > 0.5) > p <- tetrachoric.test(x, y) > p$statistic [1] -0.2616866 > p$p.value [1] 0.7935631 > p$var [,1] [1,] 0.1

[R] Spatial join between two datasets using x and y co-ordinates

dog And therefore limit the data in dataset 1 to that which has the same x and y co-ordinates as dataset 2. Any suggestions? Regards Andy Andrew McFadden MVS BVSc Incursion Investigator Investigation & Diagnostic Centres - Wallaceville Biosecurity New Zealand Ministry of Agriculture and

Re: [R] improvement of Ancova analysis

, and if so, why, b) produce a model to estimate the effect of introduction on weight, with the other variables being nuisance variables, and if so, why, c) something else (and why) because all these factors affect the position that you could adopt about the questions that you have. Andrew On Sun

Re: [R] GLMM and data manipulation (2nd try)

On Sat, May 03, 2008 at 06:43:22AM +1000, Andrew Robinson wrote: > On Fri, May 02, 2008 at 03:06:31PM -0500, Giovanni Petris wrote: > > > > Hello, > > > > I posted a question yesterday but I got no replies, so I'll try to > > reformulate it in a more conci

Re: [R] GLMM and data manipulation (2nd try)

time = rep(c(1,2), 1600), subject=rep(1:1600, each=2)) test.lmer <- lmer(response ~ time + (1|subject), data=ratings, family=binomial) but I don't know if you think that's clean or not. Andrew > Thank you in advance, > Giovanni Petris >

Re: [R] Errors using nlme's gls with autocorrelation

You're running out of RAM. Your options are 1) run code on a machine with more RAM 2) try the model on fewer observations 3) try a simpler model. Andrew On Fri, May 02, 2008 at 12:37:15PM -0700, zerfetzen wrote: > > Hi, > I am trying out a generalized least squares method

Re: [R] alternatives to RColorBrewer?

Great, thanks, that was helpful. Andrew On Thu, Apr 24, 2008 at 2:15 PM, Achim Zeileis <[EMAIL PROTECTED]> wrote: > On Thu, 24 Apr 2008, Andrew Yee wrote: > > > I've found RColorBrewer useful for its qualitative palettes, but wished > that > > it could generate m

[R] alternatives to RColorBrewer?

is doesn't generate enough distinct colors when the number of palettes is large). Thanks, Andrew [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the pos

Re: [R] using rbind() on multiple objects at once

Thanks! On 4/18/08, Katharine Mullen <[EMAIL PROTECTED]> wrote: > do.call(rbind, list.foo) > > On Fri, 18 Apr 2008, Andrew Yee wrote: > > > Is there an efficient way to use rbind() with the five dataframes > described > > in the following example: > > &

[R] using rbind() on multiple objects at once

st.foo[[5]]) #is there an easier method? For example, I naively thought you could do something like rbind(list.foo[[1:5]]) #gives an error message but that results in an error message. Thanks, Andrew [[alternative HTML version deleted]] __

[R] Instantaneous hazard of an event i.e. death

"hazard.ratio.plot" but lack the experience to know how to code this Anyone be able to help? Regards Andrew Andrew McFadden MVS BVSc Incursion Investigator Investigation & Diagnostic Centres - Wallaceville Biosecurity New Zealand Ministry of Agriculture and Forestry Phone 04 89

Re: [R] why incomplete lmer output?

Hi Maria, the output is not incomplete. See the FAQ here: http://cran.r-project.org/doc/FAQ/R-FAQ.html#Why-are-p_002dvalues-not-displayed-when-using-lmer_0028_0029_003f and the lengthy discussions about p-values and lmer (try Google). Andrew On Tue, Apr 15, 2008 at 10:51:46AM +0200, Maria

Re: [R] why is text being rasterized with text()

Thank you for looking at my code. You are absolutely right. I'm mortified that I didn't see the text() tucked inside the for loop. Once I took it out of the loop, the text() works fine. Thanks, Andrew On Thu, Apr 3, 2008 at 1:52 AM, Prof Brian Ripley <[EMAIL PROTECTED]> wro

[R] why is text being rasterized with text()

cter(seq(from=4,to=10)),">11") frame() par(usr=c(0,10,0,10) ) for (r in 1:9) { rect(r,3,r+1,3.75, border = NA, col=heat.colors(9)[r]) #this text does not appear to be rasterized text(r+0.5,2.75, value.seq[r], cex=0.5) #this text appears to be rasterized text(5.5,2.25,"express

Re: [R] peak finding

hich(c(smoothed.dx,NA) > 0 & c(NA, smoothed.dx) < 0) My experience is that this approach sometimes requires some fine-tuning, eg in fitting the smooth. I hope that this helps Andrew On Mon, Mar 24, 2008 at 11:23:57PM -0400, Research Scholar wrote: > Hi > Thanks for replyin

Re: [R] Great difference for piecewise linear function between R and SAS

th major. I think you should try Simon Wood's "Generalized Additive Models: An introduction in R." 2006 Chapman and Hall. I reviewed it recently and I think it's really very good. Cheers Andrew -- Andrew Robinson Department of Mathematics and StatisticsTel: +

[R] Updating selected variables in a data frame

ps, that I can update the main dataset with the new values from the recodes dataset? Thanks very much for your advice! Regards, Andrew C. Ward CAPE Centre Department of Chemical Engineering The University of Queensland Brisbane Qld 4072 Australia __ R-help

[R] Distances between two datasets of x and y co-ordinates

x <- c(2660156,2663703,2658165,2659303,2661531,2660914) y <- c(6476767,6475013,6475487,6479659,6477004,6476388) data2<-cbind(x,y) x <- c(266500,261) y <- c(6478767,6485013) data1<-cbind(x,y) Any suggestions on how to do this would be appreciated. Regards Andrew Phone 04 894

[R] R code for selecting places spatially and by time

,2,2,2,1,1) ) date<-as.character(rep(as.Date(c("2008-01-02","2008-01-17","2008-01-01") ,format = "%Y-%m-%d"),c(4,4,2))) cbind(x,y,date) Regards Andrew Andrew McFadden MVS BVSc Incursion Investigator Investigation & Diagnostic Centres - Wallaceville Biosecur

Re: [R] [OT] "normal" (as in "Guassian")

of "off topic", I wonder if Stigler's > >> Law is self-referential? That is, should Stigler's Law more correctly > >> be attributed to someone else? > >> > > > > No. If Stigler's Law were named after some prior person, &g

Re: [R] discrete variable

Pete, try x <- c(0, 3, 4, 3, 15, 5, 6, 5) table(x) or length(table(x)) Cheers Andrew On Sun, Mar 02, 2008 at 09:27:20PM +0100, Pete Dorothy wrote: > Hello, > > I am sorry for asking such a basic question. I could not find an answer to > it using google. > > I have a

Re: [R] A more idiomatic way to write this

How about x <- x / rep(divs, rep(1000, 1000)) ? Cheers, Andrew On Sun, Feb 24, 2008 at 10:36:23PM -0300, Andre Nathan wrote: > Hello, > > I have a vector of 1,000,000 numbers and another vector of 1,000 > divisors. What I'd like to do is to divide the first 1,000 nu

Re: [R] skip non-converging nls() in a list

Hi Marc, use try() Cheers Andrew On Mon, Feb 18, 2008 at 01:45:35PM -0500, Marc Belisle wrote: > Howdee, > > My question appears at #6 below: > > 1. I want to model the growth of each of a large number of individuals using > a 4-parameter logistic growth curve. > > 2

Re: [R] Creating a data.frame

Hi Joe, cbind coerces the data to be the same type. Much nicer is: d <- data.frame(x=x, y=y, z=z) Cheers Andrew On Wed, Feb 13, 2008 at 05:17:32PM -0500, Joe Trubisz wrote: > OK...newbie question here. > Either I'm reading the docs wrong, or I'm totally confused. >

[R] If else in R code

t;) c<-data.frame(a,b) I would like to write an if statement where if "a" is not null return "a", if "a" is null (NA) return "b" less two days. Could someone help?. Andrew Thi

[R] [OT] - standard errors for parameter estimates under ridge regression and lasso?

Dear R community, I'm curious to know how people go about estimating standard errors for parameter estimates after model selection by ridge regression and the lasso. Do you have any practical or theoretical advice? Warmly, Andrew -- Andrew Robinson Department of Mathematics and Stati

[R] changing color scheme in R in UNIX

Hi, I've been using R in Windows but am now starting to use it more often in the UNIX environment. In Windows, I'm used to the text provided by the user appearing in red, and the output from R appearing in blue. Is there a way to achieve this in UNIX? Thanks, Andrew [[altern

Re: [R] aaMI

Hi Navish, did you run require(aaMI) ? Cheers Andrew On Thu, Jan 17, 2008 at 02:17:12AM -0800, navish wrote: > > hi > i am new to R language. I want to use aaMI package which calculates the > amino acid mutual interaction for a given protein sequence. I had installed > the p

Re: [R] Grouping data

Hi Kimmo, try cut() to create a factor with levels according to the range of values, and (among other options) table() to make the table. Cheers Andrew. On Wed, Jan 16, 2008 at 10:06:23PM +0200, K. Elo wrote: > Hi, > > I am quite new to R (but like it very much!), so please apo

[R] [Off Topic] searching for a quote

Dear community, I'm trying to track down a quote, but can't recall the source or the exact structure - not very helpful, I know - something along the lines that: 80% of [applied] statistics is linear regression ... ? Does this ring a bell for anyone? Thanks, Andrew -- Andre

Re: [R] Omnibus main effects in summary.lme?

)) should avoid the error and give a conservative estimate of the significance of your interaction. see also: http://tolstoy.newcastle.edu.au/R/e2/help/06/10/3565.html and related posts. A ----- Dr. Andrew Beckerman

Re: [R] Rscript on OSX

ies on OSX. R is not designed to run from Terminal.app (actually, terminal does not do graphics well for R). X11 provides the behaviour that you would be used to in Linux etc. BTW - there is a R-Mac-help group as well for macintosh specific issues - [EMAIL PROTECTED] And

[R] survreg and cluster and frailty

od ratio test=176 on NaN df, p=NaN n= 92 Warning message: In survpenal.fit(X, Y, weights, offset, init = init, controlvals = control, : Inner loop failed to coverge for iterations 2 Many thanks, Andrew FIRST 20 Lines of data > ss[1:20,] start stop event tree hag entcir entarea depth n

Re: [R] groupedData function not found

Hi Andrea, did you try require(nlme) ? Adding a package is not necessarily the same as loading it. Cheers Andrew On Thu, Dec 27, 2007 at 10:42:27PM +, andrea previtali wrote: > > Hello, > I'm trying to use the groupedData function and R is giving me the message: > Er

Re: [R] Rbind-ing a list into one item

Hi Arjun, try do.call() Cheers Andrew On Wed, Dec 26, 2007 at 04:23:51PM -0500, Kondamani, Arjun (GMI - NY Corporate Bonds) wrote: > Hi, > > I am doing the following: > > 1. I have a list of files.. Files1=list.files("some > directory",pattern="some patt

Re: [R] Dual Core vs Quad Core

ut generally things that take a single thread will remain on a single core. As for RAM, if you're doing memory-bound work you should certainly be using a 64-bit machine and OS so you can utilize the larger memory space. ------

Re: [R] way to check if the evaluation of an expression returns an error?

?try Cheers, Andrew On Fri, Dec 14, 2007 at 02:10:52PM -0800, Thomas Pujol wrote: > Is there a recommended or "good" way to check if the evaluation of an > expression returns an error? > > e.g. > var(NA) > I wish var(NA) to return NA or "err"

[R] Calculating Rsquared values in rpart

Hi there In rpart, one can get a graph of R-squared (using rsq.rpart (fit)), in which the x axis is the number of splits, and which contains two lines - an "apparent" R squared and an Rsquared based on the x error. I would like to caclulate these R-squared values, but cannot work out from the

[R] Permutation of a distance matrix

Hi there, I would like to find a more efficient way of permuting the rows and columns of a symmetrical matrix that represents ecological or actual distances between objects in space. The permutation is of the type used in a Mantel test. Specifically, the permutation has to accomplish somethi

[R] adding in missing values in a sequence

3 3 3 5.6 4 4 6.7 5 5 NA 6 6 7.1 7 7 12.5 8 8 14.5 9 9 16.8 1010 3.4 I'm running R2.6.0 on Mac OSX. Thanks in advance, Andrew Hoskins PhD Candidate Deakin University, Australia Email: [EMAIL PROTECTED] [[alternative HTML version

Re: [R] adding in missing values in a sequence

That seems to have solved it. Thanks Mark. On 15/11/2007, at 9:57 AM, Marc Schwartz wrote: > > On Thu, 2007-11-15 at 09:41 +1100, Andrew Hoskins wrote: >> Hi, >> >> I have a data frame with two columns of data, one an indexing column >> and the other a data column

[R] cor.test() and binary sequences

with is two binary sequences. just series of 0 and 1. will the confidence intervals and p-values generated by cor.test() still apply? cheers, andrew. -- Get a free email account with anti spam protection. http://www.bluebottle.com

Re: [R] interpreting an LME regression result...

ht be other explanations for the change. I hope that this helps, Andrew On Sat, Nov 10, 2007 at 12:58:21AM -0700, Johan Jackson wrote: > Any help would be most appreciated. (Don't make me get down on my > hands and knees and beg for help, cause I'll do it!!) My boss has me > lea

[R] (no subject)

n coefficient (0.5198252). an alternative question would be: given two sequences and a calculated correlation coefficient, with what probability could i assert that the underlying processes are indeed correlated and that the calculated correlation coefficient does not simply arise by chance

[R] Dringende Aufmerksamkeit BITTE LESEN

MR. ANDREW L. J. BARENDS Standard Bank Plc. Johannesburg, South Africa1 ID Number: 0092-042. UNTERSTÜTZUNG. Dringend GESCHAEFTSANGEBOT Mein Name ist Andrew L. J. Barends und ich bin der Leiter des Corporate Affairs Committee in der Standard Bank of South-Africa PLC in Südafrika. Ich

Re: [R] Help in error of mixed models

paces in your code? Good luck Andrew On Sun, Nov 04, 2007 at 09:18:09AM -0200, Bernardo Rangel Tura wrote: > > Hi R-masters > > I read the article: Bivariate analysis of sensitivity and specificity > produces informative summary measures in diagnostic reviews. > > In this paper

Re: [R] perl module for R

http://www.omegahat.org/RSPerl/ -- Andrew J Perrin - andrew_perrin (at) unc.edu - http://perrin.socsci.unc.edu Associate Professor of Sociology; Book Review Editor, _Social Forces_ University of North Carolina - CB#3210, Chapel

Re: [R] Calculate aggregate differences

rvation that is calculated as the spurious difference between the first measurement of the group and the last measurement of the previous group. If that doesn't help, let me suggest that you construct a small worked example that will show us the kind of input and output that you are thinking about.

Re: [R] Calculate aggregate differences

Hi Cristian, instead of aggregate, how about something like: n <- dim(my_agg)[1] my_agg$vol.diff <- my_agg$Vol - c(NA, my_agg$Vol[1:(n-1)] my_agg <- my.agg[my.agg$Age > min(my.agg$Age),] (assumes same minimum age for all treatments) (not checked) Cheers, Andrew On Thu, Nov 01

[R] Deparsing part of a list argument

how I can use that to get just the functions that aren't named. Any ideas? Thanks, Andrew __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.

Re: [R] Repeated Measures/Linear Mixed Effects function

Hi Keith, it seems like a good starting position. I recommend that you spend some time studying Pinheiro and Bates's book to see where t ogo from here. Cheers Andrew On Mon, Oct 22, 2007 at 02:58:51PM -0800, Keith Cox wrote: > I have three columns of data, Xc, Trt and fish. Th

Re: [R] confusion with R syntax

Hi Mark, what's happening here is that R is applying the one-dimensional subscripting operations sequentially. Try x <- seq(1,10) x[2:4][1] Cheers Andrew On Thu, Oct 11, 2007 at 03:34:22PM -0400, Leeds, Mark (IED) wrote: > I just noticed something by accident with R syntax that

[R] multcomp and lme4

Dear list members, Can anyone please point to an example of how to use glht(multcomp) with lmer objects? I am trying: summary(glht(lmerObject, linfct = mcp(x = "Tukey"))) as I would for a glm object, but with no luck. Thank you, Andy. [[alternative HTML version deleted]] __

[R] SETAR and nonlinear model estimation

ist or are included in the estiamtion of the model. Any help provided will be greatly appreciated Andrew Phiri Post-gratuate student at the North West University (NWU), South Africa - Luggage? GPS? Comic books? [[alternative HTML versi

[R] Origin label in lattice

just below and to the left of the lower left corner of the panel. I've unsuccessfully tried to do this using the xscale.components argument, and I'm having trouble figuring out which other arguments or functions might be helpful. I would greatly appreciate any advice! Thanks, Andrew

<    3   4   5   6   7   8