Displaying 20 results from an estimated 100 matches similar to: "max/summary"
2008 Mar 04
2
paired or one-sample t-Test
Hi Guys,
I am having a real hard time trying to figure out for microarry.
Here is my code
One-Sample t-Test
dim(data.sub)
[1] 10000 140 ##there are 10000 probesets and 140 columns
hist(data.sub) ## Histogram. Identify if the probesets are normal
distributed
q<-rnorm(10000) ##generate 10000 random, normal distributed values
qqplot(data.sub,q)) ##Show the plot of the probeset
2008 Oct 15
3
request: How can we ignore a component of list having no element
Dear friends
There is a list of arrays comprising different no of rows and columns even sometimes NULL, such as [[2]] given below. How can we ignore [[2]] or others like this in the complete list. Any help in this regard is needed. Thanks
[[1]]
[,1] [,2]
[1,] 3 1
[2,] 3 1
[3,] 3 1
[[2]]
NULL
[[3]]
[,1] [,2] [,3] [,4] [,5] [,6] [,7]
[1,] 3 1
2004 Apr 22
4
RODBC installation in debian
Hello List,
I am trying to install RODBC package in a debian linux
box but getting the following message. Can anyone help
me to find what I am doing wrong here:
$ R CMD INSTALL RODBC_1.0-4.tar.gz
###
* Installing *source* package 'RODBC' ...
checking for gcc... gcc
checking for C compiler default output... a.out
checking whether the C compiler works... yes
checking whether we are cross
2009 Jun 24
1
gnls : Rho
Hello list:
How to extract the value of "Rho" from a gnls() object. I am using gnls()
function similar to
res <- gnls(y~SSmicmen(),correlation=corCompSymm(form~1|b),data=dat)
Thanks in advance,
Mahbub.
--
Mahbub Latif
School of Mathematical Sciences
Queen Mary, University of London
United Kingdom
[[alternative HTML version deleted]]
2010 Aug 11
2
help to polish plot in ggplot2
Hi,
I wanted to generate a plot which is almost like the plot generated by the
following codes.
category <- paste("Geographical Category", 1:10)
grp1 <- rnorm(10, mean=10, sd=10)
grp2 <- rnorm(10, mean=20, sd=10)
grp3 <- rnorm(10, mean=15, sd=10)
grp4 <- rnorm(10, mean=12, sd=10)
mydat <- data.frame(category,grp1,grp2,grp3,grp4)
dat.m <- melt(mydat)
p <-
2003 Sep 04
1
lattice question--- different ylim
Hi there,
I have four panels in a lattice bwplot. I want to have
two different ylim for the panels, for example panels
[1,1] and [1,2] with ylim=c(0,200) and panels [2,1]
and [2,2] with ylim=c(0,100).
Thanks for help in advance.
Mahbub.
2003 Aug 25
2
lattice question
Hi,
I want to use (similar to) las options in lattice
(bwplot) plot. Actually I want to have x-axis labels
as vertical instead of default horizontal.
Thanks in advance for your help.
Mahbub
2003 Sep 09
1
lattice plot - portrait / landscape
Hi,
How can I use portrait/landscape option in lattice
bwplot? Is there any option in trellis.device where I
can define this?
Thanks in advance,
Mahbub.
2003 Jun 09
1
installing XML in linux
Hi,
I was trying to install XML package in a linux
(dabian) machine and got the following error message.
I am not sure whether there is any error in my linux
installion. I appreciate suggestions to install XML
properly in this machine.
Thanks in advance.
Mahbub.
######Errors####
$ R CMD INSTALL -l /usr/local/lib/R/library/
XML_0.93-4.tar.gz
* Installing *source* package 'XML' ...
2003 Apr 15
5
making a dataframe out of lapply() result
Dear R-helpers,
i have a question on how to vectorize this problem:
i have a dataframe:
tester <- data.frame(groups=c("A","A","B","B","C","C"), one=c(1,1,2,2,3,3), two=c(6,6,7,7,8,8))
# i split it into a list
tester.L <- split(tester, tester$groups)
# And want to keep only the first item in each:
lapply(tester.L, function(x) x
2004 Jun 16
1
off topic: C/C++ codes for pseudo inverse
Hi,
I am looking for C/C++ codes for computing generalized
inverse of a matrix. Can anyone help me in this
regard?
Thanks,
Mahbub.
2003 Aug 21
4
anova(lme object)
Hi,
I use lme to fit models like
R> res1 <- lme(y~A+B, data=mydata, random=~1|subject)
R> res2 <- lme(y~B+A, data=mydata, random=~1|subject)
(only difference between these two models are the
sequence in which the indep variables are written in
formula)
where y is continuous and A, B, and subject are
factors. To get ANOVA table I used
R> anova(res1)
R> anova(res2)
and found
2007 Jun 22
2
multiple return
Dear User,
what's the correct way to obtain a multiple return from a function?
for example creating the simple function:
somma <- function (a, b) {
c <- a+b
return (a, b, c)
}
when I call it, it runs but returns the following output:
> somma(5, 7)
$a
[1] 5
$b
[1] 7
$c
[1] 12
Warning message:
return multi-argomento sono deprecati in: return(a, b, c)
i.e. multi-return is
2004 Oct 22
3
pgamma discontinuity (PR#7307)
Full_Name: Morten Welinder
Version: 2
OS: Solaris/space/gcc2.95.2
Submission from: (NULL) (65.213.85.217)
I changed src/nmath/standalone/test.c to read:
---------------------------------------------------------------------------------
#define MATHLIB_STANDALONE 1
#include <Rmath.h>
#include <stdio.h>
int
main()
{
double x;
for (x = 99990; x <= 100009; x++)
printf
2003 Mar 21
1
trellis plot
<html><div style=''background-color:''><P><FONT size=4>Hi there,</FONT></P>
<P><FONT size=4>I need some help about trellis plot. I have the following plot.</FONT></P>
<P><FONT size=4>x <- c(rep(LETTERS[1:4],13), rep(LETTERS[4:1],12))</FONT></P>
<P><FONT size=4>y <-
2003 Jul 22
4
greek in main title
Hello,
I have written a function that demonstrates the CLT by
generating samples following the exponential distribution,
calculating the means, plotting the histogram, and drawing
the limiting normal curve as an overlay. I have the title
of each histogram state the sample size and rate (1/theta)
for the exponential (the output is actually 4 histograms),
but I can't get the greek letter theta
2004 Nov 01
5
make apply() return a list
Hi,
I have a dataframe (say myData) and want to get a list (say myList) that
contains a matrix for each row of the dataframe myData. These matrices are
calculated based on the corresponding row of myData. Using a for()-loop to do
this is very slow. Thus, I tried to use apply(). However, afaik apply() does
only return a list if the matrices have different dimensions, while my
matrices have
2006 Apr 24
1
source compilation errors on HPUX 10.20
Hi,
I've tried to compile 3.0.22 on HPUX 10.20 and get the following
error msg: No options are used when running 'configure'.
checking configure summary
ERROR: No locking available. Running Samba would be unsafe
configure: error: summary failure. Aborting config
******
The server details:
[qhus029@HP-UX]>uname -a
HP-UX qhus029 B.10.20 A 9000/898 1684683361 two-user license
2003 Jul 21
2
bold AND italic as font in text()
Dear all,
Is it possible to somshow plot text as italic AND bold. I tried font=c(2,3)
in text(), but it doesn't work. It seems like the latter value is used.
Thanks in advance!
Sincerely,
Tord
-----------------------------------------------------------------------
Tord Sn?ll
Avd. f v?xtekologi, Evolutionsbiologiskt centrum, Uppsala universitet
Dept. of Plant Ecology, Evolutionary Biology
2004 Oct 30
2
(no subject)
Hi, there.
Does anybody know how to plot a smooth density plot for some data
simulated from certain distribution? Thanks.
Yulei
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
Yulei He
1586 Murfin Ave. Apt 37
Ann Arbor, MI 48105-3135
yuleih at umich.edu
734-647-0305(H)
734-763-0421(O)
734-763-0427(O)
734-764-8263(fax)
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$