hits=-2.6 tests=BAYES_00
X-USF-Spam-Flag: NO
On Wed, 2008-02-06 at 22:57 -0800, Michael wrote:
> Hi all,
>
> I am playing around with the S-Plus ... and I am wondering how to
> connect it back and forth with Matlab?
>
> For example, how do I call S-Plus from within Matlab and how do I call
> Mat
Hi all,
Does anybody have the source code of stochastic volatility models in R
or Matlab, for example, the Bayesian based or the simulation based SV
estimations as described by Prof Eric Zivot in the following
discussion?
https://stat.ethz.ch/pipermail/r-sig-finance/2005q4/000501.html
--
Hi all,
I am playing around with the S-Plus ... and I am wondering how to
connect it back and forth with Matlab?
For example, how do I call S-Plus from within Matlab and how do I call
Matlab from within S-Plus?
How do I transfer data back and forth?
Thanks a lot
___
Hello,
I have a data frame with a factor column, which uniquely identifies
the observations in the data frame and it looks like this:
sample1_condition1_place1
sample2_condition1_place1
sample3_condition1_place1
.
.
.
sample3_condition3_place3
I want to turn it into three separate factor columns
What about doing more to promote the R Wiki? Someone complained
to me that the R Wiki is not very useful. However, I think with the
right encouragement, it could grow like the regular Wiki, which now has
over 9 million entries in over 250 languages.
I think the R Wiki is ideal fo
Open suggestion/question:
If you in each step of an K-step iteration load/allocate a large
object, each time of a different size, followed by smaller memory
allocations (due to your analysis), you might be better of if you
could do the iteration such that the largest object is in the first
iterati
On Feb 6, 2008 11:28 AM, Tony Plate <[EMAIL PROTECTED]> wrote:
> Bert Gunter wrote:
> > I strongly suggest you collaborate with a local statistician. I can think of
> > no circumstance where multiple regression on "hundreds of thousands of
> > variables" is anything more than a fancy random number
Dear Harold,
I had the same problem some times ago. I noticed that after I run a set
commands (cleaning all non-usefull variables) for 5 times, the system
broken-down. I solved it building several scritpsNN.R and call them in a .BAT
DOS file. It worked so fine, almost in my case, and the comput
What exactly are you intending the loop to do? Why do you have the
'as.matrix' in the middle of the loop? Where was 'y' defined? Does
this do what you want?
> outer(1:5, 1:10, "+")
[,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10]
[1,]23456789 1011
[
PLEASE do read the posting guide
http://www.R-project.org/posting-guide.html and provide commented,
minimal, self-contained, reproducible code.
>
It all depends on how you constructed the x-axis. You can try abline(v=2008.25)
On Feb 6, 2008 6:41 PM, Ana Quitério <[EMAIL PROTECTED]> wrote:
> De
Dear list,
I'm trying to make a loop of a (5x10) matrix and below are my codes. Could
anybody help me figure out why my loop is not working. Thanks in advance!!
m<-1:5
n<-1:10
for(i in 1:length(m))
{ for(j in 1:length(n))
{
y[i,j]=sum(i,j)
y<-as.matrix(y[i,j])
}
}
cheer
Dear all.
I want to add a vertical line in my time series series plot.
If i had an anual data I do:
plot(data,v="2008"), but if I want to add a vertical line in fisrt quarter
of 2008, how can i do?
Thanks a lot
Ana
--
__
R-help@r-project.org mailing l
Thank all of you for your helps. They are very helpful.
But I have a further question. Suppose I have the following mixed effect
model
thetaMixed <- function(tau, i)
{
w <- 1 / (s^2 + tau^2)
mu <- sum(theta * w) / sum(w)
b <- s[i]^2 / (s[i]^2 + tau^2)
theta[i]*(1-b) + mu*b)
}
A more primitive method is about 5 times faster than Gabor's.
L <- list(
a = c("1", "2", "3"),
b = c("1"),
d = c("2", "4")
)
system.time(
for (i in 1:100)
{t1 <- unlist(L)
names(t1) <- rep(names(L), lapply(L, length))
tapply(names(t1), t1, c)
}
)
system.time(
This isn't a single command but its pretty short:
unstack(stack(L)[2:1])
On Feb 6, 2008 5:51 PM, hadley wickham <[EMAIL PROTECTED]> wrote:
> Is there a built in function to invert a list? i.e. to go from
>
> list(
> a = c("1", "2", "3"),
> b = c("1"),
> d = c("2", "4")
> )
>
> to
>
> list(
>
Is there a built in function to invert a list? i.e. to go from
list(
a = c("1", "2", "3"),
b = c("1"),
d = c("2", "4")
)
to
list(
"1" = c("a", "b"),
"2" = c("a", "d"),
"3" = "a",
"4" = "2"
)
Hadley
--
http://had.co.nz/
__
R-help@r-projec
>Tim Hesterberg wrote:
>>I'll raise a related issue - sampling with unequal probabilities,
>>without replacement. R does the wrong thing, in my opinion:
>>...
>Peter Dalgaard wrote:
>But is that the right thing? ...
(See bottom for more of the previous messages.)
First, consider the common case,
Oh yes, that's what I was looking for. I can then do a
Class1<-rbind(Informative.data, Uninformative.data)
Thanks ../Murli
-Original Message-
From: Henrique Dallazuanna [mailto:[EMAIL PROTECTED]
Sent: Wednesday, February 06, 2008 3:06 PM
To: Nair, Murlidharan T
Cc: [EMAIL PROTECTED]
Subj
Dear R users,
Does anyone know of a way to obtain approximate 95% confidence intervals
for predicted values for factor levels of fixed effects from lme? Our
goal is to use these intervals to interpret patterns across our
predicted values for certain factor levels.
Our mixed model has the foll
Perhaps:
for(i in 1:100)assign(sprintf("Informative.data.class%s", i),
rnorm(100, 0.25,1))
for(i in 1:100)assign(sprintf("Uninformative.data.class%s", i), rnorm(900))
Or working with a list:
Informative <- replicate(100, rnorm(100, 0.25,1))
Uninformative <- replicate(100, rnorm(900))
On 06/02/2
Dear Clara,
The warning is nothing to worry about. If start values aren't specified,
box.cox.powers() uses optimize() to find start values for each
transformation parameter (in this case, there's only one) prior to using
optim() to maximize the (possibly multivariate) likelihood. In the course of
On 2008-02-06, Birgit Lemcke <[EMAIL PROTECTED]> wrote:
>
> I am using R 2.6.1 on a PowerBook G4.
> I would like to perform a discriminant function analysis. I found lda
> in MASS but as far as I understood, is it only working with
> explanatory variables of the class factor.
I think you are
I am trying to generate artificial data for feature selection. Basically
trying to generate a total of 1000 features with 100 that are informative and
rest are uninformative.
Informative.data.class1<-rnorm(100,0.25,1)
Uninformative.data.class1<-rnorm(900,0,1)
Informative.data.class2<-rnorm(100,
Tim Hesterberg wrote:
>> values <- sapply(1:1000, function(i) sample(1:3, size=2, prob = c(.5, .25,
>> .25)))
>> table(values)
>>
> values
> 1 2 3
> 834 574 592
>
> The selection probabilities are not proportional to the specified
> probabilities.
>
> In contrast, in S-PLUS:
>
>>
On Wednesday 06 February 2008 (15:14:28), [EMAIL PROTECTED] wrote:
> Hi, i'm having trouble with my x and y axis. The commands i'm using are
>below. The problem is that the y axis starts at coordinate 0,1 and the x
>axis starts at coordinate 0,0. As far as I know the y axis can't start
how can i compare a cox model and a rsf?
someone can tell me how to do that in R?
Thanks
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting
This should do it for you:
x <- read.table(textConnection("V1V2 V3
1 chr1 11255 55
2 chr1 11320 29
3 chr1 11400 45
4 chr2 21680 35
5 chr2 21750 84
6 chr2 21820 29
7 chr2 31890 46
8 chr3 32100 29
9 chr3 52380 29
10 chr3 66450 46" ), header=TRUE)
outfile <- '/tempxx.txt'
cat("browser positio
On Feb 6, 2008 8:08 AM, Waterman, DG (David)
<[EMAIL PROTECTED]> wrote:
> Hi,
> I have a data frame consisting of coordinates on a 10*10 grid, i.e.
> > example
> x y
> 1 4 5
> 2 6 7
> 3 6 6
> 4 7 5
> 5 5 7
> 6 6 7
> 7 4 5
> 8 6 7
> 9 7 6
> 10 5 6
> What I would
> I want to generate different samples using the
>followindg code:
>
>g<-sample(LETTERS[1:2], 24, replace=T)
>
> How can I specify that I need 12 "A"s and 12 "B"s?
I introduced the concept of "sampling with minimal replacement" into the
S-PLUS version of sample to handle things like this:
zoo has a NAMESPACE which restricts access to the
internals of the package except for explicitly exported items.
zoo exports plot.zoo as an S3 method so that the plot generic
can find it. It is also possible to simply export a function
not specifically as an S3 method in which case it will be
visi
Gabor, maybe I am not understanding it right. I was thinking for example
something like how a call to plot actually calls plot.zoo when zoo is
loaded. And a specific call to plot.zoo would not work, etc. One would think
that plot and plot.zoo are separate and that one could call the parent as
well
Dear Rlist,
Using an example in box.cox.powers() help, I have the following warning message.
example:
library(car)
>attach(Prestige)
> box.cox.powers(income)
Box-Cox Transformation to Normality
Est.Power Std.Err. Wald(Power=0) Wald(Power=1)
0.1793 0.11081.6179 -7.4062
L.R
thanks Richie. It woks perfectly.
- Original Message
From: "[EMAIL PROTECTED]" <[EMAIL PROTECTED]>
To: joseph <[EMAIL PROTECTED]>
Cc: r-help@r-project.org; [EMAIL PROTECTED]
Sent: Wednesday, February 6, 2008 2:22:16 AM
Subject: Re: [R] inserting text lines in a dat frame
>
I
am
try
If you believe that certain changes intended only to affect the local
environment nevertheless affect the global environment please give a
code example.
On Feb 6, 2008 12:11 PM, tom soyer <[EMAIL PROTECTED]> wrote:
> Thanks Hardley. I see what you mean. You are right, I am not an expert in
> oop A
I am trying to obtain the Kruskal (1964) secondary least-squares monotonic
transformation of a rank variable given 4 categorical variables in order to
obtain optimal transformation for regression. The academic problem assigned
is to compare R, SPSS (Conjoint Analysis), and SAS' proc transreg in spe
Hi Jim
yes, this exactly want I want. However, I need one more step to get
rid of the first column so that the final file looks like this:
browser
position
chr1:1-1
browser
hide
all
track
type=wiggle_0
name=sample
description=chr1_sample
visibility=full
variableStep
chrom=chr1
span=
On 06-Feb-08 17:09:59, Peter Dalgaard wrote:
> (Ted Harding) wrote:
>> Greetings!
>>
>> I suspect that there is an error in the code for the
>> function ci.pd() in the Epi package.
>>
> Any particular reason not to include the maintainer among the
> recipients??
Oversight, and haste! Thank you
Bert Gunter wrote:
> I strongly suggest you collaborate with a local statistician. I can think of
> no circumstance where multiple regression on "hundreds of thousands of
> variables" is anything more than a fancy random number generator.
That sounds like a challenge! What is the largest regress
Thanks Hardley. I see what you mean. You are right, I am not an expert in
oop AND I don't really know how R oo works, so certainly I shouldn't be
making any sweeping statement. I was just thinking about the issue of local
vs. global, i.e. changes intended for the local environment shouldn't affect
(Ted Harding) wrote:
> Greetings!
>
> I suspect that there is an error in the code for the
> function ci.pd() in the Epi package.
>
Any particular reason not to include the maintainer among the
recipients?? Not sure whether Bendix is subscribed.
-p
> This function is for computing confide
Hi.
I have a AFTD plot. I want to add a original axis for x and y.
I did "help(plot)". But I did not see an option for that.
I need help.
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/
Greetings!
I suspect that there is an error in the code for the
function ci.pd() in the Epi package.
This function is for computing confidence intervals
for a difference of proportions between two independent
groups of 0/1 responses, and implements the Newcombe
("Nc") method and the Agrasti-Caffo
On Feb 6, 2008 10:13 AM, tom soyer <[EMAIL PROTECTED]> wrote:
> Thanks Gabor. I guess true oo encapsulation is not possible in R.
Before making such a claim, I would encourage you to actually learn
what oo means. A couple of good readings are:
Structure and Interpretation of Computer Programs.
h
Thanks,
I managed to run gls, but my problem isn't solved :)
1. I do not know the autocorrelation between errors and its estimates
change for different methods. Therefore I wasn't including it in the
corAR1 function.
2. GLS yields counteraintuitive results with one of the variables (its
sign) an
> Yes, type
>
> vignette("s3d", package="scatterplot3d")
>
> and see page 22.
which of course won't work unless you have scatterplot3d installed,
but I wonder if vignette() could be modified to try to find the
vignette on the web if the package isn't installed. It would make it
a little easy to r
> Is there anyway to produce a 3D scatterplot with error bars in the x,y,z
> directions? I have searched around and know of scatterplot3d but did
> not see any way to put error bars on the points. Any ideas?
It may be possible, but do you think that's the best solution for your
problem? Given t
On 2/6/08, Alex Brown <[EMAIL PROTECTED]> wrote:
> Here's a related problem that works for numeric but not for POSIXct
>
> I am seeing it where a panel has no at labels, but others do.
>
> This simple example only has one panel with no at labels.
>
> > baseval = 0;
> > xyplot(1:10 ~ (baseval + c(
Christopher Chizinski wrote:
> Is there anyway to produce a 3D scatterplot with error bars in the x,y,z
> directions? I have searched around and know of scatterplot3d but did
> not see any way to put error bars on the points. Any ideas?
Yes, type
vignette("s3d", package="scatterplot3d")
a
Shubha,
The calculation is:
wilderSum[1] <- x[1]
for(i in 2:NROW(x)) {
wilderSum[i] <- x[i] + wilderSum[i-1] * (n-1)/n
}
Where 'x' is the price series and 'n' is the number of periods.
You can see the calculations for all functions in TTR's source code,
which is on CRAN and r-forge
(http://r-
i returned to this problem and found that the solution was a disarmingly
straight-forward oversight on my part: 'detach' requires the version
number of the library to be specified!
thus, if you install packages '--with-package-versions' switched on,
and if you load the version of choice in an R
It seems like you didn't look at the examples in the helpfiles. See ?gls
and ?corAR1.
ir. Thierry Onkelinx
Instituut voor natuur- en bosonderzoek / Research Institute for Nature
and Forest
Cel biometrie, methodologie en
On 2/5/08, Kelvin <[EMAIL PROTECTED]> wrote:
> I am trying to create levelplot's of cpu usage for systems.
> print(levelplot(util.mean ~ x.hour * x.day, colorkey=T, cut=20,
> scales=list(x=list(at=seq(0,96,length=25),
> labels=ifelse(seq(0,24) %% 4 == 0, seq(0,24), ''))), # add
> tick
Hello R-Cracks,
I am using R 2.6.1 on a PowerBook G4.
I would like to perform a discriminant function analysis. I found lda
in MASS but as far as I understood, is it only working with
explanatory variables of the class factor. My dataset contains
variables of the classes factor and numeric.
Thanks Gabor. I guess true oo encapsulation is not possible in R.
It seems that there is an IDE for S+ in Eclipse...
On 2/6/08, Gabor Grothendieck <[EMAIL PROTECTED]> wrote:
>
> On Feb 6, 2008 9:45 AM, tom soyer <[EMAIL PROTECTED]> wrote:
> > Thanks Gabor for illustrating the basics oop in R us
Richard Pearson wrote on 02/06/2008 06:25 AM:
> Hi Thomas
[...]
> With databases, one issue that might be relevant is whether you want to
> store data in tables (e.g. one table to store one data.frame) that can
> subsequently be manipulated in the DB, or to store R objects as R
> objects (e.g. a
"tom soyer" <[EMAIL PROTECTED]> writes:
> (1) how do I encapsulate the generics? i.e., if a class has 100 methods,
> then does it mean 100 generics would be dumped in the global environment?
> Or, is it possible to define a local environment and restrict the generics
> from one class to a particula
Normally I can run an R script in batch mode with a command like this
R CMD BATCH MyScript.R MyOutput &
However I prefer to write another script containing something like
R --no-restore --save --no-readline < $1 >$2
so that I could run the original script simply on the prompt as
MyScript.R MyO
Hi,
I was wondering if someone might be willing to indulge a question about
R and the estimation of a linear regression with time-varying coefficients.
The model I am trying to estimate is of the form:
y(t) = beta(t) * x(t) + v(t)
beta(t) = gamma * beta(t-1) + w(t)
where gamma is a constant, v(
On Wed, Feb 06, 2008 at 02:34:52PM +0200, Rainer M Krug wrote:
> R objects in blobs - I never thought about that. Could you elaborate on how
> to do something like that (I am using RMySQL)?
Look at help(serialize) -- any R object can be turned into a suitable
representation, either binary (more ef
OK I got gls running but this is the error message that I got:
> gls_results <- gls(LnWRPK ~ LnWGDP+LnWYIELD,correlation = corAR1)
Error in switch(mode(x), "NULL" = structure(NULL, class = "formula"), :
invalid formula
Google didn't give any results unfortunately.
Help! :)
Thanks,
KB
> Ha
Hi
Im using the kinship package to draw a pedigree. On my data set this works fine
but when i add indivudals to the pedigree i keep getting an error i hope
someone can help me!
This is the code im using:
Data<-read.table("Tree.txt", header=T, sep=",")
attach(Data)
ped<-pedigree(id, dadid, momi
Here's a related problem that works for numeric but not for POSIXct
I am seeing it where a panel has no at labels, but others do.
This simple example only has one panel with no at labels.
> baseval = 0;
> xyplot(1:10 ~ (baseval + c(1:10)), scales=list(x=list
+ (at=list(c()), labels=list(c()),
On Feb 6, 2008 9:45 AM, tom soyer <[EMAIL PROTECTED]> wrote:
> Thanks Gabor for illustrating the basics oop in R using S3. Maybe I didn't
> have the right documents, but you example taught me more about oop in R than
> everything else I read combined! Thanks for the tip on R.oo, I plan to check
> i
Is there anyway to produce a 3D scatterplot with error bars in the x,y,z
directions? I have searched around and know of scatterplot3d but did
not see any way to put error bars on the points. Any ideas?
__
R-help@r-project.org mailing list
https://sta
See this example:
coef(summary(glm(Postwt ~ Prewt + Treat +
offset(Prewt),family=gaussian, data=anorexia)))
Is this you want?
On 06/02/2008, Munyandorero, Joseph <[EMAIL PROTECTED]> wrote:
> Dear all,
>
> After running a glm, I use the summary ( ) function to extract its
> coefficients and rela
If I understand correctly.. using the example in ?glm
> utils::data(anorexia, package="MASS")
> anorex.1 <- glm(Postwt ~ Prewt + Treat + offset(Prewt),
family = gaussian, data = anorexia)
> summ <- summary(anorex.1)
Then check
> names(summ)
for the available objects.
F
On 2/6/2008 10:07 AM, Pedro Mardones wrote:
> Dear all;
> I've been trying to change the type of line used to draw the box
> around the 3d scatterplot (package scatterplot3d) from lty=1 to lty=2
> without sucess. I would appreciate suggestions of how to do it.
Use lty.axis=2. (This is mentioned o
Dear all;
I've been trying to change the type of line used to draw the box
around the 3d scatterplot (package scatterplot3d) from lty=1 to lty=2
without sucess. I would appreciate suggestions of how to do it.
Thanks
PM
__
R-help@r-project.org mailing lis
Dear all,
After running a glm, I use the summary ( ) function to extract its
coefficients and related statistics for further use. Unfortunately, the
screen only displays a small (last) part of the results. I tried to
overcome the problem by creating/saving an object "coef" for
coefficients of the
Sorry, word wrap made that incomprehensible, I think
x y
4 5
6 7
6 6
7 5
5 7
6 7
4 5
6 7
7 6
5 6
dat <- read.table('clipboard', header=TRUE)
dat$patt <- paste(dat$x,dat$y, sep='')
mm <- as.data.frame(with(dat, table(patt)))
dat <- merge(dat, mm, by='patt')
mat <- matrix(0, ncol=10,
I think this does what you want, but there may be a more efficient way
x y
4 5
6 7
6 6
7 5
5 7
6 7
4 5
6 7
7 6
5 6
dat <- read.table('clipboard', header=TRUE) # copy sample data above
dat$patt <- paste(dat$x,dat$y, sep='')
mm <- as.data.frame(with(dat, table(patt)))
dat <- merge(dat, mm
On Wednesday 06 February 2008 14:08, Waterman, DG (David) wrote:
> Hi,
>
> I have a data frame consisting of coordinates on a 10*10 grid, i.e.
>
> > example
>
> x y
> 1 4 5
> 2 6 7
> 3 6 6
> 4 7 5
> 5 5 7
> 6 6 7
> 7 4 5
> 8 6 7
> 9 7 6
> 10 5 6
>
> What I would li
On 2008-February-06 , at 14:45 , Duncan Murdoch wrote:
> On 2/5/2008 11:12 AM, jiho wrote:
>> Dear List,
>> I am trying to get R's terminal output to a file and to the
>> terminal at the same time, so that I can walk through some tests
>> and keep a log concurrently. The function 'sink' with
Thanks Gabor for illustrating the basics oop in R using S3. Maybe I didn't
have the right documents, but you example taught me more about oop in R than
everything else I read combined! Thanks for the tip on R.oo, I plan to check
it out later.
I have a few followup questions...
(1) how do I encaps
John Fox wrote:
> Dear Peter and Iksmax,
>
> To elaborate slightly, the Rcmdr tries to figure out which menu items are
> appropriate in a given context, and as Peter says, requires that you have at
> least one factor in the active dataset before activating the pie chart menu
> item; only factors wi
Hi,
Can somebody tell me the formula for "?wilderSum" in TTR package? I mean how
are these calculated?
BR, Shubha
Shubha Karanth | Amba Research
Ph +91 80 3980 8031 | Mob +91 94 4886 4510
Bangalore * Colombo * London * New York * San José * Singapore *
www.ambaresearch.com
This e-
There may be an easier way but you can extract the
desired values from the list values in t
str(t) to see the elements in t.
test <- matrix(c(1, 1,2,2), 2,2)
tt <- apply(test, 1, t.test)
ttable <- function(tlist) {
tframe <- data.frame(NULL)
for(i in 1:length(tlist)){
t.value <- tlist[[i]]
Hi, i'm having trouble with my x and y axis. The commands i'm using are
below. The problem is that the y axis starts at coordinate 0,1 and the x
axis starts at coordinate 0,0. As far as I know the y axis can't start at 0
(because it's log scaled) ,so I would like to position the x a
Hi,
I have a data frame consisting of coordinates on a 10*10 grid, i.e.
> example
x y
1 4 5
2 6 7
3 6 6
4 7 5
5 5 7
6 6 7
7 4 5
8 6 7
9 7 6
10 5 6
What I would like to do is return an 10*10 matrix consisting of counts
at each position, so in the above example
On 2/5/2008 11:12 AM, jiho wrote:
> Dear List,
>
> I am trying to get R's terminal output to a file and to the terminal
> at the same time, so that I can walk through some tests and keep a log
> concurrently. The function 'sink' with the option split=TRUE seems to
> do just that. It works fi
There is no generally agreed upon notion of random effects for quantile
regression applications. Insofar as one is willing to accept the
idea that
random effects are just "shrunken fixed effects" one can consider
similar
schemes in the QR context; one such is described in
“Quantile Regressio
> A new user will come to the R homepage, go to CRAN via the link under
> download and from there see Packages and then be swamped by the huge
> number available. Having Task Views as a link on the R homepage would
> make these more visible.
I would think that a new user would see the download hea
Well, not exactly. package.skeleton() is very useful as a first step,
but it does *not* create a package entirely. It turns out that in
Windows, creating a package is very simple once you have downloaded all
programs needed (e.g., perl) and you have your path configured exaclty,
(exactly, exactly)
I'm not sure that this would make any difference to someone considering using R.
Would they know what CRAN stands for? Probably not unless they've used
CPAN or equivalent in the past.
Would they know what a 'Task View' is? Again probably not as its not
patently obvious what it is, it doesn't "cli
With databases, one issue that might be relevant is whether you want to
> store data in tables (e.g. one table to store one data.frame) that can
> subsequently be manipulated in the DB, or to store R objects as R
> objects (e.g. as BLOBs). My situation is likely to be the later case,
> and one of m
Hi Thomas
I'm certainly no expert but thought I'd reply as I'm likely to be in a
similar position soon.
With regards versions of R I think you should always have the latest
release version. This will mean upgrading at least every 6 months, but
this shouldn't be too much of a problem.
With OSs
hits=-2.6 tests=BAYES_00
X-USF-Spam-Flag: NO
On Wed, 2008-02-06 at 12:45 +0100, Konrad BLOCHER wrote:
> I get this message:
>
> Error: could not find function "gls" (and also)
> Error: could not find function "lm.gls"
>
> Which package is that in?
RSiteSearch("gls", restrict = "functions")
Tel
On Feb 4, 08:49 AM, Chuck Cleland wrote:
> On 2/3/2008 10:09 AM, Christoph Mathys wrote:
>> Dear R users,
>> I have a linear model of the kind
>> outcome ~ treatment + covariate
>> where 'treatment' is a factor with three levels ("0", "1", and "2"),
>> and the covariate is continuous. Treatments "1
On Feb 6, 2008, at 6:23 AM, Neil Shephard wrote:
> Charilaos Skiadas-3 wrote:
>>
>> On Feb 5, 2008, at 10:37 AM, Monica Pisica wrote:
>>
>> But perhaps I am missing something very obvious?
>
> I thought the task views were located where they are (linked from
> the page
> that lists packages) as
I get this message:
Error: could not find function "gls" (and also)
Error: could not find function "lm.gls"
Which package is that in?
Thanks,
KB
> Have you tried gls()?
>
>
>
>
> ir. Thierry Onkelinx
> Instituut voor
hits=-2.6 tests=BAYES_00
X-USF-Spam-Flag: NO
On Wed, 2008-02-06 at 03:23 -0800, Neil Shephard wrote:
>
>
> Charilaos Skiadas-3 wrote:
> >
> > On Feb 5, 2008, at 10:37 AM, Monica Pisica wrote:
> >
> > But perhaps I am missing something very obvious?
> >
> >
>
> I thought the task views were
On 2008-February-06 , at 11:25 , Alex Brown wrote:
> you could use the unix function 'script' before invoking the R
> interpreter.
Thanks for the suggestion. It would work for some cases and I did not
know about this utility. But most of the time I just put a list of
commands in a .R file
Charilaos Skiadas-3 wrote:
>
> On Feb 5, 2008, at 10:37 AM, Monica Pisica wrote:
>
> But perhaps I am missing something very obvious?
>
>
I thought the task views were located where they are (linked from the page
that lists packages) as they summarise the available packages for the given
to
Hi,
I am trying to perform Maximum Likelihood estimation of a Multivariate
model (2 independent variables + intercept) with autocorrelated errors of
1st order (ar(1)).
Does R have a function for that? I could only find an univariate option
(ar.mle function) and when writing my own I find that it
On Feb 6, 2008 2:35 AM, ONKELINX, Thierry <[EMAIL PROTECTED]> wrote:
> You'll need to transform your dataset in a long format first.
>
> library(ggplot2)
> n <- 5
> MyValues <- data.frame(Gene = factor(LETTERS[seq_len(n)]), ES =
> rnorm(n), MEF = rnorm(n), Embrio = rnorm(n), EShyp = rnorm(n))
> MyV
Try this and see if it is what you want:
x <- read.table(textConnection(" V1V2 V3
1 chr1 11255 55
2 chr1 11320 29
3 chr1 11400 45
4 chr2 21680 35
5 chr2 21750 84
6 chr2 21820 29
7 chr2 31890 46
8 chr3 32100 29
9 chr3 52380 29
10 chr3 66450 46" ), header=TRUE)
cat("browser position chr1:1-1
Is 'thr' supposed to be the mean and sd of all the values in data2_1?
If so, then
thr <- mean(data2_1, na.rm=TRUE) + sd(data2_1,na.rm=TRUE)
I am not exactly sure of "what is the problem that you are trying to
solve". You just have to make sure that the object you are creating
by precomputing has
You'll need to transform your dataset in a long format first.
library(ggplot2)
n <- 5
MyValues <- data.frame(Gene = factor(LETTERS[seq_len(n)]), ES =
rnorm(n), MEF = rnorm(n), Embrio = rnorm(n), EShyp = rnorm(n))
MyValuesMelt <- melt(MyValues, id.var = "Gene")
ggplot(MyValuesMelt, aes(x = Gene, y
you could use the unix function 'script' before invoking the R
interpreter.
example session:
$ script
Script started, file is typescript
[x86_64|[EMAIL PROTECTED]:~]
$ R --quiet --vanilla
> 1:10
[1] 1 2 3 4 5 6 7 8 9 10
> q()
[x86_64|[EMAIL PROTECTED]:~]
$ exit
exit
Script done, file
> I am trying to prepare a bed file to load as accustom track on the
> UCSC genome browser.
> I have a data frame that looks like the one below.
> > x
> V1V2 V3
> 1 chr1 11255 55
> 2 chr1 11320 29
> 3 chr1 11400 45
> 4 chr2 21680 35
> 5 chr2 21750 84
> 6 chr2 21820 29
> 7 chr2 31890 46
>
1 - 100 of 107 matches
Mail list logo
