Do you know if R has any multiv. Newton-Raphson routine implemented?
__
[EMAIL PROTECTED] mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Hi, I've been trying to find how to extract the inference info about \theta
... Is there any easy way to do it?
Thanks,
Benilton
_
Create your own personal Web page with the info you use most, at My MSN.
,
Benilton Carvalho
PhD Candidate
Department of Biostatistics
Johns Hopkins University
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide
A=matrix(1:9,3)
A[lower.tri(A)]
b
On Nov 10, 2006, at 4:50 PM, Davendra Sohal wrote:
Hi all,
Here's an interesting (for me, at least!) problem I came across:
I have a correlation matrix, let's say with 6 variables, A to F, as
column
headings and the same 6 as row headings.
The matrix is
data = read.delim(lahore.txt)
is enough for what you want to do.
b
On Nov 11, 2006, at 2:11 PM, amna khan wrote:
Respected Sir
I request you to please fill the following read.table function and
read.csvfor my understanding by assuming my data attached with this
maiL, because I
am fail to
Mea culpa. Please accept my apologies...
The problem was fixed on R-2.4.0 (final) and what I reported was on
R-2.4.0 alpha.
Good job of Martin, by the way... The problem was fixed and I didn't
notice (given that I was always using my work-around).
benilton
does
getIdx - function(tmpin, x){
tmp2 = tmpin + x
out = sort(c(tmpin, tmp2))
out = out[out=max(tmpin)]
return(out)
}
do what you want?
b
On Nov 12, 2006, at 5:20 PM, Leeds, Mark ((IED)) wrote:
I have index ( of a vector ) values of say
tempin-c(1 31 61 91 121 all the way
diff(tmp[idx])
cheers,
b
On Nov 13, 2006, at 3:06 PM, Leeds, Mark ((IED)) wrote:
I have the following set of indices, call it idx, that correspond
to the
indices of a vector say temp.
[1] 31 36 41 61 66 71 91 96 101 121 126 131 151
156 161 181 186 191 211
Hi everyone,
I have 2 environments (2 different R sessions) as described below:
Session 1:
Name of the environment: CrlmmInfo
Objects in the environment:
index1: logical index - length 238304
index2: logical index - length 238304
priors: list of 4 - (matrix 6x6, 2 vectors of length
check class()
but if all you want is to test whether it's a data.frame or not:
is.data.frame()
b
On Nov 14, 2006, at 3:07 PM, Weiwei Shi wrote:
Hi,
I am writing a generic function and need to check if an arg is a data
frame or not. I could use
is.null(dim(x)) to get what i want. But i
tb = table(df$loc, cut(df$year, seq(1970, 1985, by=5), right=F))
rs = rowSums(tb)
tb = cbind(tb, rs)
cs = colSums(tb)
tb = rbind(tb, cs)
cheers,
b
On Nov 14, 2006, at 3:20 PM, Michael Graber wrote:
Dear R List,
I am a new to R, so my question may be easy to answer for you:
I have a
How do you build your packages?
have you tried
R CMD build --binary mypkg
b
On Nov 16, 2006, at 12:32 AM, Charles Annis, P.E. wrote:
Greetings:
I've installed Rtools, MikTeX, perl, minGW, and HTML Help Workshop,
and have
succeeded in making, checking (using R CMD check mypkg) then
how about:
for (x in l) colnames(get(x)) - lower.case(colnames(get(x)))
b
On Nov 16, 2006, at 9:01 AM, Werner Wernersen wrote:
Hi,
I have a list with the names of tables, e.g.
l - c(t1,t2,t3)
and I want to change the colnames of each of the
tables in a for loop like this:
for (x in l) {
countZerosBeforeOnes - function(v)
sum(v[1:max(which(v == 1))] == 0)
apply(A, 1, countZerosBeforeOnes)
cheers,
b
On Nov 28, 2006, at 4:20 PM, Guenther, Cameron wrote:
Hi All,
If you could help me with this problem I would greatly appreciate it.
Suppose I have a matrix A:
1 1 1 1 0 1
Hi Everyone,
After searching the subject and not being successful, I was wondering
if any you could explain me the idea behind the following fact:
all(NULL == 2) ## TRUE
any(NULL == 2) ## FALSE
Thanks a lot,
Benilton
--
Benilton Carvalho
PhD Candidate
Department of Biostatistics
Johns
without breaking a lot of code, I'm sure.
Bert Gunter
Nonclinical Statistics
7-7374
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of Benilton
Carvalho
Sent: Wednesday, November 29, 2006 2:21 PM
To: R-Mailingliste
Subject: [R] tests for NULL
On Nov 29, 2006, at 7:54 PM, Martin Morgan wrote:
http://en.wikipedia.org/wiki/Empty_set
which says, in part:
Operations on the empty set
...
'any' is like 'sum', 'all' is like 'prod'.
Martin
Thank you very much Martin.
It's clear now.
Benilton
The output of sessionInfo() will be helpful here.
But I can't reproduce that on R 2.4.0 Patched r40106 + simecol (0.3-11).
Maybe you should upgrade the package?
b
On Dec 6, 2006, at 5:06 PM, Milton Cezar Ribeiro wrote:
Hi there,
I´m trying to use simecol package but I got the error
mtext(expression(beta[max]), side=1, line=2)
is it what you want?
b
On Dec 12, 2006, at 10:59 AM, javier garcia-pintado wrote:
Hi,
I'm trying to use expression() to write a text to a graphic in the
margin.
Using:
mtext(expression(beta),side=1,line=2)
writes a perfect beta greek
install.packages(plotrix)
b.
On Dec 13, 2006, at 1:42 PM, Aimin Yan wrote:
I try to type this in my R-winEdt.
but I got these. Do you know?
Aimin
install.packages('http://rh-mirror.linux.iastate.edu/CRAN/bin/
windows/contrib/2.4/plotrix_2.1-6.zip')
Warning in download.packages(pkgs,
Milton,
have you looked at the structure of your data.frame?
str(especies.aicc)
Are you sure especies.aicc is defined as numeric?
b
On Jan 9, 2007, at 10:51 AM, Milton Cezar Ribeiro wrote:
Dear Peter,
I tryed something like
head(especies.aicc)
especie aicc
1
I forgot to mention that
bitmap()
will do what you want without an X11 connection.
b
On Jan 11, 2007, at 1:09 PM, Rebecca Tagett wrote:
Hello,
I'm trying to adapt some R code that works on Windows so that it
will work
on a Linux machine.
The command :
png(myFile.png, width=600,
Hi Rebecca,
png (and also jpeg, for example) require an X11 connection.
So, assuming you're working from the command line and that your X11
server is up, you would need to do something like:
linux$ export DISPLAY=:0.0
before loading R...
and if your linux machine is remote (and you're
say your data.frame is called df
df[order(df$evidence),]
or
df[order(df$evidence, decreasing=T),] # if you want the other way
around.
b
On Jan 17, 2007, at 11:21 AM, Milton Cezar Ribeiro wrote:
Hi there,
How can I sort (order?) a data.frame using a dataframe field
(evidence) as
which.max()
b
On Jan 17, 2007, at 11:20 PM, Feng Qiu wrote:
Hi all:
A short question:
For example, a=[3,4,6,2,3], obviously the 3rd entry of
the array has the maxium value, what I want is index of the maxium
value: 3. is there a neat expression to get this index?
, D, 3, 2))
library(MASS)
myHubers - function(x)
if (length(x)1) as.numeric(hubers(x)) else c(NA, NA)
for (n in 1:N)
for (d in 1:D){
tmp - outcome[n,,d]
grp - classes[n,]
results[n, d,,] - t(sapply(split(tmp, factor(grp, levels=1:3)),
myHubers))
}
## end
--
Benilton
The null hypothesis is that the variances do not differ across
groups. The Bartlett test is sensitive to non-normality and that
might lead you to consider something more robust (eg, Levene's test).
b
On Jan 19, 2007, at 9:04 AM, Matthieu Mourroux wrote:
Bonjour,
Je voudrais tester
well, i don't use RGui, but trying:
install.packages(neural, type=source) ## from the R command line
i could not reproduce the error (which may be a problem with the copy
on the mirror you're using)...
b
On Jan 23, 2007, at 3:09 PM, Nüzhet Dalfes wrote:
Hi,
I am a total newbie to R. I am
go to R / Preferences / Startup
than you'll know what to do when you see: Default CRAN mirror.
b
On Jan 23, 2007, at 8:03 PM, Paul Gowder wrote:
Hi everyone...
Here's a good one. I'm using the R installation for mac osx. An
overly helpful friend tried to download a couple of packages for
take a look at:
?polygon
b
On Jan 25, 2007, at 8:18 AM, Mauricio Cardeal wrote:
Please, how to fill the area under the curve?
x - c(1:10)
y - c(rnorm(10))
plot(x,y)
lines(x,y)
Thanks,
Mauricio Cardeal
__
R-help@stat.math.ethz.ch mailing
In addition to Mike's comment:
x-c(0.1,0.9)
1-x[2]
[1] 0.1
x[1]==1-x[2]
[1] FALSE
all.equal(x[1], 1-x[2])
[1] TRUE
b
On Jan 26, 2007, at 11:40 AM, Mike Prager wrote:
Not at all strange, an expected property of floating-point
arithmetic and one of the most frequently asked questions
?quantile
b
On Jan 28, 2007, at 5:41 PM, Geoffrey Zhu wrote:
Hi Everyone,
I want to generate some random numbers according to some empirical
distribution. Therefore I am looking for the inverse of an empirical
cumulative distribution function. I haven't found any in R. Can anyone
give a
lattice
0.9975-11 0.9975-8 0.14-16
On Jan 29, 2007, at 7:40 AM, Michael Kubovy wrote:
On Jan 28, 2007, at 9:39 PM, Benilton Carvalho wrote:
This seems to be due to the fact that you didn't have enough
memory when running lmer2.
I might be wrong, but I think Calloc tries to get
is it sth like:
as.integer(sapply(seq(4, 22, by=9), seq, length.out=3))
you're looking for?
b
On Jan 30, 2007, at 9:29 AM, Adrian DUSA wrote:
Dear list,
This should be a simple one, I just cannot see it.
I need to generate a sequence of the form:
4 5 6 13 14 15 22 23 24
That is:
suppressWarnings(a - as.numeric(c(1, 2, pi, a, 9, z)))
b
On Jan 31, 2007, at 2:35 PM, Konrad wrote:
Hello,
Is there a way to convert a character to a number with out getting
a warning? I have a vector that has both numbers and letters in it
and I need to convert it to only numbers. At
Well, a reproducible example would be nice =)
not tested:
x = rnorm(10)
y = rnorm(20)
mymax - function(t1, t2) apply(cbind(t1, t2), 1, max)
sum(outer(x, y, mymax))
is this sth like what you need?
b
On Feb 1, 2007, at 1:18 PM, Jeffrey Racine wrote:
Greetings.
For R gurus this may be a no
The bioconductor mailing list is probably a better place to ask this
type of question.
[EMAIL PROTECTED]
But we also need to know what arrays are you working with, what the
errors are, what your sessionInfo() is
Let us know, ok?
b
On Feb 1, 2007, at 5:46 PM, Tristan Coram wrote:
Hi,
could create the linear predictor by chunks (and therefore get the
predictions). Is there a better way of doing this?
Thank you very much.
Benilton Carvalho
Department of Biostatistics
Bloomberg School of Public Health
Johns Hopkins University
__
R-help
Hi Nils,
if the server you're using is *NIX, this is what you can do:
example
ssh [EMAIL PROTECTED]
screen
R
do what you need in R
close the terminal without quitting R
ssh [EMAIL PROTECTED]
screen -r
continue working in R
%% end example
the problem is if you need X... it works
well, it's complaining because you don't have gtools installed.
how about:
install.packages(gplots, dep=T)
?
b
On Feb 20, 2007, at 1:17 PM, Randy Zelick wrote:
Hello all,
I use R on both windows and a mainframe linux installation (RedHat
enterprise 3.0, which they tell me is soon to be
well, nobody said that the density must be smaller than 1, right? :-)
it's just the value of the normal density function at the point you
asked. you may try doing that by hand and, with the correct math,
you'll get the same thing.
b
On Feb 26, 2007, at 3:03 PM, A Hailu wrote:
Hi All,
Why
is this what you mean?
tmp - combinations(3, 3, rep=TRUE)
colSums(apply(tmp, 1, duplicated))+1
b
On Mar 6, 2007, at 1:16 AM, Dylan Arena wrote:
Hi there,
I'm writing a function that calculates the probability of different
outcomes of dice rolls (e.g., the sum of the highest three rolls of
(...)
apply(tmp, 1, function(v) table(factor(v, levels=1:3)))
might be what you actually meant.
sorry for the confusion,
b
On Mar 6, 2007, at 2:00 AM, Benilton Carvalho wrote:
is this what you mean?
tmp - combinations(3, 3, rep=TRUE)
colSums(apply(tmp, 1, duplicated))+1
b
On Mar 6, 2007, at 1
say your data frame is called tmp
tmp$PRODUCTS[tmp$PRODUCTS 70] - NA
b
On Mar 26, 2007, at 12:31 PM, Sergio Della Franca wrote:
Dear R-Helpers,
I want to substitute the contents of a variable under some contitions.
I.e., I have this data set:
YEAR PRODUCTS
1 80
2
and
tmp$PRODUCTS[tmp$PRODUCTS 20] - NA.
How can i perform this double condition in the same code?
2007/3/26, Benilton Carvalho [EMAIL PROTECTED]: say your data
frame is called tmp
tmp$PRODUCTS[tmp$PRODUCTS 70] - NA
b
On Mar 26, 2007, at 12:31 PM, Sergio Della Franca wrote:
Dear R
if 'test' is your data frame...
test[, grep([tT]$, names(test))]
b
On Mar 26, 2007, at 3:06 PM, [EMAIL PROTECTED] wrote:
Dear R users
I would like to make a new dataframe from an existing dataframe,
retaining
ONLY those variables that end in the letter t
I have searched the help
have you tried replacing
.duplicate()
in your my.causality() by
vars:::.duplicate()
?
b
On Mar 28, 2007, at 3:11 PM, Leeds, Mark ((IED)) wrote:
I am using the vars package and it calls a function causality() which
then calls something
called .duplicate. I had to modify the causality
cmd = mylist = list(a = 5, b = 7)
(eval(parse(text=cmd)))
b
On Mar 28, 2007, at 11:38 AM, Brian Dolan wrote:
Hello,
I would like to take the string
mylist = list(a = 5, b = 7)
and evaluate it as a list. I have attempted to use parse and
several other functions with no success.
how about:
length(gregexpr(the, Text)[[1]])
?
b
On Apr 4, 2007, at 10:49 AM, Schmitt, Corinna wrote:
Dear R-experts,
I'm looking for an easy possibility for pattern search. I have got a
string and a special pattern. I would like to know if the pattern
appears in the string, if yes where
x=a\n
cat(x)
a
On Apr 18, 2007, at 5:31 PM, steve wrote:
If I print a sting I get an initial [1]:
xx=a
xx
[1] a
How do I get it to print just
a
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE
see
?prod
b
On May 7, 2007, at 2:25 PM, Jose Quesada wrote:
Hi,
I need to multiply all columns in a matrix so something like
apply(x,2,sum), but using multiplication should do.
I have tried apply(x,2,*)
I know this must be trivial, but I get:
Error in FUN(newX[, i], ...) : invalid unary
Hi everyone,
sorry if this was discussed before (and in this situation, could you
please point me to the discussion in the archive? My search didn't
seem to be effective).
Is there a way of getting the names of objects in a .rda file without
having to load it?
Thank you very much,
system elapsed
399.928 0.019 399.988
does not look interesting for me.
Maybe some package has some implementation of the above?
Thank you very much,
-b
--
Benilton Carvalho
PhD Candidate
Department of Biostatistics
Bloomberg School of Public Health
Johns Hopkins University
[EMAIL PROTECTED
help(toupper)
b
On May 15, 2007, at 2:41 AM, [EMAIL PROTECTED] wrote:
Dear all,
I would need a function which convert small letter into capital
letter (at
least the first letter of a character variable).
Does such a function exist in R ?
Thanks by advance
Jessica
--
Benilton
/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-
guide.html
and provide commented, minimal, self-contained, reproducible code.
--
Benilton Carvalho
PhD Candidate
Department of Biostatistics
Bloomberg School of Public Health
Johns Hopkins University
[EMAIL
install.packages(RMySQL, dep=T)
should fix it for you.
b
ps: The error says RMySQL is the problem... it is not complaining
about R itself (although it would not be a bad idea, given that the
latest R is v 2.5.0, so it would be a better idea to start by
upgrading your R)
On May 22, 2007,
AFAIK, it only depends on how much free memory you have.
b
On Jun 1, 2007, at 5:05 PM, Guanrao Chen wrote:
hi, Rers
I tried to find out the max size (# of rows, # of
columns) of a matrix that is allowed by R but failed.
Can anybody let me know?
Thanks!
Guanrao
iconv on your linux box should do the work.
b
On Jun 1, 2007, at 7:19 PM, Alberto Vieira Ferreira Monteiro wrote:
I have a windows source file.r, with the default charset of windows.
I can't use it in Linux as source(file.r), because Linux's
default is
Unicode. How can I read it?
Alberto
Well, I have my email client to organize everything by thread...
which does the work for me...
ethz: Eidgenössische Technische Hochschule Zürich
ch: Confœderatio Helvetica
best,
b
On Jun 4, 2007, at 7:25 PM, Robert Wilkins wrote:
Why does the R mailing list need such an unusual and
)
was used insted, I'd like to get NA.
thank you very much,
b
--
Benilton Carvalho
PhD Candidate
Department of Biostatistics
Bloomberg School of Public Health
Johns Hopkins University
[EMAIL PROTECTED]
__
R-help@stat.math.ethz.ch mailing list
https
sorry, I hit send before finishing my thoughts...
and as for clustering microarray data, you might want to consider the
bioconductor mailing list...
[EMAIL PROTECTED]
b
On Jun 7, 2007, at 10:42 PM, ssls sddd wrote:
Dear List,
I have another question to bother you about how to do
Hi Alex,
just in case you're trying to get genotypes from the Affymetrix 500K
set, you might want to check the oligo package available on
BioConductor.
best,
b
On Jun 7, 2007, at 10:42 PM, ssls sddd wrote:
Dear List,
I have another question to bother you about how to do clustering.
My
which(a == .4)[1]
b
On Jun 10, 2007, at 4:45 AM, gallon li wrote:
find the position of the first value who equals certain number in a
vector:
Say a=c(0,0,0,0,0.2, 0.2, 0.4,0.4,0.5)
i wish to return the index value in a for which the value in the
vector is
equal to 0.4 for the first
the recommendation was to use lower.tail=FALSE.
b
On Jun 11, 2007, at 11:21 AM, Carmen Meier wrote:
I got an answer for the other question (thank you)
But there is another question (I am afraid this is a basic
question ...)
In this tread there is a hint hwo to calculate the p-vlue of
get 1) and z=1.96 (you should get 5%)
b
On Jun 11, 2007, at 11:34 AM, Carmen Meier wrote:
Benilton Carvalho schrieb:
the recommendation was to use lower.tail=FALSE.
b
O
but then the results are significant and this does not match the
observation.
The results are matching the observations
date = c(1, 1, 1, 1, 2, 2, 3, 3, 4)
tag = c(1, 2, 3, 4, 2, 4, 1, 2, 4)
table(factor(tag, levels=1:4), factor(date, levels=1:4))
(not sure how you got Tag 1/Date 4 = 1)
On Jun 11, 2007, at 3:11 PM, [EMAIL PROTECTED] wrote:
What code can i use to convert a table like this:
Tag#Date
1
Hi,
First of all, I apologize for sending out this message, as I'm sure
the answer is on the archives, but I just can't find it.
Not long ago, there was a discussion about building Windows packages
from the source code and someone posted a link to a website to which
we could submit the
Hi everyone,
it's been a while I've been trying to save a plot created via
lattice:::xyplot
if I have a file tst.R with the following code:
y - rnorm(100)
x - rnorm(100)
z - sample(letters[1:4], 100, rep=T)
library(lattice)
bitmap(tst.png)
xyplot(y~x|z)
dev.off()
and I source it, I get the
()
}
test(100)
## end test.R
source(test.R, echo=T)
also fails in this case...
thanks a lot,
b
On Jun 15, 2007, at 8:53 PM, [EMAIL PROTECTED] wrote:
On 6/15/07, Benilton Carvalho [EMAIL PROTECTED] wrote:
Hi everyone,
it's been a while I've been trying to save a plot created via
lattice
replaced xyplot()
by plot().
Thank you very much,
b
On Jun 16, 2007, at 12:26 AM, [EMAIL PROTECTED] wrote:
On 6/15/07, Benilton Carvalho [EMAIL PROTECTED] wrote:
So, if those statements are inside a function, I have to make my
function to have an 'echo' argument/functionality? eg.:
## begin
?dir.create
b
On Jun 20, 2007, at 8:07 AM, Milton Cezar Ribeiro wrote:
Hi all,
How can I create (and check the existence of) a directory in a R
session under Windows(xp)?
Kind regards,
Miltinho
__
R-help@stat.math.ethz.ch mailing list
outer(test, fac, )
-b
On Jun 26, 2007, at 2:13 PM, Van Campenhout Bjorn wrote:
hi all, sorry for this basic question, I think I know I should use ?
apply, but it is really confusing me...
I want to create a matrix by comparing two vectors. Eg:
test-seq(1:10)
fac-c(3,6,9)
and i want
or just type:
pairs.default
b
On Jun 28, 2007, at 2:36 PM, Henrique Dallazuanna wrote:
https://svn.r-project.org/R/branches/R-2-5-branch/src/library/
graphics/R/pairs.R
--
Henrique Dallazuanna
Curitiba-Paraná-Brasil
25° 25' 40 S 49° 16' 22 O
On 28/06/07, Jiong Zhang, PhD [EMAIL
qqnorm(table[,1])
is what you want, isn't it?
and other forms would include:
par(ask=TRUE)
results = apply(table, 2, qqnorm)
par(ask=FALSE)
b
On Jun 28, 2007, at 9:50 PM, Jiong Zhang, PhD wrote:
I want to qqnorm every column in a table. When I try the first column
using
Hi,
Until now, I thought that the results of glm() and bigglm() would
coincide. Probably a naive assumption?
Anyways, I've been using bigglm() on some datasets I have available.
One of the sets has 15M observations.
I have 3 continuous predictors (A, B, C) and a binary outcome (Y).
And
[1] -93196.69
Thank you very much for the suggestion, I'll give it a try.
Best,
benilton
On Jun 29, 2007, at 11:05 AM, Peter Dalgaard wrote:
Benilton Carvalho wrote:
Hi,
Until now, I thought that the results of glm() and bigglm() would
coincide. Probably a naive assumption?
Anyways, I've
C1 - rep(-1, length(Cat))
C1[Cat == b]] - 1
b
On Jul 4, 2007, at 9:44 AM, Keith Alan Chamberlain wrote:
Dear Rhelpers,
Is there a faster way than below to set a vector based on values from
another vector? I'd like to call a pre-existing function for this,
but one
which can also handle
the smarter thing is to write a package, so you don't need source()
at all.
but your problem will be fixed if you use the full path instead:
source(/home/user/dir1/dir2/file2.R)
(and obviously you could pass the path as an argument to whatever
function you're using...)
b
On Jul 16, 2007,
The R Data Import/Export also says that this function is in the
foreign package. :-)
b
On Jul 16, 2007, at 11:38 PM, Brad Christoffersen wrote:
To Whom It May Concern:
I want to read in an S-PLUS data dump and I used
data.restore(filepath/filename)
(in R 2.5.1 for Windows 95 and
maybe:
x = c(.2, .1, .8, .3, .7, .6, .01, .2, .5, 1, 1)
breaks = seq(0, 1, .2)
LETTERS[1:(length(breaks)-1)][cut(x, breaks)]
b
On Jul 18, 2007, at 1:50 PM, Doran, Harold wrote:
Michael
Assume your data frame is called data and your variable is called
V1. Converting this to a factor is:
Check the following example for by():
require(stats)
attach(warpbreaks)
by(warpbreaks, tension, function(x) lm(breaks ~ wool, data = x))
or just type:
example(by)
b
On Jul 19, 2007, at 1:18 PM, Hongmei Jia wrote:
Dear All,
I'm trying to do 'glm' analysis by groups just
it looks like that whatever method you used to genotype the 1002
samples on the STY array gave you a transposed matrix of genotype
calls. :-)
i'd use:
genoT = read.table(yourFile, stringsAsFactors = FALSE)
as a starting point... but I don't think that would be efficient (as
you'd need to
))
##
##
user system elapsed
119.288 0.004 119.339
(for all 240K)
best,
b
ps: note that out is a list.
On Jul 20, 2007, at 2:01 AM, Latchezar Dimitrov wrote:
Hi,
-Original Message-
From: Benilton Carvalho [mailto:[EMAIL PROTECTED]
Sent: Friday, July 20, 2007 12:25 AM
jpeg(...) ##if you have X11
bitmap(..., type=jpeg) ##otherwise
b
On Jul 20, 2007, at 11:34 AM, Ding, Rebecca wrote:
Dear R users,
I used R to draw many histograms and I would like to automatically
save
them into a jpeg file. I tried the following code since I know .ps
file
could be
set.seed(123)
N = 3
K = 400
theData = matrix(rnorm(N*K), ncol=K)
theData = as.data.frame(theData)
theData = cbind(indicator = sample(0:1, N, rep=T), theData)
system.time(results - colMeans(subset(theData, indicator == 1)))
user system elapsed
2.309 1.319 3.853
b
On Jul 20,
as.integer(factor(dta[[school_id]]))
b
On Jul 20, 2007, at 9:26 PM, [EMAIL PROTECTED] wrote:
Hi All
I was wondering if I can avoid a time-consuming for loop on my
60 obs dataset.
school_id y
8 9.87
8 8.89
8 7.89
8 8.88
20 6.78
test - lapply(1:3, function(i) cbind(runif(15), rnorm(15,2)))
sapply(test, [, 16:30)
b
On Jul 22, 2007, at 3:39 PM, [EMAIL PROTECTED] wrote:
Hello,
I think I have a mental block when it comes to working with lists.
lapply and sapply appear to do some magical things, but I can't
seem
oh! and if you want to be less ad-hoc:
sapply(test, function(x) x[,2])
b
On Jul 22, 2007, at 3:50 PM, Benilton Carvalho wrote:
test - lapply(1:3, function(i) cbind(runif(15), rnorm(15,2)))
sapply(test, [, 16:30)
b
On Jul 22, 2007, at 3:39 PM, [EMAIL PROTECTED] wrote:
Hello,
I think I
are you positive that your function is doing what you expect it to do?
it looks like you want something like:
sapply(1:10, function(i) cor(mat1[i,], mat2[i,]))
b
On Jul 24, 2007, at 11:05 AM, Bernzweig, Bruce ((Consultant)) wrote:
Hi,
I've created the following two matrices (mat1 and mat2)
is calculate correlations on a row against row
basis: mat1 row1 x mat2 row1, mat1 row1 x mat2 row2, ... mat1 row1 x
mat2 row-n, mat1 row-n, mat2 row-n
- Bruce
-Original Message-
From: Benilton Carvalho [mailto:[EMAIL PROTECTED]
Sent: Tuesday, July 24, 2007 11:31 AM
To: Bernzweig
after dat.cor use:
Rmat[lower.tri(Rmat)] - dat.cor
Rmat - t(Rmat)
Rmat[lower.tri(Rmat)] - dat.cor
b
On Jul 27, 2007, at 11:28 PM, Gregory Gentlemen wrote:
Greetings,
I have a seemingly simple task which I have not been able to solve
today. I want to construct a symmetric matrix of
a - matrix(1:6, nr=1)
colnames(a) - paste(col, 1:6)
xtable(a)
On Jul 28, 2007, at 12:39 PM, Stefan Nachtnebel wrote:
Hello,
Is there a possibility to use xtable with a vector to generate a latex
table? I always get an error, that no applicable method is available.
For example:
b-1:12
maybe this is what you want?
plot(rnorm(10))
legend(topleft, A), bty=n)
?
b
On Aug 7, 2007, at 11:08 AM, Daniel Brewer wrote:
Simple question how can you position text in the top left hand
corner of
a plot? I am plotting multiple plots using par(mfrow=c(2,3)) and
all I
want to do is
As long as you keep in mind Prof. Ripley's comment, you're going to
be fine with nchar().
http://tolstoy.newcastle.edu.au/R/e2/devel/07/05/3450.html
Remember that what you want exactly is given by nchar(obj,
type=chars), which is **NOT** the default on R 2.5.1 (only on
R-2.6.0).
In your
observation %in% ID
b
On Sep 7, 2007, at 1:40 AM, Takatsugu Kobayashi wrote:
Hi RUsers,
I am wonder if I can search observations whose IDs matches any of the
values in another vector, such as in MySQL. While I am learing
MySQL for
future database management, I appreciate if anyone could
95 matches
Mail list logo