Dear all,
I do not fully understand how ggplot2 handles NAs. See the following
example:
library(ggplot2)
x - rnorm(150)
g - as.factor(c(rep(c(0,1,NA),50)))
mydf - data.frame(x,g)
m - ggplot(aes(x = x, group = g, color = g), data = mydf)
m + geom_density()
How do I get rid of the NAs (i.e.
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 0.5728534 0.66079779
hong shen hshen_1998 at 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,
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
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=NULL)
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,
Dimitris
Ken-JP wrote:
You could try 'vegdist()' in the 'vegan' package.
Greg.
On Mar 24, 11:30 pm, Wen Gu edwardg22...@hotmail.com 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
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
G'day Ken,
On Wed, 25 Mar 2009 00:13:48 -0700 (PDT)
Ken-JP kfmf...@gmail.com 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.89838968 0.2059746
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
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
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:
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
Greg chaoborid at 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
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
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
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 special
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
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
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
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.
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
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. If
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.)
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
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
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 poisson
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 holtman jholt...@gmail.com:
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
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 the
Rolf Turner wrote:
On 25/03/2009, at 12:09 PM, Frank E Harrell Jr wrote:
snip
(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.
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
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
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
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
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
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 35 42
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.com.
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
schra...@post.tau.ac.il wrote:
Thank you, your answer was extremely helpful. One last problem though: one
of the aggregate
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
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
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
ggrothendi...@gmail.com wrote:
In the Links box to the right on the sqldf home page click on
SQLite - aggregate functions and lookup
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 two
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
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 10
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
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
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
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 =
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] [,9]
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
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
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 can't
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:
On 3/10/09, Saptarshi Guha saptarshi.g...@gmail.com 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),
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
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
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 call to
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
Cable, Samuely at 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 routine
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
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
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
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(if the density function is not from stats)
sample data construction of x
density(x, some set of parameters)
...
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 greg.s...@imail.org wrote:
Personally I find those types of plots difficult to interpret. Much easier
to create, view, and interpret is to simply plot
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
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
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 Wed, Mar
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,
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:
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
completely
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)
b_r...@ml.com wrote:
It's unclear to me what your expected output
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, say
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
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
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
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 murd...@stats.uwo.ca wrote:
hong shen wrote:
Hi list,
I encountered a
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 murd...@stats.uwo.ca wrote:
hong shen wrote:
Hi list,
I encountered a
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
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
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
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 X17362438 X17362459
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
jacy.r.cro...@uth.tmc.eduwrote:
What am I missing here? Both of my column headers in the 'FREQ' table are
found
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 relative to
1 - 100 of 142 matches
Mail list logo
