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
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
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
> wan
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:
>
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
ear.
>
> What I am trying to do 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 PRO
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 mat
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, [EMA
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
>
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
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,
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
> co
me(out <- lapply(genoT, function(x) match(x, c("AA", "AB",
"BB"))-1))
##
##
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,
>
>&
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 f
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 ju
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
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 9
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,
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
2$deviance
[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
>
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 tr
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
>
> qqnorm(table$
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 <
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)
>
>
?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 lis
just fine if I 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
t;tst.png")
xyplot(y~x|z)
dev.off()
}
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
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 t
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 sou
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
ou should 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 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
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
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.
>
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 cluste
e string "x".
if,
foo(arg1)
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]
__
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 customiz
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?
>
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
_
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,
e mail please)
>
> __
> 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 commented, minimal, self-con
R ?
>
> Thanks by advance
>
> Jessica
--
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://stat.ethz.ch/m
oup(x[i,], g
[i,])))
user 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 U
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,
benilto
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], ...) : inval
> 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/li
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
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 suc
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 fu
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
p$PRODUCTS[tmp$PRODUCTS > 70] <- NA
>
> 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&
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
>
arrangement (...)"
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,
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 rol
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,
> Wh
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
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
thought is that saving the Boundary.knots and knots, I
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 Univer
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,
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 n
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 numb
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
uot;datasets"
"methods"
[7] "base"
other attached packages:
lme4 Matrix 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:
?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
> giv
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 quest
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 m
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 package
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.
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 l'homo
t(sapply(split(tmp, factor(grp, levels=1:3)),
myHubers))
}
## end
--
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
PLEAS
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?
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) a
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 connec
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", widt
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 Atti
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.packa
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 gr
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 s
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
_
d conventions, but
> even if
> you disagree, they are certainly not worth making a fuss over and
> certainly
> cannot be changed without breaking a lot of code, I'm sure.
>
> Bert Gunter
> Nonclinical Statistics
> 7-7374
>
> -Original Message-
> Fro
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
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
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:
>
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
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 da
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
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
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 21
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 th
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
__
On Nov 12, 2006, at 3:35 PM, Peter Dalgaard wrote:
>> My fix was to remap my group variable so that it would always start
>> from 1.
>
> Ugh!
>
> You shouldn't "fix" something like that, instead try your very best to
> find circumstances where you reproduce the crash as quickly as
> possible, pre
I had a similar problem, but with the silhouette() function.
I figured that the function expected that the grouping variable
started from 1, ie., if for some reason on my obs I didn't have any
observation on group 1, the function would return me random numbers
and repeated execution of the c
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 fa
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 matr
lution did not work for me.
Any suggestion?
Thank you very much,
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 post
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.
http://click.atdmt.com/AVE
97 matches
Mail list logo