similar to: Logistic regression with factorial effect

Hi everybody, If I am correct, you can compare a model with random effect with the same model without the random effect by using the nlme function, like this: no.random.model <- gls(Richness ~ NAP * fExp, method = "REML", data = RIKZ) random.model <- lme(Richness ~NAP * fExp, data = RIKZ, random = ~1 | fBeach, method = "REML")

By any chance is it possible to "add" filled.contour(...) to an existing plot like you can do using image(..., add = TRUE). I looked at the filed.contour(...) and levelplot(...) documentation, but it didn't seem to be possible. I would like to add filled.plot(...) as shown below instead of using the image(...). Filled.contour seems to provide the smoothing interpolation

hello - a question about filled.contour plots, for which i haven't found a response in previous posts - sorry if already treated. i'd like to draw several filled.contour plots (that is, maps) on the same device (a postscript file, actually). I know about layout(matrix) , split.screen or par(mfrow) : it works well for simple plots, but with filled.contour plots, i get several pages

Hi everybody, I'm interested in evaluating the effect of a continuous variable on the mean and/or the variance of my response variable. I have built functions expliciting these and used the 'mle2' function to estimate the coefficients, as follows: func.1 <- function(m=62.9, c0=8.84, c1=-1.6) { s <- c0+c1*(x) -sum(dnorm(y, mean=m, sd=s,log=T)) } m1 <- mle2(func.1,

Hello, I need to expand the membership of a custom system user depending on the availability of some group on the target system i.e. user { "logger": name => "logger", ensure => "present", groups => ["adm", "wheel", "this _group_ if it exists"], shell => "/bin/bash"; } The

Dear sirs, I'm trying to build an Opus library for an embedded application to run on a Microchip PIC32MX processor and I get a lot of errors in compilation time. My question is, what should be the minimum files and options needed to have such building, having in mind that only the very basic functionality in fixed point arithmetic is needed? Looking forward to your answer. Kind

I've tried this on both Solaris 8 and SuSE Enterprise 9 (64-bit). I get a assert-crash when using a gzipped mbox folder . OK Logged in. . SELECT test.gz * FLAGS (\Answered \Flagged \Deleted \Seen \Draft $Label4 $Label2 $Label1 $Label3 $Label5) * OK [PERMANENTFLAGS ()] Read-only mailbox. * 167 EXISTS * 0 RECENT * OK [UNSEEN 2] First unseen. * OK [UIDVALIDITY 1096038620] UIDs valid * OK

Package: xen-hypervisor-4.0-amd64 Version: 4.0.1-1 Hardware specs: Dell PowerEdge R410 32GB RAM DDR3 1TB RAID-1 HD (Perc H200 RAID controller) Processor: 2 - Intel XeonE5620 2.4Ghz 12M Cache Turbo HT 1066MHz SO: Debian Squeeze last update: 14/10/2010 uname -a: "Linux spo-r1-u4 2.6.32-5-xen-amd64 #1 SMP Fri Sep 17 22:00:48 UTC 2010 x86_64 GNU/Linux" At boot, xen-hypervisor is loaded

Hello, Below is some example code that should reproduce an error I'm encountering while trying to create a tiff plot with two subplots. If I run just the following bit of code through the R GUI the result is what I'd like to have appear in the saved tiff image: x<-seq(0:20) y<-c(1,1,2,2,3,4,5,4,3,6,7,1,1,2,2,3,4,5,4,3,6) plot(x,y,type="l",las=1,ylim=c(0,12))

I am trying to model the observed leaching of wood preservative chemicals from treated wood during an outdoor experiment where leaching is caused by rainfall events. For each rainfall event, the amount of rainfall was recorded as well as the amount of preservative chemical leached. A number of climatic variables were measured, but the most important is the amount of rainfall. I have tried a

Hi, I am trying to format some data (example matrix "m" below) for which each data point has 2 associated values separated by a comma. I want to replace values <3 with "0" to give the example output below. I have got as far as: out<-lapply(strsplit(m,","),as.numeric) Failed to identify anything along the lines of "[out<3]<-0" that works... Thank

Hello list i have an issue with my dahdi_channels.conf in span 1 when i use it like below i can do my outband calls without issue ; Span 1: TE2/0/1 "T2XXP (PCI) Card 0 Span 1" (MASTER) group=0,11 context=from-pstn switchtype = euroisdn signalling = pri_cpe channel => 17-31 context = default group = 63 but when i add in channel 1-15 like: channel => 1-15,17-31 i receive all

hello, i try to plot 3 graphs which have the same x.axis underneath each other. i'd like the plots to be aligned without margings between the boxes and draw a single x axis beneath the lowest plot. i managed to get the alignment by setting par(mar), but the middle box is stretched and i cant't figure out how to get around this. par(pin) was my guess, but this doesn't do the job

Hello all, I am an undergraduate student who is having syntax issues trying to get error bars on my graph. This is the data, which I assigned the name "Saline" to. Time Average SEM 1 -20 0.000000 0.0000000 2 3 30 0.000000 0.0000000 4 45 3.227902 0.7462524 5 60 5.066664 1.1623944 6 80 6.107491 1.5027762 7 110

Hey Everyone, I'm pretty new at R and wanted to try and make some graphs with dummy data before using it to analyze my own. I successfully made a bar graph and error bars, but I can't figure out how to align them properly (currently they are not centered on the bars and some of them aren't even close). Here's the code that I'm using: > marks <- sample(4:10, size=50,

I define the following function to convert a t-value with degrees of freedom DF to another t-value with different degrees of freedom fullDF: tConvert <- function(tval, DF, fullDF) ifelse(DF>=1, qt(pt(tval, DF), fullDF), 0) It works as expected with the following case: > tConvert(c(2,3), c(10,12), 12) [1] 1.961905 3.000000 However, it gives me warning for the example below although the

hi I'm trying to use the function odiag(x) for matrix calculations. I would like to assign new values to an off-diagonal in a matrix. When I use the diag (x) function I could write something like p<-matrix(seq(1:12),ncol=4) p.new<-matrix(rep(0,12),ncol=4) diag(p.new)<-diag(p) p.new But this won't work with odiag. How can I turn odiag (x) into something like diag (x) in order

Dear friends, I have the following dataframe: YEAR <- c(1996 , 1997, 1998, 1999, 2000, 2001, 2002, 2003, 2004 ) T_MAR <- c(2.8, 6.5, 5.4,2.4, 4, 4.1, 3, 4.4, 4.5) T_APR <- c(5.7, 7.8, 7.7, 4.6, 4.7, 6.2,5.7, 5.9, 7) T_MAY <- c(7, 8.8, 10, 6, 5.5, 7.6, 8.5, 7.3, 10.2) BUD <- c(87, 98, 93, 85, 89, 91, 87, 92, 92) BUD_SE <- c(3.6, 2, 2.4, 4, 2.4, 2.4, 4, 2.4, 3) g1 <-

Dear R-users, I'd like to plot a three-dimensional surface and at the meantime I'm using an array. I would like to have the values of my first matrix in the heights of the plot and the colors of the single facet taking into account the second matrix. I hope that the next code will help all of you to understand better my issue, Thanks in advance, Giancarlo ############################ ##

Hi, I used simple linear regression with the R software and EXCEL on the same data. Although , I find the same R2=0.84, I find different estimated values (intercept and slope). For the R software (slope =0.0009, Intercept = -0.1478), for EXCEL (slope =927.7, Intercept = 154,41). When I use the estimated values from the R software, the results seem bad, however the results of Exel seem correct.