No No wrote:
> Does this help?
> > pdf(file="lag.pdf")
> > lag.plot(a)
> > lag.plot(b)
> > dev.off()
> After that you can open each page of the "lag.pdf" file with GIMP for
> further manipulation. It gives each plot on a different page, but no
> plot is replaced.
That would be a last resort,
Ok I guess it's time to ask.
So I want to plot my data. It's my data from a frequency table, "temp". My
formula is just a Gaussian eq. I have done the nls function to get my
parameters and now I want to do the whole plot (...) and then lines(..)
This is what I have done.
> temp
bin x
1 -4
There are older versions of gridBase suitable for 2.2.1 at CRAN, e.g.
src/contrib/Archive/gridBase_0.4-2.tar.gz. Why not just install that?
On Tue, 13 Jun 2006, Jonathan Dushoff wrote:
>
> I am running R 2.2.1 on a University-supported linux installation based
> on Redhat EL3. I am sorry that i
Jonathan Dushoff wrote:
> I am running R 2.2.1 on a University-supported linux installation based
> on Redhat EL3. I am sorry that it did not occur to me to mention this
> before; I updated R very recently, with the most recent version
> available for EL3 at http://cran.cnr.berkeley.edu/bin/linux/
I am running R 2.2.1 on a University-supported linux installation based
on Redhat EL3. I am sorry that it did not occur to me to mention this
before; I updated R very recently, with the most recent version
available for EL3 at http://cran.cnr.berkeley.edu/bin/linux/redhat/el3/.
Looking at the gr
Does this help?
> pdf(file="lag.pdf")
> lag.plot(a)
> lag.plot(b)
> dev.off()
After that you can open each page of the "lag.pdf" file with GIMP for
further manipulation. It gives each plot on a different page, but no plot is
replaced.
2006/6/13, Gad Abraham <[EMAIL PROTECTED]>:
>
> Hi,
>
> I'm try
Does this help?
> pdf(file="lag.pdf")
> lag.plot(a)
> lag.plot(b)
> dev.off()
After that you can open each page of the "lag.pdf" file with GIMP for
further manipulation. It gives each plot on a different page, but no plot is
replaced.
Hope it helps
2006/6/13, Gad Abraham <[EMAIL PROTECTED]>:
>
>
Hi,
I'm trying to plot several lag.plots on a page, however the second plot
replaces the first one (although it only takes up the upper half as it
should):
par(mfrow=c(2,1))
a<-sin(1:100)
b<-cos(1:100)
lag.plot(a)
lag.plot(b)
What's the trick to this?
I'm using R 2.2.1 (2005-12-20 r36812) on
Dear all,
Is any body know the R packages for processing remote sensing images
(not the Rimages) like ERDAS Imagine does? Cheers
Yong
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting gui
Trying to follow the instructions on the page
http://tolstoy.newcastle.edu.au/R/help/05/12/17153.html (to make a
figure with an inset), I typed
> install.packages("gridBase", repos="http://cran.r-project.org";)
and was rejected.
no package 'gridBase' at the repositories in: download.pack
Thanks Jim. This seems to work for my case already.
Eric
On 6/13/06, Jim Lemon <[EMAIL PROTECTED]> wrote:
> Eric Hu wrote:
> > Hi I am trying to plot two data set in the same picture window without
> > overlapping with each other. I am using the format plot(x1,y1,x2,y2)
> > but get the following
Eric Hu wrote:
> Hi I am trying to plot two data set in the same picture window without
> overlapping with each other. I am using the format plot(x1,y1,x2,y2)
> but get the following error message:
>
>
>>plot(as.numeric(r0[,4]),as.numeric(r0[,7]),as.numeric(r0[,4]),as.numeric(r0[,7][ind[,1]]))
>
Hi I am trying to plot two data set in the same picture window without
overlapping with each other. I am using the format plot(x1,y1,x2,y2)
but get the following error message:
>plot(as.numeric(r0[,4]),as.numeric(r0[,7]),as.numeric(r0[,4]),as.numeric(r0[,7][ind[,1]]))
Error in plot.window(xlim, yl
Dear list members,
I need a dat set with follow features:
-Binomial response in longitudinal form
-overdispersion
thanks
Leidy Rocío
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide!
Thanks! Both responders understood what I was after despite my poor
explanation and came up with very helpful responses. If anyone else has
an idea, please share!
David Kling
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/list
Hi List,
Does any of you know about an implementation of the Local Ripley's K and
L functions for R? The function has been described in "Second-Order
Neigborhood Analysis of Mapped Point Patterns (Getis & Franklin, 1987,
Ecology 68:3)
Thanks.
D.
-
S
I am guessing you are starting Rcmdr using
> library(Rcmdr).
If this is the case, repeating the library command does nothing, since
Rcmdr package is already attached. Try
> Commander()
instead. Make sure you use capital C.
Cheers,
Andy
__
Andy Jaworski
518-1
Assuming that yr and conf are the two factors referred to in the
description, create a function f which calculates the ith row
of the output and use sapply like this:
attach(data)
f <- function(i) {
hsgpa <- na.omit(hsgpa[-i][conf[-i] == conf[i] & yr[-i] == yr[i]])
if (length(hsgpa
Not exactly sure what you mean, but here is something that might be close.
I used only a subset of your data to see it this is what you want. This
computes the mean of all hpgpa, excluding that row:
> data[x[['2005.e']],] # subset of your data for yr=2005, conf='e'
casehsgpa yr conf
7
It could be a little bug or a fault in my setup but I notice that after I
end an Rcmdr session I can't restart it without restarting Rgui (Windows)
Do I need to add something to my environment variables? Is it an issue with
Tcl/Tk or it's settings? Or what?
It is not a big deal but if someone has a
Hi,
I want to plot non parametric estimates of the empirical hazard function for
right censored data. I've tried many functions from different packages
(muhaz, Design, survival, eha, event), but none of them gave me what I
wanted. Am I missing something?
Here's what I want. The data below is the
Hello:
I hope none of you will mind helping a newbie. I'm a student research
assistant working with a large data set in which observations are
categorized according to two factors. I'm trying to calculate the group
mean and variance of a variable (called 'hsgpa' in the example data
presented
Hello,
is there any function/package in R for multivariate quantile estimation?
(I couldn't find any for the moment)
Thank you!
Anna
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! h
Before I reinvent the wheel, has anyone set up scripts to read in
array CGH data files from NimbleGen, for use in the aCGH package or
a similar package?
I have normalized NimbleGen files
([array_id]_[identifier]_normalized.txt in the NimbleGen naming scheme)
several to a directory, to read in an
Hi,
I have some problems I have a Cauchy distribution with density function
f(x) = sigma / (pi * (sigma^2 + (x- miu)^2) ),
where sigma = scale and miu = location (in my case sigma = 3, miu = 0), and I
have to find with bootstrap
E | sigma_estimated^3 sigma^3 | (#),
In using Nelder-Mead outside R, I find it critical to restart the
algorithm (repeatedly) after it thinks it's found a solution, to see if
it can do better. I can't say whether the R and Matlab implementations
do this automatically or not.
on 6/12/2006 3:00 PM Anthony Bishara said the followin
Thanks for the feedback. I should've mentioned before that the
function is non-smooth. Also, it has a 3-element free parameter vector, and
I've been using a grid of 27 vectors of starting parameters.
Anthony
-Original Message-
From: Prof Brian Ripley [mailto:[EMAIL PROTECTED]
Try this:
score[tapply(rownames(score), score$id, tail, 1),]
On 6/12/06, Tony Chu <[EMAIL PROTECTED]> wrote:
> Dear R users:
>
> I have a small test dataframe as the follows :
>
> math = c(80,75,70,65,65,70)
> reading = c(65,70,88,NA,90,NA)
> id = c('001','001','001','002','003','003')
> score
one approach is the following:
math <- c(80, 75, 70, 65, 65, 70)
reading <- c(65, 70, 88, NA, 90, NA)
id <- c('001', '001', '001', '002', '003', '003')
score <- data.frame(id, reading, math)
#
ind <- cumsum(tapply(score$id, score$id, length))
score[ind, ]
I hope it helps.
Best,
Dimi
Have you also tried 'nlminb'? A year or so ago, Doug Bates switched
from optim to nlminb for mixed-effects estimation. I'm not certain I
know why, but I think 'nlminb' may automatically adjust the 'scale'
parameter by default, while "optim" does not automatically adjust the
comparab
Unless you know the function to be non-smooth, I suggest you use
method="BFGS" in R.
BTW, all such algorithms are only designed to find local minima, and so
the choice of starting point may be crucial.
On Mon, 12 Jun 2006, Anthony Bishara wrote:
> Hi,
> I'm having a problem converting a Matlab
Dear R users:
I have a small test dataframe as the follows :
math = c(80,75,70,65,65,70)
reading = c(65,70,88,NA,90,NA)
id = c('001','001','001','002','003','003')
score = data.frame(id, reading, math)
> score
id reading math
1 001 65 80
2 001 70 75
3 001 88 70
4 002
Hi,
I'm having a problem converting a Matlab program into R. The R code works
almost all the time, but about 4% of the time R's optim function gets stuck
on a local minimum whereas matlab's fminsearch function does not (or at
least fminsearch finds a better minimum than optim). My understanding i
Version 0.9-99 of package "subselect" is now on CRAN.
The subselect package has functions that search for k-variable subsets
of a p-variable (p>k) data set that are optimal under some
criterion. Search algorithms include a full search algorithm "leaps",
based on Furnival and Wilson's leaps and b
Using Perl or R shouldn't be an either/or issue.
We can call Perl code from R (via RSPerl at least)
and then we get the benefit of well tested, fast code
that might exist in Perl to read the log files
and the ability to do interactive, exploratory analysis in R.
And reading the file in Perl does g
"Gabriel Diaz" <[EMAIL PROTECTED]> writes:
> I'm taking an overview to the project documentation, and seems the
> database is the way to go to handle log files of GB order (normally
> between 2 and 4 GB each 15 day dump).
> In this document http://cran.r-project.org/doc/manuals/R-data.html,
> say
Dear Eric
Do you really have habitats nested within lagoons or are they
partially crossed (meaning that you have the same habitats in
different lagoons)?
If you have them perfectly nested, I think that you cannot
calculate both a fixed effect for habitats and a random effect
for lagoon (see the e
On Sat, 10 Jun 2006, Gregory Pierce wrote:
> Hello friends and fellow R users,
>
> I have successfully tabulated and entered my survival data into R and
> have generated survival curves. But I would like to be able to determine
> what the survival rates are now at one month, three months, six mont
I'm working with lm for some time, so I know the function.
The problem was that I purchased the book to learn more about linear models
and couldn't find the equivalent for the maximum likelihood in R, I used lm,
mle, and a few others, but it never allow you to set a variance parameter.
But after re
I wouldn't use a DBMS at all -- it is not necessary and I don't see
what you would get in return. Instead I would split very large log
files into a number of pieces so that each piece fits in memory (see
below for an example), then process them in a loop. See the list and
the documentation if you h
Thanks.
Reading the discussion of differences between S and R, I had almost come to
this conclusion but wanted to be sure. Since R keeps its working data in
RAM, the use of chapters seemed unlikely, but I wanted to be sure.
I am, however, surprised this has not come up before. The online
document
Dear All,
The issue I reported earlier seems to be restricted at least to Mac
OSX GUI using 2.3.1. It does not occur in the Terminal in Mac OSX (R
version 2.3.1) or in Linux (R version 2.2.1). I will separately send
this to r-sig-mac.
Best,
On 12 Jun 2006, at 11:26,
Hello
Thanks all for the answers.
I'm taking an overview to the project documentation, and seems the
database is the way to go to handle log files of GB order (normally
between 2 and 4 GB each 15 day dump).
In this document http://cran.r-project.org/doc/manuals/R-data.html,
says R will load all
On Mon, 2006-06-12 at 18:18 +0530, Arun Kumar Saha wrote:
> Dear all R users,
>
> I am wondering whether there is any way to replace all "NA" values in a data
> frame by some numerical value, suppose 1000?
>
> Thanks and Regards
Hi, Arun,
sapply(dat, function(x) {x[is.na(x)] <- 1000; x})
e.g.:
On 6/12/2006 8:48 AM, Arun Kumar Saha wrote:
> Dear all R users,
>
> I am wondering whether there is any way to replace all "NA" values in a data
> frame by some numerical value, suppose 1000?
In a vector it is easy, e.g. x[is.na(x)] <- 1000. A dataframe is a list
of vectors, so you could itera
Hi
yes, but I do not recommend you to do it
> df<-data.frame(a=1:10, b=11:20)
> df[5,2]<-NA
> df
a b
1 1 11
2 2 12
3 3 13
4 4 14
5 5 NA
6 6 16
7 7 17
8 8 18
9 9 19
10 10 20
> df[is.na(df)]<-1000
> df
ab
1 1 11
2 2 12
3 3 13
4 4 14
5 5 1000
Hello. I am working with a matrix 612x618 and when I transform this data
in a geodata, and when, for example, I want to calculate the mean with the
funtion summary, the program say that the vector is too big and don't
calculate. I want know if I do anything wrong or the program, with a geo
On Mon, 12 Jun 2006, Stefano Calza wrote:
> Hi all.
>
> I may missing something here, but if I do summary.default(1:), I get:
>
> Min. 1st Qu. MedianMean 3rd Qu.Max.
> 12500500050007500
>
>
> but if I do summary.default(1:10001) I get:
>
> Min. 1st Qu.
Dear all R users,
I am wondering whether there is any way to replace all "NA" values in a data
frame by some numerical value, suppose 1000?
Thanks and Regards
[[alternative HTML version deleted]]
__
R-help@stat.math.ethz.ch mailing list
https:
Try
summary(1:10001, digits = 5)
or
old <- options(digits = 8)
summary(1:10001)
Note default for digits= in:
args(summary.default)
On 6/12/06, Stefano Calza <[EMAIL PROTECTED]> wrote:
> Hi all.
>
> I may missing something here, but if I do summary.default(1:), I get:
>
> Min. 1st Qu.
Hi all.
I may missing something here, but if I do summary.default(1:), I get:
Min. 1st Qu. MedianMean 3rd Qu.Max.
12500500050007500
but if I do summary.default(1:10001) I get:
Min. 1st Qu. MedianMean 3rd Qu.Max.
1250150
On Mon, 12 Jun 2006, Jeff Miller wrote:
> Cameron and Trivedi in their 1998 Regression Analysis of Count Data refer to
> NB1 and NB2
>
> NB1 is the negative binomial model with variance = mu + (alpha * mu^1)
> yielding (1+alpha)*mu
>
> NB2 sets the power to 2; hence, variance = mu + (alpha*mu^2)
>
ZhanWu Dai wrote:
> I am an initial user of R. Could you give me some explanations or examples on
> how to solve the first order differential equations by the first-order
> Runge-Kutta method?
>
> Thank you very much
>
> Kind regards
>
>
>
> -
On Mon, 12 Jun 2006, Lorenzo Isella wrote:
> Dear All,
>
> A simple question: packages like fitdistr should be ideal to analyze
> samples of data taken from a univariate distribution, but what if
> rather than the raw data of the observations you are given directly
> and only a histogram?
(It is
Dear All,
A simple question: packages like fitdistr should be ideal to analyze
samples of data taken from a univariate distribution, but what if
rather than the raw data of the observations you are given directly
and only a histogram?
I was thinking about generating artificially a set of data
corr
The fourth-order Runge-Kutta method is in the odesolve package on CRAN,
but they manual says, only for didactial purposes. First order
differential equations are generally conveniently solved by the
lsoda function in this package. I hope this helps, even if it does not
directly answer your questi
Dear all,
I am having a problem using grid when rotating a viewport. It seems
to plot everything on a grey background colour which I am not able to
get rid of. Even book examples such as that that plot figure 5.10 in
P. Murrell's R Graphics book show the same behaviour.
-BEGIN PGP SIGNED MESSAGE-
Hash: SHA1
Allen S. Rout wrote:
>
>
> Don't expect a warm welcome. This community is like all open-source
> communities, sharply focused on its' own concerns and expertise. And,
> in an unusual experience for computer types, our core competencies
> hold little
Hi
I would recommend to use list instead
how do you know that the result shall not be
AB
1 NA
3 5
4 6
5 3
2 5
> A<-1:5
> B<-4:7
> L<-list(A,B)
> L
[[1]]
[1] 1 2 3 4 5
[[2]]
[1] 4 5 6 7
If you insist on data frame you has to tell the program which cells
are to be empty or
I am an initial user of R. Could you give me some explanations or examples on
how to solve the first order differential equations by the first-order
Runge-Kutta method?
Thank you very much
Kind regards
__
R-help@stat.math.ethz.ch mailing list
https
On Mon, 2006-06-12 at 14:44 +0530, Arun Kumar Saha wrote:
> Dear all r-users,
>
> Suppose I have two data frame:
>
> A
> 1
> 3
> 4
> 5
> 2
>
> and
>
> B
> 5
> 6
> 3
> 5
>
> Now I want combine this two data frames without losing any value from either
> data frame. More precisely I want to see
>
Dear all r-users,
Suppose I have two data frame:
A
1
3
4
5
2
and
B
5
6
3
5
Now I want combine this two data frames without losing any value from either
data frame. More precisely I want to see
AB
1 5
3 6
4 3
5 5
2 NA
I tried with cbind function but failed, as it only works whe
> De : [EMAIL PROTECTED]
> [mailto:[EMAIL PROTECTED] De la part de Wouter
> Envoyé : vendredi 9 juin 2006 19:15
> Hi list,
>
> I'm looking for a way to calculate quantiles (such as in
> quantile()), but with the ability to assign a different
> weight to each member of the sample vector.
> [...
Hi
see ?rect
e.g.
> image(1:12, 1:12, outer(1:12, 1:12,"+"))
> for (i in 1:12) rect(10+(i/12-1/12), 8, 10+i/12,
10,col=heat.colors(12)[i])
>
HTH
Petr
On 9 Jun 2006 at 17:19, Yves Magliulo wrote:
From: Yves Magliulo <[EMAIL PROTECTED]>
To: r-help@stat.math
On Sat, 10 Jun 2006, Prof Brian Ripley wrote:
> ?regex does describe this:
>
> A range of characters may be specified by giving the first and last
> characters, separated by a hyphen. (Character ranges are
> interpreted in the collation order of the current locale.)
>
> You did not te
Hi
not sure what you want but if you just would like to filter numeric
columns
file[,sapply(file, is.numeric)]
HTH
Petr
On 10 Jun 2006 at 22:26, Milton Cezar wrote:
Date sent: Sat, 10 Jun 2006 22:26:35 -0300 (ART)
From: Milton Cezar <[EMAIL PROTECTED]>
To:
66 matches
Mail list logo