Displaying 13 results from an estimated 13 matches similar to: "How to plot multiple variables on the same graph"
2007 Oct 24
1
Error in nls model.frame
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,
2008 Jun 12
4
How to increase the for() loop speed?
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
2010 Nov 07
3
Computing ergodic mean with CODA
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
2012 Jan 03
6
calculate quantiles of a custom function
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
2012 Nov 02
2
Merge data frame with mispelling characters
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 =
2011 Jul 11
1
problem finding p-value for entropy in reldist package
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
2011 Dec 28
2
Contingency tables example
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
--
2012 Feb 28
1
Error in solve.default(res$hessian * n.used) :Lapack routine dgesv: system is exactly singular
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 <-
2011 Jul 08
4
How to label specific points on a scatterplot
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:
2011 Dec 28
2
Gale-Shapley Algorithm for R
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
2011 Jan 27
1
Saving log file in R and display run time
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
2010 Dec 23
1
Running sweave automatically using cygwin
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
2011 Jan 31
0
Function rearrange (quantreg)
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 <-