similar to: Computing ergodic mean with CODA

Hi all, I have two basic questions, hope you should help me: 1. How do I save a log file in R with the results? For example, in Stata it can be done by using "log using c:\...\test.txt" 2. How do I display the execution time of one function like the one below: CalculaCorrelacao <- function(construto, n) { library(polycor) for (i in 2:n){ for (j in i:n+1){ x1

Hi all, Hope someone could help me. I am trying to run automatically the conversion of an Rwn file to a tex file. I am using windows 7, and cygwin. I tried to run automatically the Sweave.sh script, in its the most recent version available at R webpage: http://cran.r-project.org/contrib/extra/scripts/Sweave.sh Unfortunately, I got this error message: =========================== Raquel at

How to plot multiple variables on the same graph Dear R users, I want to plot the following variables (a, b, c) on the same graph. The x-axis must be the variable mat and the graph must have the type = "l". How can I do that?? > a [1] 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 > b [1] 0.90665296 0.82473871 0.75269217 0.68917606 0.63304639 0.58332339 [7]

May I have a question on how to solve the following problem by R code? Mainly we want to solve the equation show in the attached image. The equation is a continuous version of Markov process. In the equation, we have been able to achieve two things using R code: [1] From year-2009 sample data, we can estimate the marginal density ?f(x ; 2009)? by using R function ?density()? [2] From

Dear R users, I would like to know if there is a way to increase the for() loop speed because in my routine the calculations are too slow. Best regards. Rafael Barros de Rezende Cedeplar - Center for Development and Regional Planning Face, UFMG ([1]http://www.cedeplar.ufmg.br) -- Esta mensagem foi verificada pelo sistema de antiv??rus e acredita-se estar livre de

Dear all, In many papers regarding time series analysis of acquired data, the authors analyze 'marginal distribution' (i.e. marginal with respect to time) of their data by for example checking 'cdf heavy tail' hypothesis. For i.i.d data this is ok, but what if samples are correlated, nonstationary etc.? Are there limit theorems which for example allow us to claim that

I''ve run across something that would save me days of trouble. Situation, the contents of one ZFS file system needs to be moved to another ZFS file system. The destination can be the same Zpool, even a brand new ZFS file system. A command to move the data from one ZFS file system to another, WITHOUT COPYING, would be nice. At present, the data is almost 1TB. Ideally a "zmv" or

Is there any implicit function in R to minimize a function with many variables? Many thanks Raquel ---------------------- R Granell, Rheumatology Unit & School of Mathematics University of Bristol Raquel.Granell at bristol.ac.uk Tel:07968079410

Error in model.frame When I run the following nls model an error message appears and I dont know how to solve that. Could you help me?? mat = c(1,2,3,4,5,6,7,8,9,12,16,24,36,48,60) for (i in 1:length(j30)) { bliss = nls(c(j[i,1:length(mat)]) ~ b0 + b1*((1-exp(-k1*mat))/(k1*mat)) + b2*(((1-exp(-k2*mat))/(k2*mat))-exp(-k2*mat)), start = list(k1=0.1993, k2=0.1993, b0= 22.0046,

Hi, I'm running a tutorial ("Meta-analyses of data from two (or more) microarray data sets"), which use wgcna package. I have an error in the function modulePreservation (it is below). I'm using R2.13 Can you help me? Do you know, what is happens? Thanks Raquel multiExpr = list(A = list(data=t(badea)),B = list(data=t(mayo))) # two independent datasets (dim = 13447 x 36) mp =

Hello, I want to use optim with simulated data. My program generates a sample and then calls optim to optimise a function that depends on the generated sample. I want to do this with 100 samples. But optim does not 'work' for the 100 samples, sometimes I get a message saying that sigma matrix was not invertible and therefore the optimisation was not carried out, and then my program

Hi, I guess that my problem has an obvious answer, but I have not been able to find it. Suppose I create a custom function, consisting of two beta-distributions: myfunction <- function(x) { dbeta(x,2,6) + dbeta(x,6,2) } How can I calculate the quantiles of myfunction? I have not seen any continous function treated in the docs, and applying the "quantile function" gives me an

Hello, <http://tolstoy.newcastle.edu.au/R/help/04/07/0312.html#0313qlink1> Could anybody tell me how to generate discrete Laplacian distribution? I need to sample uma discretised Laplacian density like this: J( g -> g´) ~ exp (-lambda | g´ - g |) g in {0,…, gmax} Thanks, Nicolette [[alternative HTML version deleted]]

Hello dear R-helpers, I'm working with R-2.15.2 on Windows 7 OS. I'm stucked with a merge of two data frames by characters. In each data frame I got two different list of names, that is my main-key to be merged. To figure out what I'm saying, I build up a modified "?merge" example, with errors by purpose: # Data for authors: authors <- data.frame( surname =

Hola, gracias por la respuesta, No me funcionó debido a que los nombres de las variables no están seriadas, es decir, los nombres de las variables son del tipo: x23 x25 x30, x301 x320, x80. Entonces me da este error: Error in eval(predvars, data, env) : object 'pot24' not found. Debido a que pot24 no existe, ya que de pot23 se brinca a pot30. En Stata es algo muy simple de hacer, solo

Hi, I am using the reldist package and having problems determining the p-value for the entropy value from the reldist function. I am able to properly determine the entropy value, but cannot figure out what function to use to find the p-value. I have tried using rpy, rpluy (which provides p-values for the polarization values) and investing the results from reldist(). Thus, far I cannot find the

I have this question: http://r.789695.n4.nabble.com/file/n4240145/data%5D.jpg This table and data represents a sample of social housing residents in Copenhagen and I'm asked to analyse the data with a view to seeing what might be related to residents satisfaction levels. And I just don't know where to start, if anyone could help this would be greatly appreciated Regards Sam --

Hi there! I´m a noob when it comes to R and I´m using it to run statisc analysis. With the code for ARIMA below I´m getting this error: Error in solve.default(res$hessian * n.used) :Lapack routine dgesv: system is exactly singular The code is: > s.ts <- ts(x[,7], start = 2004, fre=12) > get.best.arima <- function (x.ts, maxord=c(1,1,1,1,1,1)) + { + best.aic <- 1e8 + n <-

Command I am using for the plot: plot(Raw[][Plate==101]~well[][Plate==101], xlab="Well", ylab="Raw", main="Plate 101") I only want to label points on the graph where well equals B13, G13, K14 and N14 with the name of the well. Thank you for your help. -- View this message in context:

Dear R-helpers, I'm not a speciallist in writing complex functions, and the function still very rusty (any kind of suggestions are very welcome). I want to implement Gale-Shapley algorithm for R Language. It is based on http://www.jstor.org/stable/10.2307/2312726 Gale and Shapley (1962) , and it has evolved to