similar to: How to plot multiple variables on the same graph

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,

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

Hi all, I would like to compute ergodic mean using MCMC output from WinBUGS. I tried using CODA package, but it seems that it is not implemented yet. Could anyone help me to compute this? Attached to this email are my output and index files. Kind regards, Raquel -- Raquel Rangel de Meireles Guimar?es Doutoranda em Demografia raquel at cedeplar.ufmg.br

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 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 =

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

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

Dear all How can I obtain the data from the function "rearrange" in package quantreg More especifically, based on the example below (available in the help of the rearrange function), how can I access the data generated by "rearrange(zp)" ? data(engel) z <- rq(foodexp ~ income, tau = -1,data =engel) zp <-