similar to: nlme and spatially correlated errors

Hola a todos, Estoy intentando realizar un análisis de medidas repetidas en R y me gustaría pediros opinión porque estoy un poco atascada. Estoy utilizando las funciones lme y lmer. Tengo datos de censos de grillos (efi) realizados semanalmente (week) en tres parcelas (plot) diferentes, cada una de las cuales tiene seis transectos (transect) donde se realizan dichos censos. Los datos son de dos

I have a data frame with two columns, one of which is a factor (Species) and the other is numeric (BA, which stands for basal area). Here's a sample: Species BA ACSA 55.7632696 FRAM 122.9933524 ACSA 67.54424205 ACSA 89.22123136 ACSA 82.46680716 ACSA 22.46238747 ACSA 19.94911335 ACSA 20.42035225 ACSA 19.00663555 ACSA 21.67698931 ACSA 57.80530483 ACSA 30.31636911 Dead 43.98229715 Dead

I've been trying to figure out how to accomplish the following... I've got a list (returned from a function) and I would like to "cbind()" the lists together to create a "cross tab" report or simply bind them together somehow the function returns a list that looks like the following: > all$BM $species [1] "BM" $vbar.nobs [1] 3 $vbar.sum [1] 54.05435

http://procontrol.fi/test/ - Crashfix for multipart messages - More optimized mbox rewriting, should be "good enough" for most uses I think. I've been using it for a while now, seems to work. Still some problems: imap(cras): Apr 15 01:48:48 Warning: Our dotlock file /home/cras/mail//full-disclosure.lock was overridden And I'm sure it wasn't. Or at least shouldn't have

Hi everyone, I have been working some days in a nested loop in R but I can't find the solution. I have a data.frame with an unique ID for individuals and unique ID for different stands, for each indiviadual I have a dbh record and a SBA (stand basal area) field. Pma<-rep (1:40) P<-seq(1,4, 1) Plot<-rep(P,10) dbh2<-rnorm(40, mean=200, sd=5) SBA2<-rnorm(40, mean=10, sd=1) As

Hi Dovecot: 2.2.19.rc2 OS: CentOS Linux release 7.1.1503 (Core) FS: XFS Problem: Dovecot it does not include line breaks in the letter. Bug or feautre? Message size in fs: 1125 bytes Message: From cras at irccrew.org Tue Jul 23 19:39:23 2002 Received: with ECARTIS (v1.0.0; list dovecot); Tue, 23 Jul 2002 19:39:23 +0300 (EEST) Return-Path: <cras at irccrew.org> Delivered-To: dovecot

Hey everyone, I am new to R, and I'm making a scatter plot graph where i have a bunch of plots/points that fall into 9 unique categories. I want each category to have a unique colour, however, with the coding I have (below), the colour black is repeated for two of my plot types. Does anyone know a quick way to get 9 unique colours?? Coding: plotba = plot (predictedba ~ actualba,

Dear List, I'm an [R] novice starting analysis of an ecological dataset containing the basal areas of different tree species in a number of research plots. Example data follow: > Trees<-data.frame(SppID=as.factor(c(rep('QUEELL',2), rep('QUEALB',3), 'CORAME', 'ACENEG', 'TILAME')), BA=c(907.9, 1104.4, 113.0, 143.1, 452.3, 638.7, 791.7, 804.3),

Dear R-users, I want to fit a linear model with Y as response variable and X a categorical variable (with 4 categories), with the aim of comparing the basal category of X (category=1) with category 4. Unfortunately, there is another categorical variable with 2 categories which interact with x and I have to include it, so my model is s "reg3: Y=x*x3". Using fit.contrast to make the

Dear R-users. I'm struggeling to fix the superscript of a label of a figure axis. For some reason R doesn't recognize the "hat" symbol. plot(1,1,xlab="ligth intensity (PAR)",ylab=expression("mass Pteridium rhizomes" (gr/0.25m^2))) A very similiar scriptline does not give any problem at all: plot(1,1,xlab=expression("balsa plot basal

Hello R Help! I would like to make a legible boxplot of tree growth rates for each of seven tree species at each of seven different sites. It's a lot of data to put on one figure, I know. I made a beautiful, interpretable figure using color, but my target journal can't deal with color figures. I can use seven shades of grey to fill the boxes, but the figure then becomes uninterpretable -

Hi all, I am using the vegan package to run a prcincipal components analysis on forest structural variables (tree density, basal area, average height, regeneration density) in R. However, I could not find out how to extract factor loadings (correlations of each variable with each pca axis), as is straightforwar in princomp. Do anyone know how to do that? Moreover, do anyone knows

Hi all, I am totally new to new to wine, after installing it and configuring wine.conf I'm trying to get matlab works, but the program a second after starting cras, here's the output: wine matlab Warning: /usr/X11R6/bin/winereal not accessible from a DOS drive Warning: /usr/X11R6/bin/winereal not accessible from a DOS drive Warning: /usr/X11R6/bin/winereal not accessible from a DOS drive

L.S. I've some problems on a couple of my samba servers. The servers running CentOS release 5 with samba-3.0.25b-1.el5_1.4 The config is pretty standard, not much special in it. But the log files are filling up with error messages like this one: [2008/03/12 17:08:47, 0] smbd/service.c:set_current_service(184) chdir (/var/lib/samba/profiles) failed On the clients (all windows XP) are

Are there some known problems related to ctime changing by itself? I didn't find anything with google at least. tar seems to be complaining about it all the time: [cras@hurina] ~% tar cf test.tar /mnt/cygwin/home tar: Removing leading `/' from member names tar: /mnt/cygwin/home/Timo Sirainen/xxx/xxxxx.xxx: file changed as we read it This happens with my Debian/unstable (samba

Dear R users, I am looking for some advice on the proper construction of a mixed model in R, using the difference in means as the response and treating within-means residuals as a random effect. I have a dataframe (my own, a snippet of which is given below) that is composed of observations of pollen viability in flowers along tree branches. Flowers (1 to 3 per position) were collected from

Dear All, here a problem I think many of you can solve in few minutes. I have a dataframe which contains values of plot id, diameters, heigths and basal area of trees, thus columns names are: id | dbh | h | g head(ipso, n=10) id dbh h g 1 FPE0164 36 13.62 0.10178760 2 FPE0164 31 12.70 0.07547676 21 FPE1127 57 18.85 0.25517586 13 FPE1127 39 15.54 0.11945906 12 FPE1127

Dear All, I am not able to run the random forest with my dataset.. X<- 280 records with satellite data(28 columns) - B1min, b1max, b1std etc.. y<- 280 records with 3 columns - TotBasal Area, Stem density and Volume yref <- y[1:230,] #Keeping 1st 230 records as reference records want to set 0 to y values for records 231 to 280.. yimp <- y[231:280,] #records for which we want

Hi I am conducting 2-sample Kolmogorov Smirnov tests for my Masters project to determine if two independant tree populations have the same size-class distribution or not. The trees have been placed into size-class categories based on their basal diameters. Once I started running the stats on my data, I got confused with the results. Just to show an example of what I was testing I ran stats

I've been trying to figure out how to prevent a column that is the result of an aggregate function call so that I can use it in further calculations. For example, I would like to aggregate the expf for the data.frame by sp (character) and dbh (double d=rounded to integer) using the command: > st2 <- aggregate( ntrs$expf, by=list(sp=ntrs$sp,dbh=ntrs$dbh), sum ) > st2$expf <- st2$x