Say I have:
> set.seed( 1 )
> m <- matrix( runif(5^2), nrow=5, dimnames = list( c("A","B","C","D","E"),
> c("O","P","Q","R","S") ) )
> m
O P Q R S
A 0.2655087 0.89838968 0.2059746 0.4976992 0.9347052
B 0.3721239 0.94467527 0.1765568 0.7176185 0.2121425
C
hong shen yahoo.com> writes:
>
> I encountered a situation that a data frame is defined by two packages. Both
of them are loaded by library().
> 2. If I want to reference the data frame from package A insted of B, how can I
do it?
Either change the loading sequence of library(). Or, if you wan
Is it possible to interpolate a value for x with knowledge of y?
For example, approx(x, y, xout) will give me y's given a set of x's,
which is opposite to what I'm after. I've tried switching x and y,
e.g., approx(y, x, xout), but in a real data set it is possible to
have more than one y for a gi
Dear readers,
I would like to install a local binary package (.tar file) on a Mac (OSX
10.5.6) with R 2.8.0...
I tried :
mypkgdir = "/Users/philippesaner/Desktop/Rpkgs"
install.packages("R2WinBUGS_2.1-12.tar", destdir = mypkgdir, lib =
"/Library/Frameworks/R.framework/Resources/library", repos=
try this:
set.seed(1)
m <- matrix(runif(5^2), nrow = 5,
dimnames = list(c("A","B","C","D","E"), c("O","P","Q","R","S")))
v <- c(m)
names(v) <- paste(rownames(m), colnames(m)[col(m)], sep = ".")
# or
# names(v) <- outer(rownames(m), colnames(m), paste, sep = ".")
v
I hope it helps.
Best,
D
You could try 'vegdist()' in the 'vegan' package.
Greg.
On Mar 24, 11:30 pm, Wen Gu wrote:
> Does anyone have a good method for calculating Jaccard coefficients now that
> the dissimilarity() function is no longer an option?
>
> Wen Gu
>
> John Jay College of Criminal Justice445 West 59 StreetN
Hello!
Does exist a command to start a foreign program (Fortran language)? I wrote
a few r-scripts in which data are prepared and several outputfiles are
generated that are used as inputfiles for further calculation. For these
calculations another program is needed, so it would be practically if I
G'day Ken,
On Wed, 25 Mar 2009 00:13:48 -0700 (PDT)
Ken-JP wrote:
>
> Say I have:
>
> > set.seed( 1 )
> > m <- matrix( runif(5^2), nrow=5, dimnames =
> > list( c("A","B","C","D","E"), c("O","P","Q","R","S") ) )
> > m
> O P Q R S
> A 0.2655087 0.898389
Perfect. In conjunction with Jorge's contrib that works a treat. Thanks.
Dan
On Tue, 2009-03-24 at 19:00 -0400, David Winsemius wrote:
> ?try
> ?tryCatch
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read t
Dear all,
I would like to ask information about the package mvtnorm in R.
It is very useful for me the "qmvnorm" comand, but I see that it can
compute only quantile for equicoordinate. Is it possible to have a
curve (or a set of values) if we don't want equicoordinates?
Thank you
Best regards
An
Hi,
I am searching for a way to calculate the power of a Spearman or Kendall
correlation analysis over a range of effect values. Is there any way I can
do this with R?
thanks for the help.
--
View this message in context:
http://www.nabble.com/calculate-the-power-of-kendall-or-spearman-correl
Dear Bernd,
Omitting the NA values from the dataset will work.
ggplot(aes(x = x, color = g), data = na.omit(mydf)) + geom_density()
Notice that I have omitted the group argument. It is redundant in this
case.
HTH,
Thierry
-
so you want to find a needle in a haystack, not an easy task. You should
account for multiple tests, which is as far as I can see not done in the
code yet - or you have to accept that you find a bunch of hay which
accidentally looks pretty much like a needle.
There are some solutions in doing su
Greg gmail.com> writes:
>
> You could try 'vegdist()' in the 'vegan' package.
>
>
Yep, you could. However, if you want to have Jaccard for binary data although
your data are not binary, you must set binary = TRUE in vegan::vegdist.
Indeed, the greatest problem with recommending Jaccard is tha
I'm not really looking for a needle in a haystack, there are a small
number of the 60 tests (about 20) that are likely to concord with other
experiments I have, and in a particular pattern. Since I already have
the data in tables for graphic depiction, I was hoping to have a
reasonably easy way to
Hi Joris,
glm() handles proportions but will give you a warning (and not an error)
about non-integer values. So if you get an error then there should be
something wrong with the syntax, model or data. Can you provide us with
a reproducible example?
Cheers,
Thierry
-
Thank you, Bill, for your answer!
I am also at a total loss when looking for an explanation. I just can't
remember what I did differently...
At least the errors are confined to a rather small dataset so the
repetition of all the glm.nb() analyses won't take much time. The only
thing I found
Hi
I've just upgraded my version of R. I'm used to having the help pages turn up
on a simple text window, and not html. Now it seems that when I call a help
page it starts up a browser, and says it can't find the manual page(probably
since I didn't install html help).
I guess there is some op
On Tue, 2009-03-24 at 11:45 -0700, CE.KA wrote:
> Hi R users,
> I have a big program
> So in Rgui I can t see all the execution of it
> Is there a way to ask R if there is "Errors" in my program
> Sincerely yours
Hi
Normally i use 3 functions in debug R routines: trace, browser and
debug.
In spec
I am using R in C#. I installed R 2.8.1 and R-(D)Com R scilab
DCOM3.0-1B5.exe. I got an error when i run the web application at the line
where R is initilized as
"System.Runtime.InteropServices.COMException was unhandled by user code
Message="Exception from HRESULT: 0x80040013"
Source="Intero
I am trying to carry out discriminant analysis using the syntax:
mydata.lda<-lda(Y~X+Z,data=mydata) and result that I am getting is
Call:
lda(Y ~ X + Z, data = mydata)
Prior probabilities of groups:
0 1
0.3636364 0.6363636
Group means:
XZ
0 17726.750 4020.00
ONKELINX, Thierry schrieb:
Dear Bernd,
Omitting the NA values from the dataset will work.
ggplot(aes(x = x, color = g), data = na.omit(mydf)) + geom_density()
Dear Thierry,
thanks for your reply! Of course, regarding my little toy example
na.omit() perfectly works. However, my real data se
Dear R users
I would like to calculate the Cumulative incidence for an event
adjusting for competing risks and adjusting for covariates. One way to
do this in R is to use the cmprsk package, function crr. This uses the
Fine & Gray regression model. However, a simpler and more classical
appro
Good Morning,
I have a stupid problem about inclusion of greek letters in a
plot built with function pairs().
First of all, I have a matrix with 3 columns
and 1000 rows and I would like to use pairs() in order to plot points in the
upper panel and to compute correlation in the lower panel.
In
Dear Bernd,
If only one variable is causing this problem you could make a subset on
that variable.
data = mydf[!is.na(mydf$g), ]
HTH,
Thierry
ir. Thierry Onkelinx
Instituut voor natuur- en bosonderzoek / Research Inst
Greg,
it seems an obvious behavior to me
y=c(2,2,2,3,3,3,1)
x=1:length(y)
plot(x,y)
lines(x,approxfun(x,y)(x)) # for every x it exists one only value of y
plot(y,x)
lines(sort(y),approxfun(y,x)(sort(y))) # for some y it exists more
than one value of x!
approxfun return a function. By definition
hong shen wrote:
Hi list,
I encountered a situation that a data frame is defined by two packages. Both of
them are loaded by library(). My questions are
1. How could I tell the data frame is from which package?
find("theName")
will tell you where it found a variable called theName.
2. I
PhilippeR wrote:
Dear readers,
I would like to install a local binary package (.tar file) on a Mac (OSX
10.5.6) with R 2.8.0...
How did the package get created? .tar is not an extension used by R.
(Source packages are .tar.gz on all platforms, binary packages are .tgz
on Mac OS.)
Dunca
One further version: this one with a header and with NA's replacing
the -'s that apparently has not deleted any cases with missing data:
http://www.stat.auckland.ac.nz/~wild/764/s764data/prostatic.tab
--
David Winsemius
On Mar 24, 2009, at 11:51 PM, Ravi Varadhan wrote:
Fine detective wor
If you needed to restrict yourself to mvtnorm you could simulate and
get an estimate.
If you were willing to look at other packages you could try the pmnorm
function in mnormt.
On Mar 25, 2009, at 4:00 AM, Antonio Lucadamo wrote:
Dear all,
I would like to ask information about the package
Arafat;
You should listen to Umesh rather than myself on this point, and i
should refrain from posting after my bedtime. I should have typed "glm".
--
David Winsemius
On Mar 25, 2009, at 12:37 AM, Umesh Srinivasan wrote:
> Hi,
>
> Not exactly sure what you mean, but if you want to run a pois
Use the 'reshape' package:
> library(reshape)
> melt(m)
X1 X2 value
1 A O 0.26550866
2 B O 0.37212390
3 C O 0.57285336
4 D O 0.90820779
5 E O 0.20168193
6 A P 0.89838968
7 B P 0.94467527
8 C P 0.66079779
9 D P 0.62911404
10 E P 0.06178627
11 A Q 0.20597457
12
try also
m <- matrix( runif(5^2), nrow=5, dimnames = Names<- list(
c("A","B","C","D","E"),
c("O","P","Q","R","S") ) )
data.frame(expand.grid(Names[[1]],Names[[2]]),as.numeric(m))
data.frame(code=outer(Names[[1]],Names[[2]],paste,sep=".")[1:25],num=as.numeric(m))
Patrizio
2009/3/25 jim holtma
Hi Enrico,
see ?plotmath and the examples
example(plotmath)
there you will find all you need for your task.
hth.
enrico.fosco...@libero.it schrieb:
Good Morning,
I have a stupid problem about inclusion of greek letters in a
plot built with function pairs().
First of all, I have a matrix with
Forgive me for not being more clear.
> Would you expect that one y
> value returns more than one x? I don't
No, I don't either. I want to know the value of that x, however. For
example:
x <- c(2.743, 3.019, 3.329, 3.583, 4.017)
y <- c(0.000, 0.025, 0.025, 0.158, 1.000)
I would like to know th
Rolf Turner wrote:
On 25/03/2009, at 12:09 PM, Frank E Harrell Jr wrote:
(2) Scrolling down to ``Byar and Green prostate cancer data'' appeared
to get
me to the right place. But I couldn't see any signs of any ``R binary
files''.
Please look again. It's under the heading "R". Unfort
Thank you so much both for the answer. I think I have a better handle
on this now. Yes, Loblolly$Seed is an ordered factor, but I didn't
realize that the default for ordered factor is contr.poly.
And then I was further confused because I didn't realize the
coefficient names generated (not
enrico.fosco...@libero.it wrote:
Good Morning,
I have a stupid problem about inclusion of greek letters in a
plot built with function pairs().
First of all, I have a matrix with 3 columns
and 1000 rows and I would like to use pairs() in order to plot points in the
upper panel and to compu
simon bond wrote:
Hi
I've just upgraded my version of R. I'm used to having the help pages turn up on a simple text window, and not html. Now it seems that when I call a help page it starts up a browser, and says it can't find the manual page(probably since I didn't install html help).
I gu
If you want to integrate Fortran into R, see the manual Writing R
Extensions.
If you want to just call some external program (independent of the
language it is written in) via the operating system, see ?system (and
?shell on Windows).
Uwe Ligges
thoeb wrote:
Hello!
Does exist a command t
Bogdan Tanasa wrote:
Hi all,
I am sorry to fill your inbox mail with a naive question, although I would
really appreciate your opinion.
I am looking for a R function that takes a set of aligned sequences, and
transforms a position frequency
matrix into a position weight matrix. Thanks a lot.
Thank you, your answer was extremely helpful. One last problem though: one
of the aggregate functions I'd like to apply on the columns is
concatentation (equivalent to the paste() function). So if I have a given
character column in three separate rows sharing the same ids with the value
"apple" in
- Does it still happen with R-2.9.0 alpha?
- If so, can you please send us reproducible code in order to make able
to see the error on our systems?
Best wishes,
uwe Ligges
Andrew Yee wrote:
Hi, I was wondering if someone in the mailing list has any insight into this
segfault error that I con
Ravi Varadhan wrote:
Fine detective work, David. Now, you can see the reasons for my frustration - multiplicity of data sets combined with non-existent documentation of the source of data in journal articles (e.g. Kay 1986; Lunn and McNeil 1995).
Best,
Ravi.
Yes that is a big frustration
Okay, this one is hard to put into words:
> x <- 1:9; names(x) <- x
> y <- x %o% x
> y
1 2 3 4 5 6 7 8 9
1 1 2 3 4 5 6 7 8 9
2 2 4 6 8 10 12 14 16 18
3 3 6 9 12 15 18 21 24 27
4 4 8 12 16 20 24 28 32 36
5 5 10 15 20 25 30 35 40 45
6 6 12 18 24 30 36 42 48 54
7 7 14 21 28
CORRECTION:
> my.foo <- function( a, b ) { c <- a - b; } #really more complex, but just
> to illustrate
--
View this message in context:
http://www.nabble.com/Doing--o--that-operates-on-columns-instead-of-atomics-tp22701363p22701400.html
Sent from the R help mailing list archive at Nabble.c
In the Links box to the right on the sqldf home page click on
"SQLite - aggregate functions" and lookup group_concat.
On Wed, Mar 25, 2009 at 9:05 AM, Schragi Schwartz
wrote:
> Thank you, your answer was extremely helpful. One last problem though: one
> of the aggregate functions I'd like to appl
Hi,
I have started very recently with R in order to get excellent Box and
Whisker plots. I could plot my data nicely. However, I can't figure
out from R-mailing list archive or google search either, how to place
an Angstrom sign/symbol on the y-axis (any axis in principle), after a
usual y
hacking up on gabor's solution, i've created a trivial function that
will allow you to access a file given a path relative to the path of the
file calling the function.
to be concrete, suppose you have two files -- one library and one
executable -- located in two sibling directories, and you want
I've added an example to FAQ 3 on the home page that
illustrates group_concat.
http://sqldf.googlecode.com
On Wed, Mar 25, 2009 at 9:06 AM, Gabor Grothendieck
wrote:
> In the Links box to the right on the sqldf home page click on
> "SQLite - aggregate functions" and lookup group_concat.
>
> On We
Wacek Kusnierczyk wrote:
> hacking up on gabor's solution, i've created a trivial function that
> will allow you to access a file given a path relative to the path of the
> file calling the function.
>
> to be concrete, suppose you have two files -- one library and one
> executable -- located in tw
Hello, this is a request for assistance that I submitted earlier, this time
with the dataset. My mistake for taking up bandwidth. I've also rephrased
the question to address an additional concern.
I'm working on a windows XP machine with R 2.8.1
1). I'd like a barchart (or other lattice type d
Maybe something like the code below?
There is surely a more elegant way of handling the indices, maybe with plyr?
-s
> xx <- array(1:6,c(3,2)); xx
[,1] [,2]
[1,]14
[2,]25
[3,]36
> yy <- array(1:12,c(3,4)); yy
[,1] [,2] [,3] [,4]
[1,]147 1
Dear R users,
I would like to draw some histograms as seen in the page whose address I
wrote below. I searched through the web a lot and I found a page which
describes how I can do it for older versions of R. For newer versions they
recommend to install the package R.basics in R.clusters but this
Dear all:Request your indulgence. The econophysics gurus do this stuff all
the time: all their PDFs are smooth, with neat log x axis.
1. The kernel density estimate (KDE) function returns the empirical
probability density at 2^n points (min: 512). The big question is how do I
scale back the x-value
Pooja Jain wrote:
Hi,
I have started very recently with R in order to get excellent Box and
Whisker plots. I could plot my data nicely. However, I can't figure
out from R-mailing list archive or google search either, how to place
an Angstrom sign/symbol on the y-axis (any axis in princip
Well I located the error and it had to do with an existing .RData file
that had a print method defined for a class which is part of a package
that wasn't loaded. Starting an interactive session, removing the
objects of the specialized class from the top-level environment, and
re-saving the environ
Dear Evrim,
That is easy to do with the ggplot2 package. You only need the data in a
"long" format. melt() is very usefull to convert data from a wide to a
long format.
library(ggplot2)
n <- 100
Wide <- data.frame(X1 = rnorm(n, mean = -0.5), X2 = rnorm(n, mean = 0,
sd = 2), X3 = rnorm(n, mean = 0
Is this what you want:
> my.foo <- function( a, b ) a - b
> z <- combn(ncol(y), 2) # get all the combinations of 2 columns
> result <- do.call(cbind, lapply(seq(ncol(z)), function(.cols){
+ my.foo(y[,z[1,.cols]], y[, z[2, .cols]])
+ }))
>
>
> result
[,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8]
Google for
R.classes bundle
which will get you to the appropriate page on the author's site where
you can download and install it. Installing the bundle will install
a number of packages including R.basic which contains plot.histogram.
On Wed, Mar 25, 2009 at 10:40 AM, evrim akar wrote:
> Dear
Jose Antonio wrote:
Dear R users,
I have some environmental variables and I need to find the best combination
of them in order to separate two main groups (coded 1 and 2). I have
performed a discriminant analysis using the stepclass function as a method
for selecting the most relevant environm
Thanks Uwe.
Did the following:
ylab=expression(paste("Absolute Error ", (ring(A
Kanu
On 25 Mar 2009, at 14:57, Uwe Ligges wrote:
>
>
> Pooja Jain wrote:
>> Hi,
>> I have started very recently with R in order to get excellent Box
>> and Whisker plots. I could plot my data nicely. However,
I want to do a series of contour plots, 4 in all. The data is coming
from a data frame named "nd.frame", which has elements "xdf", "ydf",
"zdf", and "pndt". I am treating "pndt" as a factor, and it has four
levels. I make a call to the lattice graphics routine "contourplot"
like so:
contourp
On 3/10/09, Saptarshi Guha wrote:
> Hello,
> I have an example of a 2 paneled plot, with two different aspect
> ratios displayed on one page.
> An example would help
>
> n=20
> x1 <- cumsum(runif(n))
> x2 <- cumsum(runif(n))
> d <- data.frame(val=c(x1,x2),id=c(1:n,1:n), nt=c(rep("A",n),rep(
This is not so much a question as a contribution, but comments are welcome.
Comments:
1) thank you very much to Paul Smith in the post
https://stat.ethz.ch/pipermail/r-help/2008-March/157249.html
This is intended to build on that example with something more complex
than
a 2x2 set of
For the first question, add a groups argument. E.g.
barchart(HSI ~ Scenario | Region, Wbirdsm, groups = HydroState)
Also note that using Wbirdsm$HSI makes your call less readable, so I
added the data argument.
For your second question, setting the key does not set the color
theme. You want to se
On 3/25/2009 11:15 AM, Cable, Samuel B Civ USAF AFMC AFRL/RVBXI wrote:
I want to do a series of contour plots, 4 in all. The data is coming
from a data frame named "nd.frame", which has elements "xdf", "ydf",
"zdf", and "pndt". I am treating "pndt" as a factor, and it has four
levels. I make a
Personally I find those types of plots difficult to interpret. Much easier to
create, view, and interpret is to simply plot the lines from density estimates.
See the density function or the logspline package.
--
Gregory (Greg) L. Snow Ph.D.
Statistical Data Center
Intermountain Healthcare
gre
Cable, Samuely hanscom.af.mil>
> I want to do a series of contour plots, 4 in all. The data is coming
> from a data frame named "nd.frame", which has elements "xdf", "ydf",
> "zdf", and "pndt". I am treating "pndt" as a factor, and it has four
> levels. I make a call to the lattice graphics r
Hi
I am plotting a set of data inside a for loop.
Is it possible to use plot in for loop without redrawing the whole
plot? Am using par(new=TRUE) but that draws on top of the previous plot.
Couldn't find any threads about the topic.
Thanks
Mohan
--
___
Dear all,
I am trying to manually re-sort rows in a number of tables. The rows aren't
sorted on any particular values but are simply ordered by user choice (as shown
by the row numbers in the code). I have been able to carry out each
re-arrangement without the use of the 'for' loop, but cannot
Dear R users,
I am running R-2.8.1 with wavethresh-2.2-11
I have a problem with the function wr, used for reconstructing wd objects. The
problem is that I cannot set start.level to any integer value greater than 1.
Whenever I set it to any value greater than 1, I get:
Building level: 1 re
Can't make sense of calculated results and hope I'll find help here.
I've collected answers from about 600 persons concerning three
variables. I hypothesise those three variables to be components (or
indicators) of one latent factor. In order to reduce data (vars), I
had the following idea:
I am afraid your notion of a "concrete idea" is less concrete than
what I would need to understand what you are requesting. Your first
lines of example code should be:
library()
density(x, )
... since stats::density() by default would return 512 y estimates,
even if the length of x we
Or use frequency polygons, if you want to stay with the
interpretability of a histogram.
Hadley
On Wed, Mar 25, 2009 at 12:07 PM, Greg Snow wrote:
> Personally I find those types of plots difficult to interpret. Much easier
> to create, view, and interpret is to simply plot the lines from densi
It's unclear to me what your expected output is. If you are trying to
add additional data sets to an existing plot then ?lines should be
sufficient.
-Original Message-
From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org]
On Behalf Of Mohan Singh
Sent: Wednesday, March 2
Aren't you missing a sep='' in your last call to paste?
-Original Message-
From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org]
On Behalf Of Steve Murray
Sent: Wednesday, March 25, 2009 1:58 PM
To: r-help@r-project.org
Subject: [R] Manual sort in a for loop
Dear all,
On 03/25/09 19:06, soeren.vo...@eawag.ch wrote:
> Can't make sense of calculated results and hope I'll find help here.
>
> I've collected answers from about 600 persons concerning three
> variables. I hypothesise those three variables to be components (or
> indicators) of one latent factor. In
I assume you need to use 'get' to retrieve the value:
table_year=1951
for (i in (paste("arunoff_",year,"_temp",sep=""))) {
assign(paste("arunoff_",table_year,
sep=""),get(paste("arunoff_",table_year,"_temp"))[c(10,7,9,5,4,12,1,3,2,8,11,6),])
table_year = table_year+1
}
On Mar 25, 2009, at 8:37 AM, Greg wrote:
Forgive me for not being more clear.
Would you expect that one y
value returns more than one x? I don't
No, I don't either. I want to know the value of that x, however. For
example:
x <- c(2.743, 3.019, 3.329, 3.583, 4.017)
y <- c(0.000, 0.025, 0.
Hii,
Is it possible in R to write out to the graph the exact values of the
X-Axis, even if the values are high ?
I will avoid to get the values in this Form for example: 3e+06
Can anybody help ?
--
View this message in context:
http://www.nabble.com/diversify-the-values-of-the-x-Axis...-tp227
well, the literal answer is that paste("arunoff_",table_year,"_temp")
is a character vector of length 1 so your indexing cannot work. What
you want is to index the data that corresponds to this variable name,
?get
But I should stress that this manipulation with assign and get seems
complete
Hi, Brain,
I have a similar question. lines() can only add lines to plot(). If I use
xyplot(), how do I add lines from a different dataset? Thanks.
Jun
On Wed, Mar 25, 2009 at 1:18 PM, Rowe, Brian Lee Yung (Portfolio Analytics)
wrote:
> It's unclear to me what your expected output is. If you a
Hi,
I have the following type of data:
myData <- data.frame(x = 1:5, y = letters[1:5], xDate =
seq(as.Date("2001/2/1"), as.Date("2005/2/1"), by="year") )
> myData
x y xDate
1 a 2001-02-01
2 b 2002-02-01
3 c 2003-02-01
4 d 2004-02-01
5 e 2005-02-01
What I need is a new column, s
I understood this to mean you want to open a new plotting window on each
iteration of the loop. If this is correct, I usually go about it by
using x11()
If you're looking to add additional lines or points, then you may want
to look at the aptly named functions lines() and points().
If neither of
Dear all,
Trying to extract a few rows for each element of a list of
data.frames, I'm puzzled by the following behaviour,
d <- lapply(1:4, function(i) data.frame(x=rnorm(5), y=rnorm(5)))
str(d)
lapply(d, "[", i= c(1)) # fine, this extracts the first columns
lapply(d, "[", j= c(1, 3)) #
Hello R experts,
I went through R mailing,Nabble R.I could not find solution.Can someone
help me.
I have undirected Graph.
Here is an example of spreadsheet I have( Unique 3559 Nodes)
snippet of 4 rows.
Node1 Node2 Weights
1 2 5
2 3 30
2 4 30
1
basically the for loop goes something like this
setCounters(135);
for(i in 1:query) {
plot(c2data[start1:count,],c3data[start1:count,],xlim=c(-30,100),
ylim=c(-30,90), sub=i);
points(..);
#count increments
}
so if i use windows() or x11(), i get different plot windows
if I use par(),
This might not be particularly elegant, but if you put your initial
plot() outside of the loop and then use the loop to place your points,
you might get what you want. If others have better solutions, I'd be
interested as well. Note that if you take this approach, you might want
to specify xaxt="
Hello Nathan,
a couple of points.
1) You did not write it, but it seems that you use the igraph package.
It might be worth pointing this out, and also, to post to the
igraph-help mailing list if you don't get any answer here.
2) igraph can easily create graphs from edge lists, there is no need
t
Hi Simon,
the help-pages for ?options don't list the properties chmhelp and
htmlhelp if I remember correctly. But these are the two options you
should set to FALSE if you don't want HTML-Help (nor CHM-based help on
Windows).
options(htmlhelp=FALSE)
options(chmhelp=FALSE)
Regards,
Jannis
__
Thanks for your reply.
But for some packages, it works. For others, you may get
Error: 'theName' is not an exported object from 'namespace:A'
What does "exported object" mean?
On Mar 25, 6:48 am, Duncan Murdoch wrote:
> hong shen wrote:
> > Hi list,
>
> > I encountered a situation that a data
A::thename works only when thename is exported by the name space A.
A:::thename works only when thename is defined in the name space A.
getAnywhere(thename)$objs$"package:A" always works.
On Mar 25, 6:48 am, Duncan Murdoch wrote:
> hong shen wrote:
> > Hi list,
>
> > I encountered a situation th
Hi Folks,
Some of you (since you have contributed) are aware of a quite
vigorous discussion currently in progress on the MedStats
Google Group. Others, who possibly could contribute usefully,
may not be.
For the moment it is at the top of the "Discussions" list at
http://groups.google.com/group
not really, as
plot(NA,NA,xlim=c(-30,100), ylim=c(-30,90), ...);
is the main plotting statement inside the loop, draws different
scatterplots for each loop
i can draw points or lines inside for loop, but not the
plot(NA,NA,xlim=c(-30,100), ylim=c(-30,90), ...);
statement,
it draws on top of pr
Hi Sören,
>> (1) Is there an easy example, which explains the differences between
>> pca and pfa? (2) Which R procedure should I use to get what I want?
There are a number of fundamental differences between PCA and FA (Factor
Analysis), which unfortunately are quite widely ignored. FA is expli
What am I missing here? Both of my column headers in the 'FREQ' table are found
in the 'genotype'; however, they aren't being recognized.
> colnames(FREQ)
[1] "X17362526" "X17362627"
> colnames(genotype)
[1] "X17362311" "X17362316" "X17362346" "X17362420" "X17362421" "X17362422"
"X17362435" "X
Dear Jacy,
Try this:
colnames(FREQ) %in% colnames(genotype)
See ?"%in%" for details. Also, take a look at also ?match.
HTH,
Jorge
On Wed, Mar 25, 2009 at 3:49 PM, Crosby, Jacy R
wrote:
> What am I missing here? Both of my column headers in the 'FREQ' table are
> found in the 'genotype'; how
Wacek, this work for me. Takk!
Mvh.
Marie
On Wed, Mar 25, 2009 at 2:59 PM, Wacek Kusnierczyk <
waclaw.marcin.kusnierc...@idi.ntnu.no> wrote:
> Wacek Kusnierczyk wrote:
> > hacking up on gabor's solution, i've created a trivial function that
> > will allow you to access a file given a path relat
Hi,
I am a biologist (relatively new to R) analyzing data which we predict
to fit a power function. I was wondering if anyone knew a way to model
piecewise functions in R, where across a range of values (0-x) the data
is modeled as a power function, and across another range (x-inf) it is a
li
1 - 100 of 142 matches
Mail list logo