similar to: Calculate monthly means

Dear list, I've been learning how to make a 2x2 paneled dotplot in lattice without any previous experience using lattice. my code thusfar is: nut<-read.table("/Users/colinwahl/Desktop/nutsimp_noerror.csv", T, sep= ",") attach(nut) nut1<-data.frame(Nitrate, Total_Nitrogen, Phosphate, Total_Phosphorus) nut1<-as.matrix(nut1) rownames(nut1)<-group

Hi all, I've got some image plots on which I'd like to include some gene information (in the margins using mtext). Unfortunately, the description is rather long and will need to be wrapped to fit on several lines. From what I know about mtext, it's really only meant for single-line labels, not paragraphs. Here's some sample code of the idea I'm trying to accomplish:

Dear list users, I have trouble with covariance analysis. I measured nitrate concentrations in the soil (NO3) and the percentage of legumes (LEG, continuous), affected by 2 different CO2 concentrations (CO2, discrete). I suspect that CO2 has an effect on LEG and NO3, but also that LEG has an effect on NO3, so this is the formula I wrote to test this: NO3 ~ CO2 + LEG + CO2:LEG Will LEG be

I have a collection of datasets in separate data frames which have 3 independent test parameters (w, x, y) and one dependent variable (z) , together with some additional static test data on each row. What I want is a data frame which contains the test data, the parameters (w, x, y) and the mean value of all (z)s in the Z column. Each datasets has around 6000 rows and around 7 columns, which

Hi, how are you? My name is Irucka Embry and I would like assistance on properly subsetting some data obtained from a .csv file. > file <-"info_n.csv" > INFO<- getMetaDataFromFile(file) INFO is 92 observations of 14 variables > INFO<- subset(INFO,site.no==02169570) INFO is 0 observations of 14 variables When the site.no is used I should have 1 observation of 14

Hello, I'm trying to make a lattice plot (using xyplot()). I have included a "layout=c(3,4)" statement, giving me 12 plots per page and an "as.table=TRUE" statement, directing the way the plots are laid out. I have 18 plots altogether and so 6 of them end up on the second page. Everything looks fine for the first page, but the x-axis labels (e.g. 1993, 1994...) are all

I have a question concerning the Wilcoxon signed-rank test, and specifically, which R subroutine I should use for my particular dataset. There are three different commands in R (that I'm aware of) that calculate the Wilcoxon signed-rank test; wilcox.test, wilcox.exact, and wilcoxsign_test. When I run the three commands on the same dataset, I get different p-values. I'm hoping that

Hi, #make example data dat <- data.frame(matrix(rnorm(15),ncol=5)) colnames(dat) <- c("ab","cd","ef","gh","ij") If I want to get a subset of the data for the middle 3 columns, and I know the names of the start column and the end column, I can do this: mysub <- subset(dat,select=c(cd:gh)) If I wanted to do this just on the column names,

#is there a way to get NA in the table of descriptive statistics instead of the function stopping Thank you in advance #data x.f <- structure(list(Site = structure(c(9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L), .Label = c("BC", "HC", "RM119", "RM148", "RM179", "RM185",

Hello, I am trying to make a heatmap in R and am having some trouble. I am very new to the world of R, but have been told that what I am trying to do should be possible. I want to make a heat map that looks like a gene expression heatmap (see http://en.wikipedia.org/wiki/Heat_map). I have 43 samples and 900 genes (yes I know this will be a huge map). I also have copy numbers associated with

Hi, I am having trouble understanding how to approach a simulation: I have a sample of n=250 from a population of N=2,000 individuals, and I would like to use either permutation test or bootstrap to test whether this particular sample is significantly different from the values of any other random samples of the same population. I thought I needed to take random samples (but I am not sure how

Dear list, I have imported a spss data file in R, where a comma is used to separate the decimal numbers, e.g. 3,567 instead of 3.567 as in R. How can I replace the comma with a dot for all values in the data frame >kk a b c d e f 1 ,0925199910320613 82523 8 ,855 ,803 ,69 2 ,278314161923372 91657 26 1,285 1,032 ,823 3

Hi, all I didn't get any response from swig for my question. see if I can get some help here Thanks ---------- Forwarded message ---------- From: charlie <charlie.xia.fdu@gmail.com> Date: Tue, Nov 4, 2008 at 1:55 PM Subject: SWIG with R and C++ STL To: swig-user@lists.sourceforge.net Hi all, I am new to SWIG. I encountered some problem when I try to SWIG to R some C++ modules. Here

I measured nitrate concentration and primary production (PP) biweekly for 23 months in one headwater stream. I would like to use linear regression to determine if PP is related to nitrate concentration. My dataframe is called "data" and consists of the vectors Rdate, PP, and nitrate. Rdate is the observation date in class "date" and PP is primary production. I first

Dear group, i have 114 files (all have 5 columns) and variable number of rows. I want to read all the files and select the first column and 4th column and fill the values into a big matrix. I have a pre-made matrix that would have all rownames (1st column) and coloumn names (4th column from file to be read). mymat <- matrix(data = 1, nrow = length(gomfs), ncol = length(x))

Hello, my name is Toby, and i live in Sydney Australia. i would like to buy a 25Kg bag of strontium nitrate and have it delivered to an address in sydney. i understand that long distance postage can sometimes be very costly, so can you please give a a quote for the 25Kg bag, plus the postage cost to my address in sydney. thankyou for your time

Dear Rich, It depends how the data is generated. Although I am not an expert in ecology, I can explain it based on a biomedical example. Certain variables are generated geometrically (exponentially), e.g. MIC or Titer. MIC = Minimum Inhibitory Concentration for bacterial resistance Titer = dilution which still has an effect, e.g. serially diluting blood samples; Obviously, diluting the

Hello, I would like to perform an NMDS on the following: I have two independent variables, which are sites and treatments. I have 6 sites which are peatlands. I collected 5 replicates (at the same time) from each of the sites. I used each of the replicates in a treatment. There were 4 treatments. 2 were controls and the other two were a simulated 30 day drought and a simulated 60 day drought.

Dear Forum, I have a data frame as  mydat = data.frame(basel_asset_class = c(2, 8, 8 ,8), defa_frequency = c(0.15, 0.07, 0.03, 0.001)) > mydat   basel_asset_class defa_frequency 1                 2          0.150 2                 8          0.070 3                 8          0.030 4                 8          0.001 I need to get the subset of this data.frame where no of records for the

Dear Forum, I have a data frame as  mydat = data.frame(basel_asset_class = c(2, 8, 8 ,8), defa_frequency = c(0.15, 0.07, 0.03, 0.001)) > mydat   basel_asset_class defa_frequency 1                 2          0.150 2                 8          0.070 3                 8          0.030 4                 8          0.001 I need to get the subset of this data.frame where no of records for the