Gabor,
thanks for the clarification, now I understand the expression.
Thanks to everyone
Bart
>From: "Gabor Grothendieck" <[EMAIL PROTECTED]>
>To: "Bart Joosen" <[EMAIL PROTECTED]>
>CC: r-help@stat.math.ethz.ch
>Subject: Re: [R] How to read in this data format?
>Date: Thu, 1 Mar 2007 16:46:21
Dear R list:
I have a data like this
> head(data.19.pr.2)
prAveSd M#Re Aa
1 1A24 57.66 33.50 20 ALA_1 ALA
2 1A24 72.16 19.75 20 GLN_2 GLN
3 1A24 103.52 8.64 20 TYR_3 TYR
4 1A24 38.67 15.51 20 GLU_4 GLU
5 1A24 54.56 16.43 20 ASP_5 ASP
6 1A24 999.00 0.00 20 GLY_6 GLY
> leve
On Mar 1, 2007, at 6:28 PM, <[EMAIL PROTECTED]> wrote:
First off, if you are working in perl you might want to be aware of
ruby and the "r for ruby" project:
http://rubyforge.org/projects/r4ruby/
> Hello,
>
> $R->send(qq (xVal <- c(1,2,3,4,5,6)));
> $R->send(qq (yVal <- c(3,5,2,6,1,5)));
> $R->
Hi Bert,
On Thursday 01 March 2007 23:59, Bert Gunter wrote:
> I am not sure exactly what you are asking for below,
I guess I'm really asking for some advanced comfort / elegance. See below.
> but I wonder if your
> query could be satisfied by the judicious use of the ... argument in a
> wrappe
On Thu, 2007-03-01 at 15:33 -0800, Bricklemyer, Ross S wrote:
> All,
>
> Apparently the attachment did not make it through the list serve.
>
> I am having difficulty installing R-2.4.1 on Mandriva Linux 2007
> Discovery. Thanks to the help of Doug Bates I got further along the
> path (Thanks Dou
On Mar 1, 2007, at 4:36 PM, Robert Duval wrote:
> You might want to start looking at the FAQ's
>
> http://cran.r-project.org/faqs.html
>
> in particular
>
> http://cran.r-project.org/doc/FAQ/R-FAQ.html#R-and-S
>
Thanks I must admit that I had not looked at the FAQ's, but I have
now and though it
All,
Apparently the attachment did not make it through the list serve.
I am having difficulty installing R-2.4.1 on Mandriva Linux 2007
Discovery. Thanks to the help of Doug Bates I got further along the
path (Thanks Doug!!). I now cannot figure out which line of code needs
modified in order fo
On 01-Mar-07 Paul Lynch wrote:
> I'm not sure this is the most efficient, but how about:
>diag(m[a,b])
> ?
m[cbind(a,b)] will also do it:
m
[,1] [,2] [,3] [,4] [,5] [,6] [,7]
[1,] 1.1 1.2 1.3 1.4 1.5 1.6 1.7
[2,] 2.1 2.2 2.3 2.4 2.5 2.6 2.7
[3,] 3.1 3.2 3.3 3.4 3.5 3
On Windows XP it worked for me on both 2.4.1 and 2.5.0. I did notice
that on 2.4.1 it says "using Synchronous WinInet calls" but does not
say this on 2.5.0. See below for the two transcripts.
> ftp <- "ftp://anonymous:[EMAIL PROTECTED]/edgar/full-index/company.idx"
> download.file(url=ftp, destf
Hello,
I am in the process of writing a Perl program to carry out and analyze
large numbers of regressions using R as the regression engine, and I am
using Statistics::R to create the communication bridge between the two
programs. Statistics::R creates a pipe between R and Perl and uses
Rterm.exe
try:
> a = c(1, 4, 5)
> b = c(2, 6, 7)
>
> m <- matrix(1:49,7,7)
> m[cbind(a,b)] # 'array' indexing
[1] 8 39 47
> m
[,1] [,2] [,3] [,4] [,5] [,6] [,7]
[1,]18 15 22 29 36 43
[2,]29 16 23 30 37 44
[3,]3 10 17 24 31 38 45
[4,]4 11 18
On Thu, 2007-03-01 at 12:35 -0800, lalitha viswanath wrote:
> Hi
> Thanks much for that input.It was extremely helpful.
>
> I am seeking some input about another stumbling block
> using RODBC; SQLQuery et al with large BLOB values.
>
> Although the following query
> dataFromDB <- sqlQuery(channe
Thomas:
I am not sure exactly what you are asking for below, but I wonder if your
query could be satisfied by the judicious use of the ... argument in a
wrapper function to par(), like
myPar=function(bg="lightgray", pch=19,...)par(bg=bg,pch=pch,...)
or perhaps
myX11 <- function(width=10, bg="l
All,
I have managed to download files from web sites and ftp sites using R,
so just for fun I tried to do so from the SEC's ftp site using the
following code:
ftp <- "ftp://anonymous:[EMAIL PROTECTED]/edgar/full-index/company.idx"
download.file(url=ftp, destfile="test.txt")
And it does not work.
I'm not sure this is the most efficient, but how about:
diag(m[a,b])
?
On 3/1/07, yoo <[EMAIL PROTECTED]> wrote:
>
> Let's say i have
>
> a = c(1, 4, 5)
> b = c(2, 6, 7)
>
> and i have matrix m, what's an efficient way of access
> m[1, 2], m[4, 6], m[5, 7]
> like of course m[a, b] = is not
I would like to use R for submissions to FDA/CDRH (the medical device
company I work for currently uses only SAS). Previous postings to the list
regarding R and 21 CFR 11 compliance have been very helpful. However,
reluctance to using open source software for statistical analyses and
reporting r
Hi,
I am writing a function and I want to pass arguments to a function I
will call internally without having to specify all the possible
arguments to pass. Usually I would use '...' but I want to do this for
two functions that I will call, and the two functions do not take the
same arguments. S
This is a covariance calculation question so nothing to do with R but
maybe someone could help me anyway.
Suppose, I have two random variables X and Y whose means are both known
to be zero and I want to get an estimate of their covariance.
I have n sample pairs
(X1,Y1)
(X2,Y2)
.
.
.
.
.
(Xn,Yn)
Hi Petr,
On Thursday 01 March 2007 21:50, Petr Klasterecky wrote:
> I am no expert on these topics but currently I am solving a similar
> issue using the .Rprofile file and the .First function. So maybe it's
> enough to put
> .First <- function(){
> par(whatever you want)
> further instructions if
>> Also, if x is a data frame, latex(x) contains the row numbers.
>> Can I get rid of them here as well?
I think you are asking for the rowname=NULL argument.
latex(x, rowname=NULL)
See ?latex to confirm if that is what you are looking for.
Rich
__
Dear list,
I have this huge array of numbers, of dimensions 67 x 33 x 51 x 6, the 6 being
the replications. I wanted to test for evidence of autocorrelation between the
6 replications, marginally. I can calculate the first-order autocorrelation
very easily using an appropriately defined functi
You might want to start looking at the FAQ's
http://cran.r-project.org/faqs.html
in particular
http://cran.r-project.org/doc/FAQ/R-FAQ.html#R-and-S
robert
On 3/1/07, Charilaos Skiadas <[EMAIL PROTECTED]> wrote:
> I just acquired a copy of "Statistical Models in S", I guess most
> commonly know
Let's say i have
a = c(1, 4, 5)
b = c(2, 6, 7)
and i have matrix m, what's an efficient way of access
m[1, 2], m[4, 6], m[5, 7]
like of course m[a, b] = is not going to do, but what's an expression that
will allow me to have that list?
Thanks!
--
View this message in context:
http://www.nabb
The White Book provides the original S Language Specification. This was what
existed at Bell labs way back then. Subsequent implementations, both S-Plus
and R, will differ on details.
Also, a lot of development effort has flowed over the dam since publication,
so both implementations contain lots
On 3/1/07, Bart Joosen <[EMAIL PROTECTED]> wrote:
> Dear All,
>
> thanks for the replies, Jim Holtman has given a solution which fits my
> needs, but Gabor Grothendieck did the same thing,
> but it looks like the coding will allow faster processing (should check this
> out tomorrow on a big datafil
Hi.
Here's an odd request that actually arose out of my own bad planning.
Suppose I do the following:
which(v1 %in% v2)
I will get a set of indices for v1, and they will be ordered in the same
order that v1 is in. I want the indices of the intersection for v1 ordered
according to v2.
I do t
I am no expert on these topics but currently I am solving a similar
issue using the .Rprofile file and the .First function. So maybe it's
enough to put
.First <- function(){
par(whatever you want)
further instructions if neccessary
}
Petr
Thomas Friedrichsmeier napsal(a):
> The following questi
Ranjan Maitra napsal(a):
> Here's one suggestion: convert the matrix into a three-dimensional array and
> use apply on it.
Converting to 3 dims should not be neccessary:
m <- matrix(rnorm(110),ncol=22)
t.list <- apply(m,1,function(x){t.test(x[1:11],x[12:22],paired=TRUE)})
However, I have no ide
I just acquired a copy of "Statistical Models in S", I guess most
commonly known as the "white book", and realized to my dismay that
most of the code is not directly executable in R, and I was wondering
if there was a source discussing the things that are different and
what the new ways of
Hi
Thanks much for that input.It was extremely helpful.
I am seeking some input about another stumbling block
using RODBC; SQLQuery et al with large BLOB values.
Although the following query
dataFromDB <- sqlQuery(channel, "select
uncompress(columnName) from tableName where Id=id ");
returns jus
Dear All,
thanks for the replies, Jim Holtman has given a solution which fits my
needs, but Gabor Grothendieck did the same thing,
but it looks like the coding will allow faster processing (should check this
out tomorrow on a big datafile).
@gabor: I don't understand the use of the grep command
Sorry, meant for this to go to the whole list.
On Mar 1, 2007, at 3:29 PM, Charilaos Skiadas wrote:
> On Mar 1, 2007, at 2:52 PM, steve wrote:
>
>> Unfortunately, this applies to print.xtable, and not to latex. I
>> want to
>> know how to eliminate them using latex()
>
> 1) Why do you need to u
Use 'strsplit' to create a vector and then 'rbind' into a matrix or create a
dataframe:
> x <- c("id1\tid2\tid3\tid4", "1\t2\t3\t4", "5\t6\t7\t8")
> x
[1] "id1\tid2\tid3\tid4" "1\t2\t3\t4" "5\t6\t7\t8"
> (y <- strsplit(x, "\t"))
[[1]]
[1] "id1" "id2" "id3" "id4"
[[2]]
[1] "1" "2"
A new version of X2R has just been uploaded to CRAN. It should be
available at mirrors within a few days.
This contains revisions to the For2R component to fix a bug in which
data frames were not written correctly if the user did not pass row labels.
The new version is supplied as files X2R.zi
> set.seed(1)
> student.cop <- ellipCopula("t", param = c(0.5, 0.6, 0.7), dim = 3, dispstr =
> "un",df=5)
> x<-rcopula(student.cop,1000)
> fit <- fitCopula(x, student.cop, c(0.5,0.5,0.5,5))
Warning messages:
1: NaNs produced in: qt(p, df, lower.tail, log.p)
2: NaNs produced in: sqrt((df + d)/(Q +
Unfortunately, this applies to print.xtable, and not to latex. I want to
know how to eliminate them using latex()
Greg Johnson wrote:
> steve bowdoin.edu> writes:
>
>> I would like to get rid of the row numbers using xtable and latex.
>>
>> but I don't want the row numbers. Is it possible to g
The following question/idea came up on the RKWard development mailing list,
but might be of general interest:
Is there a nice way to customize the default look of all graphs on all
devices? I.e. a way to - for instance - set the following options before each
plot:
par(bg="light gray", las=2, p
You can split() the original vector according to its elements into a list,
then use lapply() to count the lengths of the list:
> x <- round(runif(100,1,10))
> unlist(lapply(split(x,f=x),length))
1 2 3 4 5 6 7 8 9 10
5 14 9 13 11 8 8 14 12 6
This can be useful if you want to use so
You didn't' actually tell us what forfact is. You could try something like
(untested)
My.matrix <- my.matrix[!(my.matrix[, 1] %in% remove.values), , drop = FALSE]
Also, reading some of the introductory R manuals you can find on CRAN in the
documentation section will really give you a leg up on
I'd suggest to use the function mt.teststat in the package multtest or
rowttests in the package genefilter. Both can be found athe the bioconductor
webpage (www.bioconductor.org)
Stef
On Thu, Mar 01, 2007 at 12:48:29PM -0600, Ranjan Maitra wrote:
Here's one suggestion: convert the matrix i
I decided to run an experiment: just reading in a file which is 78MB in binary
format (of ints). It takes less than 30s using a laptop with 512 MB RAM, 2.3
GHz Intel-4 single processor. At that point, I did not notice that Ramzi was
talking about a .RData file.
For huge files, I usually do not
Dear all,
again i have a problem with special case of dropping lines out of a
matrix.
i tried the following, where throwout is a vector of length 18 with
simple id values that should be compared to any
for (k in 1:length(throwout))
{
mymatrix=matrix(mymatrix[-(forfact[k]),],ncol=4)
On Thu, 2007-03-01 at 09:22 -0800, ramzi abboud wrote:
> Is R file IO slow in general or am I missing
> something? It takes me 5 minutes to do a load(MYFILE)
> where MYFILE is a 27 MB Rdata file. Is there any way
> to speed this up?
>
> The one idea I have is having R call a C or Perl
> routin
Daniel Jimenez wrote:
> Dear R users I'm new in R management and maybe It's a silly question. I'm
> working with GLM to obtain predictive models. I have some problesm
> with the
> prediction instruction:
>
>
>
>> DatosTotal <- read.csv("Var_perdizcsv.csv", sep =";")
>> edvariable <- edit(DatosTota
Here's one suggestion: convert the matrix into a three-dimensional array and
use apply on it.
Ranjan
On Thu, 1 Mar 2007 11:51:29 -0600 (CST) Nameeta Lobo <[EMAIL PROTECTED]> wrote:
> Hello all,
>
> I am trying to run a two sample t-test on a matrix which is a
> 196002*22 matrix. I want to ru
Here is one way of doing it using the reshape package:
> # test data from email
> x <- "$$ Experiment Number:
+ $$ Associated Data:
+
+ FUNCTION 1
+
+ Scan 1
+ Retention Time 0.017
+
+ 399.8112 184
+ 399.8742 0
+ 399.9372 152
+
+
+ Scan 2
+ Retention Time 0.021
+
+ 399.8112 181
+ 399.8
On Thu, 2007-03-01 at 18:20 +0100, Bruno C. wrote:
> Hello,
> I have several questions around matrices size limitation:
> can i create a matrix with 70 rows?
> does it depends on the number of columns?
> if so can Is there a way to ask to R, given the number of
> columns, the max
Hi
Mark Allen wrote:
> Dear Reader
>
> I am trying to change the font in a plot and after several trials finally
> came up with the following code.
>
> plot(var_a, var_b, pch = 16, font.lab = 10, font = 10)
> points(var_a, var_c, pch = 3, font = 10)
> legend(0.1, 0.8, legend = c(“var_b”, “var_
On Thu, 2007-03-01 at 23:48 +1000, Mark Allen wrote:
> Dear Reader
>
> I am trying to change the font in a plot and after several trials finally
> came up with the following code.
>
> plot(var_a, var_b, pch = 16, font.lab = 10, font = 10)
> points(var_a, var_c, pch = 3, font = 10)
> legend(0.1,
It is not slow on my system. The file was 34MB on disk and took about 37
seconds to write out (probably mostly disk I/O on my laptop) and 12 seconds to
read in after I flushed the system cache.
> x <- runif(27e6/4) # creates a 34MB file on disk
> object.size(x)
[1] 5424
> system.time(save.
I have temp list which have 19 data.frame
I want to get number of levels for pr in the first dat.frame
I do this like this:
temp[[1]]$pr just has "1A24"
after I do nlevels(temp[[1]]$pr)
I expect to get 1, but I get 19
anyone know why?
> tail(temp[[1]]$pr)
[1] 1A24 1A24 1A24 1A24 1A24 1A24
19 Lev
Hi Radek,
Radek Blatny <[EMAIL PROTECTED]> writes:
> Dear colleagues,
> I've just come across a problem with the following command which is a
> part of the "metaOverview.R" code file provided as an monography-
> accompanying file at
> http://www.bioconductor.org/docs/mogr/metadata:
The approp
A 27MB .RData file is relatively big, in may experience. What do you think is
slow? Maybe it's your computer that is slow?
-roger
ramzi abboud wrote:
> Is R file IO slow in general or am I missing
> something? It takes me 5 minutes to do a load(MYFILE)
> where MYFILE is a 27 MB Rdata file. I
Dear R users I'm new in R management and maybe It's a silly question. I'm
working with GLM to obtain predictive models. I have some problesm with the
prediction instruction:
DatosTotal <- read.csv("Var_perdizcsv.csv", sep =";")
edvariable <- edit(DatosTotal)
pre <- predict(rlfinal, DatosTotal,
Just an idea: Two things that can slow down save()/load() is if you
save() in ASCII format or a compressed binary format. If this is your
case for MYFILE, try to resave in a non-compressed binary format. See
?save for details.
/HB
On 3/1/07, ramzi abboud <[EMAIL PROTECTED]> wrote:
> Is R file I
Read in the data using readLines, extract out
all desired lines (namely those containing only
numbers, dots and spaces or those with the
word Time) and remove Retention from all
lines so that all remaining lines have two
fields. Now that we have desired lines
and all lines have two fields read the
I agree with Jason: !2, prefer 1, can accept 3.
Clint BowmanINTERNET: [EMAIL PROTECTED]
Air Dispersion Modeler INTERNET: [EMAIL PROTECTED]
Air Quality Program VOICE: (360) 407-6815
Department of Ecology FAX:(360)
Is R file IO slow in general or am I missing
something? It takes me 5 minutes to do a load(MYFILE)
where MYFILE is a 27 MB Rdata file. Is there any way
to speed this up?
The one idea I have is having R call a C or Perl
routine, reading the file in that language, converting
the data in to R obj
Hello all,
I am trying to run a two sample t-test on a matrix which is a
196002*22 matrix. I want to run the t-test, row-wise, with the
first 11 columns being a part of the first group and columns
12-22 being a part of the second group.
I tried running something like (temp.matrix being my 196002
Hello,
I have several questions around matrices size limitation:
can i create a matrix with 70 rows?
does it depends on the number of columns?
if so can Is there a way to ask to R, given the number of columns, the max
number of rows I can have?
I need in fact to do regression p
Hello
I am looking for ROracle - Module but found only a version that can
installed with R Version > 2.30.
I am running a production system where I can't change the currently
installed R. That's why I need a ROracle - Module that is compatible
with R version 2.2.1
Thanks
Reinhard
__
Dear Reader
I am trying to change the font in a plot and after several trials finally
came up with the following code.
plot(var_a, var_b, pch = 16, font.lab = 10, font = 10)
points(var_a, var_c, pch = 3, font = 10)
legend(0.1, 0.8, legend = c(var_b, var_c), pch = c(16,3))
It does change t
Definitely not #2. Prefer #1 but #3 is ok as well.
Thanks for contributing and inquiring.
- Original Message -
From: "hadley wickham" <[EMAIL PROTECTED]>
To: <[EMAIL PROTECTED]>
Sent: Sunday, February 25, 2007 7:44 AM
Subject: [R] Double-banger function names: preferences and suggesti
On Thu, 2007-03-01 at 08:34 -0800, lalitha viswanath wrote:
> Hi
> Thanks much for the prompt response to my earlier
> enquiry on packages for regression analyses.
> Along the same topic(?), I have another question about
> which I could use some input.
>
> I am retreiving data from a MySQL databas
I think it will be better if I explained what I want to do.
I have the abundance of Penilia (zooplankter) as y and I have the
temperature, salinity and chlorophyll concentration at different depths. I
want to know which how the abundance of Penilia is modulate by those
environmetal variables and w
Hello everybody,
I have a big problem that I do not manage to solve !
I will be very grateful if you can solve this !
I want to fit a t Copula with the copula package :
> student.cop <- ellipCopula("t", param = c(0.5, 0.6, 0.7), dim = 3, dispstr =
> "un",df=5)
> x<-rcopula(student.cop,1000)
> fi
Hi
Thanks much for the prompt response to my earlier
enquiry on packages for regression analyses.
Along the same topic(?), I have another question about
which I could use some input.
I am retreiving data from a MySQL database using
RODBC.
The table has many BLOB columns and each BLOB column
has d
You can also have a look at structable() in vcd, especially the indexing
functions. The problem with OLAP in R is, that you will have to create
sth. like a hierarchical factor to handle rollup/drilldown correctly.
And you will have to decide whether it's purely memory-based (fast
calculations,
Dear R users,
I am having a little trouble with grouping data.
---Detailed explanation (summary below)
A small sample of my data is below (which has already been rounded and
grouped a little from the raw data for clarity).
I am sampling data from an unknown gam
You can't expect general-purpose tools like read.table in R to be able
to deal with highly specialized file format. Here's what I'd start. It
doesn't put data in the format you specified exactly, but I doubt you'll
need that. This might be sufficient for your purpose:
dat <- readLines(file("you
for the size of your data file, I think R can handle it. of course, it
also depends on your hardware. however, it might not be a good idea to
do heavy data manipulation work in R.
stata has very good routine for survey analysis. i am not sure if R is
as good as stata in terms of survey analysis.
steve bowdoin.edu> writes:
>
> I would like to get rid of the row numbers using xtable and latex.
>
> but I don't want the row numbers. Is it possible to get rid of them?
>
> Also, if x is a data frame, latex(x) contains the row numbers.
> Can I get rid of them here as well?
>
> Steve
Steve
Oh, Boy. This might result in a data dump since each of us has a personal
library. Here are the top dozen or so from mine:
1. Agresti, Alan, Categorical Data Analysis, 2nd ed., Wiley, 2002
2. Box, George E. P., William G. Hunter, and J. Stewart Hunter,
Statistics for Experimenters
On 01-Mar-07 bunny , lautloscrew.com wrote:
> Hi there,
>
> is there a possibility to count the number of appearances of an
> element in a vector ?
> i mean of any given element.. deliver all elements which are exactly
> xtimes in this vector ?
>
> thx in advance !!
If it is a specific eleme
I hope this question is sufficiently different from the other requests
for book recommendations that it's not repetitious. If not, I apologize
in advance.
I'm curious what standard reference books working statisticians, or
biostatisticians, have within easy reach of their desk. I'm a computer
syst
bunny , lautloscrew.com wrote:
> Hi there,
>
> is there a possibility to count the number of appearances of an
> element in a vector ?
> i mean of any given element.. deliver all elements which are exactly
> xtimes in this vector ?
>
> thx in advance !!
> myvec <- c(0, rep(1:6, each = 2), 7,
From: Chuck Cleland
>
> ahimsa campos-arceiz wrote:
> > Dear useRs,
> >
> > In a data.frame (df) I have several columns (x1, x2, x3xn)
> > containing data as a continuous numerical response:
> >
> > df
> > var x1x2 x3
> > 1143 147 137
> > 2 9393 117
> >
On Thu, 2007-03-01 at 15:20 +0100, bunny , lautloscrew.com wrote:
> Hi there,
>
> is there a possibility to count the number of appearances of an
> element in a vector ?
> i mean of any given element.. deliver all elements which are exactly
> xtimes in this vector ?
>
> thx in advance !!
Vec
I would like to get rid of the row numbers using xtable and latex.
The commands
d = cbind(1:10,rep(1:2,5))
ans = xtable(d)
latex(ans)
gives output containing
\begin{tabular}{rrr}
\hline
& 1 & 2 \\
\hline
1 & 1.00 & 1.00 \\
2 & 2.00 & 2.00 \\
3 & 3.00 & 1.00 \\
4 & 4.00 & 2
Take a look at "table".
Hope this helps,
Matt
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of bunny ,
lautloscrew.com
Sent: Thursday, March 01, 2007 9:20 AM
To: r-help@stat.math.ethz.ch
Subject: [R] count the # of appearances... [Broadcast]
Hi there,
is
hello,
I have plots http://uosis.mif.vu.lt/~roka5178/time series.JPG
I must to do almost the same. I do not have any data.
but i know that this is time series data with parabolic changeable mean.
I can not generate tie series data with parabolic mean.
how can i do this?
and have http://uosis.mif.v
Dear Bunny,
? table
might be what you wish.
Best wishes
Wolfgang Huber
EBI
bunny , lautloscrew.com wrote:
> Hi there,
>
> is there a possibility to count the number of appearances of an
> element in a vector ?
> i mean of any given element.. deliver all elements which are exactly
> x
?table
Best,
Dimitris
Dimitris Rizopoulos
Ph.D. Student
Biostatistical Centre
School of Public Health
Catholic University of Leuven
Address: Kapucijnenvoer 35, Leuven, Belgium
Tel: +32/(0)16/336899
Fax: +32/(0)16/337015
Web: http://med.kuleuven.be/biostat/
http://www.student.kuleuven.b
Hello Bunny/Lautioscrew,
sum(your vector == your chosen element) should do what you want...
HTH,
Ranjan
On Thu, 1 Mar 2007 15:20:19 +0100 "bunny , lautloscrew.com" <[EMAIL PROTECTED]>
wrote:
> Hi there,
>
> is there a possibility to count the number of appearances of an
> element in a vecto
On 01-Mar-07 R Heberto Ghezzo, Dr wrote:
>
> Hello,
> I do not know if I am completely out of it but . . .
> if x,y,z is a point in a sphere and [u,v,w]'A[u,v,w] = 1 is the
> equation of an ellipsoid and A = T'T (cholesky) then
> T.[x,y,z] should be a point in the ellipsoid ? isn't it?
> Heberto G
Dear all
I am fitting and analyzing linear mixed-effects models using the
R command 'lme'. The following is the results:
dental.fit <- lme(fixed = distance~age, random = ~age + cluster
= ~subject, data = dental)
> summary(dental.fit)
Variance/Covariance Components Estimates:
St
Hi there,
is there a possibility to count the number of appearances of an
element in a vector ?
i mean of any given element.. deliver all elements which are exactly
xtimes in this vector ?
thx in advance !!
__
R-help@stat.math.ethz.ch mailing list
R Heberto Ghezzo, Dr wrote:
>
> I do not know if I am completely out of it but . . .
> if x,y,z is a point in a sphere and [u,v,w]'A[u,v,w] = 1 is the
> equation of an ellipsoid and A = T'T (cholesky) then
> T.[x,y,z] should be a point in the ellipsoid ? isn't it?
>
Yes, it's a point _on_ the ell
On Thu, 2007-03-01 at 14:56 +0100, Marc A. Rohling wrote:
> Hello,
>
> I try to handle a simple bar-plot, but it turns out to be not as simple
> as I thought.
>
> 1) I have created a .dat-File, e.g. test.dat:
>
> DATA DATA-SEM
> 2.2 0.32
> 6.2 1.30
> 12.7 1.61
> 48.6 3.08
> 4.1 0.86
>
Hello,
I do not know if I am completely out of it but . . .
if x,y,z is a point in a sphere and [u,v,w]'A[u,v,w] = 1 is the equation of an
ellipsoid and A = T'T (cholesky) then
T.[x,y,z] should be a point in the ellipsoid ? isn't it?
Heberto Ghezzo
McGill University
Montreal Canada
_
Dear all,
i have a problem which deals with indexing again, i guess.
I am a beginner and though i read the S Poetry stuff it is not that
easy for me
to use the examples for my problems, maybe all this needs some
practice...
Here´s the question. As i want to perform a factor analysis i need
? par
it is the xpd you're looking for.
x <- seq(-pi, pi, len = 65)
par(xpd=TRUE)
plot(x, sin(x), type="l", col = 2,xpd=NA)
legend(x = 0, y = -1.5, "legend text", pch = 1, xjust
= 0.5)
--- Jenny Barnes <[EMAIL PROTECTED]> wrote:
> Hi folks,
>
> Do you mind if I ask a related question that I
On Thursday 01 March 2007 08:38, Dacha Atienza wrote:
> I want to do this model using mgcv
>
> caco1.gam<-gam(caco1$Pa~s(caco1$T10)+s(caco1$S10)+s(caco1$C10), data=caco1)
> But I want to indicate the degrees of freedom of each term, how I have to
> do it?
e.g.
s(caco1$T10,k=6,fx=TRUE)
will give yo
Hello,
I try to handle a simple bar-plot, but it turns out to be not as simple
as I thought.
1) I have created a .dat-File, e.g. test.dat:
DATADATA-SEM
2.2 0.32
6.2 1.30
12.71.61
48.63.08
4.1 0.86
4.5 0.32
1.5 1.13
1.2 1.08
The first row is the data represen
Don't know if this helps, but...
gam in package mgcv will let you set up smooths that interact with factors
using the `by' variable mechanism. See ?gam.models, particularly the last
example. Prediction is supported.
Simon
--
> Simon Wood, Mathematical Sciences, University of Bath, Bath, BA2 7
Dear colleagues,
I've just come across a problem with the following command which is a
part of the "metaOverview.R" code file provided as an monography-
accompanying file at
http://www.bioconductor.org/docs/mogr/metadata:
##
R> hasChr <- eapply(GOTERM, function(
Well, not extremely elegant, but should work:
1) open your file in some ascii text editor, delete the rubbish at the
beginning up to line Scan 1, and replace all spaces in names - e.g. make
a mass replace of 'Retention Time' by let say 'RetentionTime'.
2) Use read.table(), matrix() and data.fram
?readline
Note that this gives you a character vector, so you might want to use
as.numeric() after that.
Petr
d. sarthi maheshwari napsal(a):
> Hi,
>
> I am trying to read some value from keyboard in one of my R program. I want
> to do manipulation on data based on this value. But I am facing a
Dear R user,
I have a simple question for you: I have created a global function that
evoke other subsidiary functions. when I run the global function I want
to store the outcomes of the of each subsidiary function into a
variables. an examples
my.fun=function(vector, index){
a=fun.a(vector,
Hi,
I am trying to read some value from keyboard in one of my R program. I want
to do manipulation on data based on this value. But I am facing a problem
which is described below with the help of an example:
My code is similar to:
:
:
cat("\n","Enter value:: ","\n")
y<-scan(n=1)
a <- b*y
cat("\n
1 - 100 of 111 matches
Mail list logo