similar to: Crash with seasonal ARIMA

Hi, is there a way I can calculate a summary statistics for a columns matrix let say we have this matrix x <- matrix(sample(1:8000),nrow=100) colnames(x)<- paste("Col",1:ncol(x),sep="") if I used summary summary(x) i get the output for each column but I need the output to be in matrix with rownames and all the columns beside it this how I want it

Dear all, I am pretty new to R only having an introduction course, so please bare with me. I am doing my PhD at The Max Planck Institute of Immunobiology where I am analyzing some calorimetry data from some mice. I have a spreadsheet consisting of measurements of the respiratory exchange rate at different time points measured every 9 minutes over some days. My goal is "bin"/average the

Hi, I want to do a meta-analysis with count data for treatement/control cases. Mi problem is that I need to use zero values (an informative zero value) for the mean and standard deviation for one of the treatement, but R has a problem: "Studies with zero values for sd.e or sd.c get no weight in meta-analysis". I can agroup the case by Family (byvar=Family). ¿Can you help me? Thanks!

Hi! I want to do a meta-analysis with count data for treatement/control cases. Mi problem is that I need to use zero values (an informative value) for the mean and standard deviation for one of the treatement, but R has a problem: "Studies with zero values for sd.e or sd.c get no weight in meta-analysis". I can agroup the case by Family (byvar=Family). ¿anybody help me? Thanks! >

My goal is simple: calcuate GC content of each sequence in a list of nucleotide sequences. I have figured out how to vectorize, but all my attempts at memoization failed. Can you show me how to properly memoize my function? There is a StackOverflow post on the subject of memoization, but it does not help me: http://stackoverflow.com/questions/7262485/options-for-caching-memoization-hashing-in-r

This error keeps appearing when i perform a paired t-test in R Error in t.test.default(payoff, paired = T) : 'y' is missing for paired test This is the method i have used > read.table("MeanPayoff.txt",header=T) Open Closed1 47.50000 42.37502 49.25000 50.00003 50.00000 49.80004 33.50000 20.00005 34.75000 33.88006 35.50000 20.50007 33.35000 12.87508 50.00000

Hello, I'm trying to do bar plot where 'sex' will be the category axis and 'occupation' will represent the bars and the clusters will represent the mean 'income'. sex occupation income 1 female j 12 2 male b 34 3 male j 22 4 female j 54 5 male b 33 6

Dear R-helper: Suppose we have a matrix: Gene sample1 sample2 Gcnt1 12.0000 52.80000 Max 8.8000 39.10000 Tmem176b 67.9000 304.70000 Shmt2 8.6000 42.40000 Rtn4 11.5000 57.70000 Il17re 7.6000 38.80000 Bclp2 6.2000 32.10000 Mobkl3 4.4000 32.20000 Akr1b10 3.4000 30.10000

Hi R-users: [R 2.2 on OSX 10.4.3] I have a (sparse) vegetation data frame with 500 rows (sampling units) and 177 columns (plant species) where the data represent % cover. I need to summarize the cover data by returning the names of the most dominant and the second most dominant species per plot. I reduced the data frame to omit cover below 5%; this is what it looks like stacked. I have

I have many text files in the format below and in certain rare instances such as below there can be nothing in one of the fields so a double comma is written but I won't know this because I am reading in many,many files sequentially. # TEXT FILE 2004-02-10 00:01:31.00000,,105.60000000 2004-02-10 00:01:32.00001,,105.60000000 2004-02-10 00:01:45.00000,,105.60000000 2004-02-10

Hi I am following the model building strategy that is outlined in the Pinheiro and Bates book wrt including covariates but am having a problem with the plot. Basically I am using 4 covariates (1 of them is continuous) and 3 of them are fine but the 4th one is being shown as a scatterplot despite the fact that it is a factor. I have explicitly declared this to be a factor (pcat<-as.factor(pcat))

I have a time series for two variables measured simultaneously, but at different intervals. The variables are not independent, so the patterns in the times series are very similar (one variable goes up when the other goes up, etc). For example, let's pretend my data look something like this (but with more noise)... var1_times<-c(0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10) var1_values<-c(0,

I have a .txt file obtained by saving a data frame in which the first four columns are factors (but represented as 1,2,3 etc). The first four lines are "Pig" "Evit" "Cu" "Litter" "Start" "Weight" "Feed" "Time" "4601" "1" "1" "1" 26.5 26.5 NA 1 "4601" "1"

Hi there, This seems like it should be simple. I have a data frame of climate data sampled every 10 min. I want to average the entire data frame into 30 min values (i.e., one value for each half hour). Functions like running.mean give me a moving average but I want to reduce the size of the entire frame.. Any ideas how? Cheers! Example of my data timestamp Voltage LwTempDownelling

Hi, I'm a bit confused how julian() works. If I understand right, it returns the number of days since the origin. I have a vector: > SLDATX[1:10] [1] "1986-01-06" "1986-01-17" "1986-02-02" "1986-02-04" [5] "1986-02-04" "1986-02-21" "1986-03-06" "1986-03-25" [9] "1986-04-06"

Dear UseRs,Extremely sorry for a basic question. I have a matrix of 19 rows and 365 columns. what i want to do is the following...First i want to leave out column number 1 and want to calculate the row wise mean of the remaining columns, which will obviously give me 365 values in one column, and then subtracting these values from the column i left out i.e. col=1 then i want to leave out column 2

Hi Everyone, I have a dataset like this mean sd 0% 25% 50% 75% 100% n BODY TEMPERATURE 36.41099 0.4015699 35.1 36.22222 36.5 36.66667 37.1 89 DIASTOLIC BLOOD PRESSURE 73.60079 9.4656186 50.0 67.00000 73.0 80.00000 95.0 253 HEIGHT 171.94000 9.2011670 153.5 166.50000 173.0 176.25000 190.0

I have data that looks like below and I would like to re-order the values within the "date" column. I would like to have each year organized like so: Sep-71 Oct-71 Nov-71 Dec-71 Jan-71 Feb-71 Mar-71 Apr-71 May-71 Jun-71 Jul-71 Aug-71 Sep-72 Oct-72 etc... Is there any way I can order the column in my own fashion and just move Sep-Dec to the beginning of each year? I am planning on doing

Hi 2 all ! I have problems with concurrent filesystem actions on a ocfs2 filesystem which is mounted by 2 nodes. OS=RH5ES and OCFS2=1.2.6 F.e.: If I have a LV called testlv which is mounted on /mnt on both servers and I do a "dd if=/dev/zero of=/mnt/test.a bs=1024 count=1000000" on server 1 and do at the same time a du -hs /mnt/test.a it takes about 5 seconds for du -hs to execute: 270M