Max ha scritto:
> Hi everyone!
>
> I've recently begun to learn R for my job as the IT department suffers
> from lack of funding for new software. I was talking to the guy in
> charge of Requisitions and have found out the budget for books is in
> great shape.
>
> So, I'm curious what books peop
Dear all:
"predict.glm" provides an example to perform logistic regression when the
response variable is a tow-columned matrix. I find some paradox about the
degree of freedom .
> summary(budworm.lg)
Call:
glm(formula = SF ~ sex * ldose, family = binomial)
Deviance Residuals:
Min
Hi
Well, R does exactly what it says. From help page.
"Otherwise, x and y must be vectors or factors of the same length"
I do not know SAS but I presume that
> tables bloodtype*state
gives you something like
tab <- table(bloodtype, state)
and
chisq.test(tab)
shall give you the expected res
정 태훈 gmail.com> writes:
> I got a reply for my previous several postings saying that I was
> spamming the r-help mailing list.
> I am very sorry to all subscribers if I did that.
> But I've been reposting my message to the mailing list several times
> because I didn't know whether my help pos
Hi all
I am trying to perform the follwing:
fit<-pvclust(wq, method.hclust="ward", method.dist="euclidean")
but get a strange error message that I just cant figure out.
Has anyone come across this? Any help would be most appricieated
Error in hclust(distance, method = method.hclust) :
NA/NaN/
Hi
No error message? It seems to me that it is OK. There is of course limited
width for printout on screen (see GUI preferences/console columns), and
therefore everything which exceeds specified limits is printed on second
and subsequent row.
try
dim(form1) or str(form1) for evaluation of you
Dear all:
I did the following because I was not sure what a logistic regression
fits ."small.glm" is a glm fit , setting the family to be binomial.
>fitted(small.glm)->p
>log(p/(1-p))->left
> for(i in 1:length(coef(small.glm))){
+ coef(small.glm)[i]*model.matrix(small.glm)[,i]->res[,i]
+
The message seems to be messed up but if your data has one more column
than the header then it will assume the first column is the row names.
Use the R count.fields function to diagnose this.
On Dec 7, 2007 1:24 AM, <[EMAIL PROTECTED]> wrote:
> I am using read.table and keep on getting this mess
I am using read.table and keep on getting this message.
The function is confusing my first column with a row.names column. I have
checked the table carefully using excel and it seems quite symetric and
with a name at the top of each column.
This is what I am using:
form1 = read.table("c:/horses
Try this:
matplot(xxx[,1], xxx[2:4], type = "l")
with(xxx, {
segments(x, z1, x, z2)
points(x, z1)
points(x, z2)
})
Omit the two points commands if you don't want circles at the ends of
the segments.
On Dec 6, 2007 6:30 PM, David Rees <[EMAIL PROTECTED]> wrote:
> Hi,
>
> I
On 06/12/2007, Loren Engrav <[EMAIL PROTECTED]> wrote:
> As for news readers
> I found R and R.mac and R.Bio on the sites you recommend, thank you very
> much, they would avoid the individual emails, but then I would have to go
> look at them, which might be Ok
>
> Deluge? Well, there are from R an
Ryan Briscoe Runquist wrote:
>
>
> Hello,
>
> I am using a simple linear model and I would like to get an AIC value. I
> came across both AIC() and extractAIC() and I am not sure which is best to
> use. I assumed that I should use AIC for a glm and extractAIC() for lm,
> but if I run my mod
Peter Dalgaard <[EMAIL PROTECTED]> wrote in
news:[EMAIL PROTECTED]:
> The passage you cite from the manual could do with a
> rephrasing, although it probably isn't technically incorrect. As it
> stands, it reminds me a bit of the old Monty Python sketch:
>
> "Our *three* weapons are fear, surpris
Hello Vito,
Thanks for your reply and apologies for not being clearer.
I'd like to fit a three-segmented relationship to each level but have only 3
unique breakpoints. The result would be 9 slopes, one of which would be zero.
I achieved this by finding the 3 breakpoint with:
init.bp <- c(297.
Hi,
I'm using natural cubic splines from splines::ns() in survival
regression (regressing inter-arrival times of patients to a queue on
queue size). The queue size fluctuates between 3600 and 3900.
I would like to be able to run predict.survreg() for sizes <3600 and
>3900 by assuming that the
Hi, all;
I got a reply for my previous several postings saying that I was
spamming the r-help mailing list.
I am very sorry to all subscribers if I did that.
But I've been reposting my message to the mailing list several times
because I didn't know whether my help post was actually posted or n
As for news readers
I found R and R.mac and R.Bio on the sites you recommend, thank you very
much, they would avoid the individual emails, but then I would have to go
look at them, which might be Ok
Deluge? Well, there are from R and Bio and R-Mac every morning 30 or 35, and
10-15 more during the
Hello,
I am using a simple linear model and I would like to get an AIC value. I
came across both AIC() and extractAIC() and I am not sure which is best to
use. I assumed that I should use AIC for a glm and extractAIC() for lm,
but if I run my model in glm the AIC value is the same if I use AIC(
You can use the socket implementation of snow, with
library(snow)
cl <- makeCluster(rep("localhost", 2), type="SOCK")
to start up a cluster of 2 R processes.
luke
On Thu, 6 Dec 2007, Saeed Abu Nimeh wrote:
> Hi,
> I installed the snow package on a unix box that has multiple cores. To
On Dec 6, 2007 6:13 PM, Emmanuel Charpentier <[EMAIL PROTECTED]> wrote:
> Gabor Grothendieck a écrit :
> > On Dec 6, 2007 3:26 PM, marciarr <[EMAIL PROTECTED]> wrote:
>
> [ Snip ]
>
> >> 2- how do I solve a simple equation? Considering the equation y=
> >> exp(-x)^12,
> >> I would like to find
Thomas Pujol a écrit :
> R-help users,
> Thanks in advance for any assistance ... I truly appreciate your expertise.
> I searched help and could not figure this out, and think you can probably
> offer some helpful tips. I apologize if I missed something, which I'm sure I
> probably did.
>
Gabor Grothendieck a écrit :
> On Dec 6, 2007 3:26 PM, marciarr <[EMAIL PROTECTED]> wrote:
[ Snip ]
>> 2- how do I solve a simple equation? Considering the equation y= exp(-x)^12,
>> I would like to find the values of x for, for example, y=0.01, so
>> exp(-x)^12=0.01. How do I do that using
Hi,
I am having trouble plotting the graph I need given the follow kind of data
> xxx <- data.frame(
"x"=c(1,2,3,4,5),
"y1"=c(2,4,3,5,6),
"y2"=c(3,4,6,3,1),
"y3"=c(1,3,5,7,3),
"z1"=c(1,NA,3,5,NA),
"z2"=c(2,N
If you are using namespaces see:
?assignInNamespace
On Dec 6, 2007 5:38 PM, Tao Shi <[EMAIL PROTECTED]> wrote:
>
> Hi list,
>
> Sorry for the vague title, but here is the scenario.
>
> I'm writing an R package, let's say, 'pkg1', which contains 3 functions: f1,
> f2, f3. f2 and f3 are helper fu
On Dec 6, 2007 3:26 PM, marciarr <[EMAIL PROTECTED]> wrote:
>
> Hello R users,
> I have been looking through Help files and Nabble list for the answers for
> these simple questions, but it seems to be fruitless.
> 1- in a data frame with two columns, x and y, how do I get the corresponding
> value
Hi list,
Sorry for the vague title, but here is the scenario.
I’m writing an R package, let’s say, ‘pkg1’, which contains 3 functions: f1,
f2, f3. f2 and f3 are helper functions for f1, i.e. f1 calls f2 which in turn
calls f3.
f1 <- function(…) {
….
f2()
…
}
f2 <-
The chi-square does not need your two categorical variables to have equal
levels, nor limitation for the number of levels.
The Chi-square procedure is as follow:
χ^2=∑_(All Cells)▒〖(Observed-Expected)〗^2/Expected
Expected Cell= E_ij=n((i^th RowTotal)/n)((j^th RowTotal)/n)
Degree of Freedom=df=
The chi-square does not need your two categorical variables to have equal
levels, nor limitation for the number of levels.
The Chi-square procedure is as follow:
χ^2=∑_(All Cells)▒〖(Observed-Expected)〗^2/Expected
Expected Cell= E_ij=n((i^th RowTotal)/n)((j^th RowTotal)/n)
Degree of Freedom=df=
The chi-square does not need your two categorical variables to have equal
levels, nor limitation for the number of levels.
The Chi-square procedure is as follow:
χ^2=∑_(All Cells)▒〖(Observed-Expected)〗^2/Expected
Expected Cell= E_ij=n((i^th RowTotal)/n)((j^th RowTotal)/n)
Degree of Freedom=df=
Dear list,
I am wondering if there are any parameters in heatmap.2
to be able to adjust the size of row labels. Or the size
of the plot can be made smaller. The plot I got is big
and cannot see the full row labels completely.
And what is the maximum of rows which could allow row
labels visible in
marciarr wrote:
> Hello R users,
> I have been looking through Help files and Nabble list for the answers for
> these simple questions, but it seems to be fruitless.
> 1- in a data frame with two columns, x and y, how do I get the corresponding
> value of x to, let's say, the minimum value of the y
Is the cross validation procedure implemented in lda() from the MASS
package internal or external?
Thanks,
Bijun Tan
Cooper Union
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://
Hello:
I would like to run a Cox proportional hazards regression on crayfish
dislodgement at different water velocities by crayfish size class and
substrate (rock) type. Additionally, there is a covariate variable, rock
movement that may be influencing crayfish dislodgment. So...
I have crayfi
Hi,
Could somebody help me with the following. I want to calculate the integral
over an implicit function. I thought to integrate over a function depending on
uniroot. In previous topics I found a thread about finding the root of an
integral. And that works. But the other way around, does not w
Dear David,
You may also consider using the biocep project' tools and frameworks. they
provide an advanced bridge that allow you to exchange between R and Java
any standard R Object and any mapped S4 object.
the object extracted to Java (an RList for you data) can be serialized to
a file (saved a
I'd say that was pretty clearly a bug in stats4 (which as I recall was
needed to get around the scoping awkwardnesses of S4).
But could you not write a logLik method for your class? E.g.
> logLik.bar <- function(object, ...) structure(pi, class="logLik", df=1)
> AIC(bar.tmp)
[1] -4.283185
> lib
I have learned something new - thanks for the strwrap info.
The problem with posting from Nabble is that by the time your post actually
gets to the list (2 hours after you posted it in this case) and you've
written some line like "Knowing the R list, someone can probably reduce this
function to 2
Happy Holidays Folks!
XLSolutions Corporation has scheduled 2 great R/Splus Advanced
programming courses. More course can be viewed on our website.
R/Splus Advanced Programming* San Francisco January 21-22,
2008
R/Splus Advanced Programming* New York City January 24-25,
2
Hi,
Is there an R library that has the same functionalities of Splus7.0+ library
correlatedData?
I'd appreciate any input.
Hakan Demirtas
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/
Hi,
I am trying to use stepAIC to do forward stepwise selection for a logistic
model.
I did
hi <- glm(melanoma ~ ., family = binomial, data=data)
lo <- glm(melanoma ~ 1, family = binomial, data=data)
stepAIC(lo, scope=list(upper = hi , lower = lo), direction = "forward" )
But the full m
On Thu, 2007-12-06 at 14:00 -0500, Gerard Tromp wrote:
> Greetings!
>
> Recently when I tried to use lm.influence I get the following error:
> Error in .Fortran("lminfl", model$qr$qr, n, n, k, as.integer(do.coef), :
> Fortran symbol name "lminfl" not in DLL for package "base"
>
> This occurs
On Thu, 2007-12-06 at 09:45 -0800, mogra wrote:
> I have huge data file, and I would like randomize just one column at a time ,
> is there any easy way?
>
> Thanks a lot.
If you just want to randomly sample <\hint> from a single column,
_independent of the other columns_, you can use:
DF$Col
On Thu, 2007-12-06 at 09:28 -0800, Jim Price wrote:
> I wrote a little utility function for exactly this reason, which I use with
> long titles. You may want to add calls to par to adjust the upper margin if
> you are using raw graphical functionality (plot et al) - but lattice adjusts
> the upper
?table should work
table (DMS, ITS)
I am not clear on what kind of chart you want.
will
plot(DMS, ITS) do what you want?
--- christopher snow <[EMAIL PROTECTED]> wrote:
> I have a dataset with two variables that are
> factors:
>
> 1) Decision Making Satisfaction (DMS), values = A -
> Compl
Hello R users,
I have been looking through Help files and Nabble list for the answers for
these simple questions, but it seems to be fruitless.
1- in a data frame with two columns, x and y, how do I get the corresponding
value of x to, let's say, the minimum value of the y column (min (data$y)) ?
On 12/5/07 12:15 PM, "Prof Brian Ripley" <[EMAIL PROTECTED]> wrote:
On Wed, 5 Dec 2007, David Coppit wrote:
> Or, given that I'm dealing with just a single array, would it be better to
> roll my own I/O using write.table or write.matrix from the MASS package?
It would be much easier. The save(
On Thu, 2007-12-06 at 12:51 +, christopher snow wrote:
> I have a dataset with two variables that are factors:
>
> 1) Decision Making Satisfaction (DMS), values = A - Completely, B -
> Mostly, C - Partly, D - Not at all
> 2) IT Satisfaction values (ITS), values = A - Completely, B - Mostly,
Hello:
How can I work around the conflict between the S3 and S4
illustrated in the example below? I'm developing a package that
requires a function in 'stats4', but when 'stats4' is attached, it
breaks my AIC function. I could give my AIC function another name so it
no longer uses the
The 'table' function will give you the simple counts. Plotting a table
with the 'plot' function gives common charts for this. The 'CrossTable'
function in the gmodels package creates the table along with additional
information. There are a lot of other functions for creating/working
with tables
Hello,
I am not only interested in finding out which genes are the most highly up-
or down-regulated (which I have done using the linear models and Bayesian
statistics in Limma), but I also want to know which genes are consistently
highly transcribed (ie. they have a high intensity in the chann
\n to start a new line
plot(1:10,main="This is my really long title \n and
it's so long that I can see \n just about half of
it.")
--- Svempa <[EMAIL PROTECTED]> wrote:
>
> I want to fit a fairly long main title for a plot,
> supposedly by changing row
> after a while. As for now it starts way
Try This:
plot(A,main=paste("This is my really long title and","\n","it's so long that
I can see
just about half of it.", sep = " "))
-- Bert Gunter
Genentech
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do
Try this:
df <- data.frame(DMS=factor(rep(LETTERS[1:4], 10)),
ITS=factor(rep(LETTERS[1:4], 10)))
table(df)
plot(table(df))
--
Henrique Dallazuanna
Curitiba-Paraná-Brasil
25° 25' 40" S 49° 16' 22" O
On 06/12/2007, christopher snow <[EMAIL PROTECTED]> wrote:
> I have a dataset with two variables
Greetings!
Recently when I tried to use lm.influence I get the following error:
Error in .Fortran("lminfl", model$qr$qr, n, n, k, as.integer(do.coef), :
Fortran symbol name "lminfl" not in DLL for package "base"
This occurs on both Linux and Windows platforms (details below).
Searching the ma
op <- par(mar=(c(5, 4, 4, 4) + 0.1))
plot(1:25,runif(25,0,1),ylab="First Y-axis
label",xaxt="n")
axis(4,at=seq(0.2,1,.2), labels=1:5)
mtext("Second Y-axis label",side=4, line=2)
par(op)
--- Marcin Kozak <[EMAIL PROTECTED]> wrote:
> Hi,
>
> Consider this simple plot:
> > plot(1:25,runif(25,0,1)
Michael Kubovy virginia.edu> writes:
>
> This is one of the cases where I wish there were a function that
> "exercised" all the arguments of a graphics function by visualizing
> the effect of changing two or three levels of each argument (one by
> one, of course). This might have the side
I have huge data file, and I would like randomize just one column at a time ,
is there any easy way?
Thanks a lot.
--
View this message in context:
http://www.nabble.com/Randomizing-one-column-in-the-dataMatrix-tf4957535.html#a14197423
Sent from the R help mailing list archive at Nabble.com.
_
I wrote a little utility function for exactly this reason, which I use with
long titles. You may want to add calls to par to adjust the upper margin if
you are using raw graphical functionality (plot et al) - but lattice adjusts
the upper margin automatically so you wouldn't need to add anything e
On Thu, 2007-12-06 at 09:27 -0800, qian z wrote:
> I used merge() in a function, but it doesn't return correct data frame.
>
> add.name <- function(data, x)
>
> {
>
> ...
> ...
>
> newfile <- merge(data, resid, by =0, all.x=TRUE, all.y= FALSE)
> newfile
>
>
>
try this:
plot(0, main=paste(strwrap("This is my really long title and it's so
long that I can see just about half of it.", width=50),
collapse="\n"))
On Dec 6, 2007 7:16 AM, Svempa <[EMAIL PROTECTED]> wrote:
>
> I want to fit a fairly long main title for a plot, supposedly by changing row
> aft
On Thu, 2007-12-06 at 07:16 -0800, Svempa wrote:
> I want to fit a fairly long main title for a plot, supposedly by changing row
> after a while. As for now it starts way outside the picture margin at the
> left and continues way out right passed the right margins.
>
> >plot(A,main="This is my re
Try inserting \n into the title where you would like it to start on a
new line, e.g.:
plot(A,main="This is my really long title\nand it's so long\nthat I can
see just about half of it.")
You may need to give yourself more room in the margin for multi-line
titles (see ?par and the mar entry) and u
I used merge() in a function, but it doesn't return correct data frame.
add.name <- function(data, x)
{
...
...
newfile <- merge(data, resid, by =0, all.x=TRUE, all.y= FALSE)
newfile
}
-
[[alternative H
This will give you the percents in the same order as your original data (as
this is what your original code did)
apply(tdat, 2,
function(x) {
o <- order(x)
oldo <- order(o)
prc <- cumsum(x[o]) / sum(x)
prc[oldo]
})
Jason Law
Statistician
City of Portland, Bureau of
Hi,
I need to fit a general generalized linear model, for observations
on a response Y which is Gamma-distributed and observed in clusters. So, if
E(Y) = mu,
then for a suitable link function f, f(mu) is a (linear) function of
both fixed and random effects. Are there R packages or fu
I want to fit a fairly long main title for a plot, supposedly by changing row
after a while. As for now it starts way outside the picture margin at the
left and continues way out right passed the right margins.
>plot(A,main="This is my really long title and it's so long that I can see
just about
I have a dataset with two variables that are factors:
1) Decision Making Satisfaction (DMS), values = A - Completely, B -
Mostly, C - Partly, D - Not at all
2) IT Satisfaction values (ITS), values = A - Completely, B - Mostly, C
- Partly, D - Not at all
I would like to produce a table (matrix)
Todd Remund hotmail.com> writes:
> I know that there is a function, (spectro3D), that produces the Power
Spectrogram. Are there R functions
> that produce the Frequency Spectrogram and the Phase Spectrogram? Thank you
for your time.
fft in stats gives you all you need, possibly combined with s
Thanks Ben, that fixed the display within R. However, even after changing
the display settings, the matrix elements still appear to be exported in
single precision. The matrix object is being passed into my C routines as
an SEXP Numeric type, and somewhere along the way, some of the digits are
g
Hi,
I am trying to use stepAIC to do forward stepwise selection for a logistic
model.
I did
hi <- glm(melanoma ~ ., family = binomial, data=data)
lo <- glm(melanoma ~ 1, family = binomial, data=data)
stepAIC(lo, scope=list(upper = hi , lower = lo), direction = "forward" )
But the full mod
Thank you all. That's exactly what I was looking for.
TLowe wrote:
>
> Hey Folks,
>
> Could somebody help me rewrite the following code?
>
> I am looping through all records across 5 fields to calculate the
> cumulative
> percentage of each record (relative to each individual field).
>
> I
Is this basically what you want to do? (Please include commented,
minimal, self-contained, reproducible code so we don't have to guess
at what you want)
> x <- data.frame(a=runif(10), b=runif(10))
# do for one column
> cumsum(x$a)/sum(x$a)
[1] 0.05892073 0.08129611 0.11067218 0.28640268 0.2896982
How about this example?
## sample data frame with two columns
df <- data.frame(x = abs(rnorm(20)), y=abs(rnorm(20,2)))
## create new variables in df with an lapply call
df[c("cpctx","cptcty")] <- lapply(df, function(x) cumsum(x)/sum(x))
A possible improvement would be to construct the new column
On Wed, 2007-12-05 at 22:33 -0800, Bin Yue wrote:
> Dear all:
> By comparing glmresult$y and model.response(model.frame(glmresult)), I
> have found out which one is
> set to be "TRUE" and which "FALSE".But it seems that to fit a logistic
> regression , logit (or logistic) transformation has
On 12/6/2007 8:01 AM, vito muggeo wrote:
> Dear all,
> Is there *any* reason explaining what I describe below?
Probably, but you're unlikely to get a correct explanation if you don't
give us code to reproduce the problem.
From the look of things, your real example is quite complex, which is
pr
Thanks a lot for the quick pointer, Richie. I will take a close look of the
RTisean package.
Gang
From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED]
Sent: Thu 12/6/2007 11:46 AM
To: Chen, Gang (NIH/NIMH) [C]
Cc: r-help@r-project.org
Subject: Re: [R] Any package
> -Original Message-
> From: [EMAIL PROTECTED]
> [mailto:[EMAIL PROTECTED] On Behalf Of Michael Kubovy
> Sent: Thursday, December 06, 2007 2:12 AM
> To: Dieter Menne
> Cc: [EMAIL PROTECTED]
> Subject: Re: [R] Using expression in Hmisc Key()
>
[snip]
> This is one of the cases where I wish
R-help users,
Thanks in advance for any assistance ... I truly appreciate your expertise.
I searched help and could not figure this out, and think you can probably offer
some helpful tips. I apologize if I missed something, which I'm sure I probably
did.
I have data for many "samples".
It sounds like you should sample x and y together using the
block bootstrap. If you have the usual situation, x and y in columns
and observations in rows, then sample blocks of rows.
Even though observations in y are independent, you would take
advantage of that only for bootstrapping statistics
The error message says you have duplicate row names and that
is not allowed. Make sure you have the same number of elements
on each line of data as in the header. If you have one more on each line
than on the header then the first data item on each line will be regarded
as the row name. See ?cou
Hey Folks,
Could somebody help me rewrite the following code?
I am looping through all records across 5 fields to calculate the cumulative
percentage of each record (relative to each individual field).
Is there a way to rewrite it so I don't have to loop through each individual
record?
# t
On 12/6/07, Christopher Oezbek <[EMAIL PROTECTED]> wrote:
> But is there a technical reason for making this distinction? It seemed to
> me that most panel functions can deal with x and y parameters and that
> panel.densityplot could make use of the same mechanisms as panel.violin to
> subdivide ba
--- begin included message
I am running R 2.6.1 on windows xp
I am trying to fit a cox proportional hazard model with a shared
Gaussian frailty term using coxme My model is specified as:
nofit1<-coxme(Surv(Age,cen1new)~ Sex+bo2+bo3,random=~1|isl,data=mydat)
With x1-x3 being dummy variables, and i
hi gabor,
i was able to get your suggestion to work. i have been going through
the R help tools to figure out what each step actually does because i
have something similar but hours 2,5,8,11,14,17 and 20 are missing. i
haven't had any luck. each "mean value" that is calculated is the
same.
The RTisean package has wiener filter functions (wiener1 and wiener2).
Regards,
Richie.
Mathematical Sciences Unit
HSL
"Chen, Gang (NIH/NIMH) [C]" <[EMAIL PROTECTED]>
Sent by: [EMAIL PROTECTED]
06/12/2007 16:20
To
cc
Subject
[R] Any package for deconvolution?
I want to run deconvolu
Hi,
I installed the snow package on a unix box that has multiple cores. To be
able to exploit the multiple cores (on one pc) do I still need to install
the rmpi package (or rpvm). Another question, if i run a bayesian simulation
on the multiple core after setting them up correctly (using snow), wou
I want to run deconvolution of a time series by an impulse or point-spread
function through Wiener filter, regularized filter, Lucy-Richardson method, or
any other approaches. I searched the CRAN website and the mailing list archive,
but could not find any package for such a deconvolution analys
Thanks James!
Allen
On Dec 6, 2007 9:21 AM, James W. MacDonald <[EMAIL PROTECTED]> wrote:
>
>
> affy snp wrote:
> > Dear list,
> >
> > I am using heatmap.2(x) to draw a heatmap. Ideally, I want to the matrix
> > x clustered only by columns and keep the original order of rows
> unchanged.
> > Is
Thanks Ashoka!
Allen
On Dec 6, 2007 12:27 AM, Ashoka Polpitiya <[EMAIL PROTECTED]>
wrote:
> Check the Rowv, Colv options to heatmap.2
>
> data(mtcars)
> x <- as.matrix(mtcars)
> heatmap.2(x, Rowv=FALSE, dendrogram="column")
>
> -Ashoka
>
> Scientist - Pacific Northwest National Lab
>
> On
I know that there is a function, (spectro3D), that produces the Power
Spectrogram. Are there R functions that produce the Frequency Spectrogram and
the Phase Spectrogram? Thank you for your time.
[[alternative HTML version deleted]]
__
R-hel
Eric,
your code work well for my need,
i'm not skillful in html environment, but your functions in R2HTML give
simple many
output requirements
thanks
Roberto
Eric Lecoutre wrote:
>
> Hi Roberto,
>
> here is a way that presumes you know some (basic) HTML tags:
>
> library(R2HTML)
> direct
There is an example in the example section of
library(zoo)
?plot.zoo
starting at:
## plot with left and right axes
On Dec 6, 2007 10:16 AM, Marcin Kozak <[EMAIL PROTECTED]> wrote:
> Hi,
>
> Consider this simple plot:
> > plot(1:25,runif(25,0,1),ylab="First Y-axis label",xaxt="n")
>
> I want to
Peter Dalgaard wrote:
>> Answering to myself and staying with the same example:
>> system.file("data/stroke.csv",package="ISwR")
>> allows direct access to the example file (name).
>>
>>
> Yes, but...
>
> This works right until you turn on LazyData for your package, then you
> end up w
huang min wrote:
>
> Hi,
>
> I have searched the lists but still can not solve the problem. I am using
> a
> windows machine. After I sweave some Rnw file, I got a tex file. However,
> the tex file can not be compiled. I know the problem is in the line
> \usepackage{C:/PROGRA~1/R/R-25~1.1/shar
An alternative solution that allows you to break it down by categories is in
the functions below. Comments/suggestions welcome and encouraged.
Note that much thanks is due to those who responded to an earlier post of
mine on a similar topic.
To use (assuming you have a data.frame falled pathDist
Hi, All;
I've compiled and installed successfully "presumably" 64-bit R on
Leopard.
But when I tried to run R, I got the following error:
/usr/bin/R: line 179: /Library/Frameworks/R.framework/Resources/etc/
i386/ldpaths: No such file or directory
Here is the configuration I used:
CONFIG_SHEL
Marcin Kozak wrote:
> Hi,
>
> Consider this simple plot:
>
>> plot(1:25,runif(25,0,1),ylab="First Y-axis label",xaxt="n")
>>
>
> I want to add an additional axis as
>
>> axis(4,at=seq(0.2,1,.2), labels=1:5)
>>
>
> I have no idea how to add now the title of the new axis as "Second
>
Hi,
Consider this simple plot:
> plot(1:25,runif(25,0,1),ylab="First Y-axis label",xaxt="n")
I want to add an additional axis as
> axis(4,at=seq(0.2,1,.2), labels=1:5)
I have no idea how to add now the title of the new axis as "Second
Y-axis label". I want this text to be vertically directed fro
Loren Engrav wrote:
>
> Thank you
>
> As per advice from several R users I have set
>
> r-project.org, stat.math.ethz.ch,fhcrc.org, stat.ethz.ch, math.ethz.ch,
> hypatia.math.ethz.ch
>
> all to be "safe domains"
>
> But still some R emails go to Junk and require to be found manually
>
> I
Johannes Graumann wrote:
> Johannes Graumann wrote:
>
>
>> On Thursday 06 December 2007 11:52:46 Peter Dalgaard wrote:
>>
>>> Johannes Graumann wrote:
>>>
Hello,
I'm trying to integrate example data in the shape of a tab delimited
ASCII file into my package and th
1 - 100 of 136 matches
Mail list logo
