Hello
Well I am relatively new so some of these issues may not fall under the subject
that I have used.
1. How do I do a Pareto. Following is the approach I took.
My data looks like this
df2_9
Reaason.for.failure Frequency
1 Phy Conn1
2Power failure3
On Sat, 2007-12-08 at 22:51 +0100, Irene Mantzouni wrote:
> Hi all!
>
> I am fitting a (mixed) model with a factor (F) and continuous response and
> predictor:
> y~F+F:x
>
> (How) can I check the significance of the model at each factor level (i.e.
> the model could be significant only at on
Prof Brian Ripley wrote:
> Hmm, S'S is numerically singular. crossprod(S) would be a better way to
> compute it than crossprod(S,S) (it does use a different algorithm), but
> look at the singular values of S, which I suspect will show that S is
> numerically singular.
>
> Looks like the answer
Hi all,
Since I was having several problems trying to work with a large matrix I
started to use the package R.huge but I'm having the following problem
> x<-FileByteMatrix("covtype.data",nrow=581012,ncol=55)
Error: cannot allocate vector of size 770.8 Mb
In addition: Warning messages:
1: Reach
What you say about mixture models is true in general, however this fit was
the best of 100 random EM starts. Unbounded likelihoods I believe are only
a problem for continuous data mixture models and mine was discrete. Anyway
it's nearly midnight now here so I'd better sleep. Before I go, here are
t
Hi,
I'm sorry, but I will most likely withdraw R.huge from CRAN anytime
soon. The File***Matrix classes had some problems, which when solved
made it too slow. See BufferedMatrix package in Bioconductor instead.
/Henrik
On 09/12/2007, [EMAIL PROTECTED]
<[EMAIL PROTECTED]> wrote:
>
>
> Hi all,
>
Dear all,
I'm having difficulty getting access to data generated by survfit and
print.survfit when they are using with a Cox model (survfit.coxph).
I would like to programmatically access the median survival time for
each strata together with the 95% confidence interval. I can get it on
screen, b
Mark Wardle wardle.org> writes:
> I'm having difficulty getting access to data generated by survfit and
> print.survfit when they are using with a Cox model (survfit.coxph).
>
> I would like to programmatically access the median survival time for
> each strata together with the 95% confidence in
So S is numerically singular: the last 10 or 11 values are effectively 0.
As Peter Dalgaard said, scaling may help but only if the parameters are on
drastically different scales.
On Sun, 9 Dec 2007, [EMAIL PROTECTED] wrote:
> What you say about mixture models is true in general, however this fit
Marc Schwartz wrote:
> On Sun, 2007-12-09 at 12:26 +1100, Tim Churches wrote:
>> Marc Schwartz wrote:
>>> Perhaps the most notable format that is lacking is the SAS proprietary
>>> format (not the Transport format), which is not openly published and to
>>> my knowledge, has not been independently r
Marc Schwartz wrote:
> On Sun, 2007-12-09 at 01:04 +, Jonas Stein wrote:
>> Hi,
>>
>> what is actually the best method to include R-plots into LaTeX
>> documents?
>> At the moment i use
>>
>> postscript("myplot.eps", width = 12.0, height = 9.0, horizontal =
>> FALSE,
>> onefile = TRUE, paper
Dieter Menne wrote:
> Mark Wardle wardle.org> writes:
>
>> I'm having difficulty getting access to data generated by survfit and
>> print.survfit when they are using with a Cox model (survfit.coxph).
>>
>> I would like to programmatically access the median survival time for
>> each strata togethe
Hi,
I would like to obtain the names of all objects that are provided as further
arguments ("...") in a function call. I have created a minimal example that
illustrates my wish (but isn't useful otherwise):
f1 <- function( ... ) return( deparse( substitute( ... ) ) )
x1 <- 1
x2 <- 2
x3 <- 3
f
On 09/12/2007 8:28 AM, Arne Henningsen wrote:
> Hi,
>
> I would like to obtain the names of all objects that are provided as further
> arguments ("...") in a function call. I have created a minimal example that
> illustrates my wish (but isn't useful otherwise):
>
> f1 <- function( ... ) return
I haven't read about R's math-notation facilities in this thread, hence
my question for Jonas Stein: Have you already looked into ?plotmath ?
Uwe Ligges
Frank E Harrell Jr wrote:
> Marc Schwartz wrote:
>> On Sun, 2007-12-09 at 01:04 +, Jonas Stein wrote:
>>> Hi,
>>>
>>> what is actually the
Dear Everybody,
Enclosed are the tiny beginning of a program and its output.
As a pity, if I load the image the elements of
are non numerical; returns .
How can I keep the numbers numerical?
Thank you in advance.
hopefully,
Mag. Ferri Leberl
#Programm zur Erstellung eines Wasserwegenetzplanes
#
You program contains:
de.dd<-c("Dresden",51,2,13,44)
and since you are mixing characters with numerics and putting them
into a single vector, the data is coerced to a common data type, which
in this case is numeric. What are you trying to do? Do you want to
use a list?
> de.dd
[1] "Dresden" "5
Meant to say coerced to "character".
See "?c"
On Dec 9, 2007 6:14 AM, Mag. Ferri Leberl <[EMAIL PROTECTED]> wrote:
> Dear Everybody,
> Enclosed are the tiny beginning of a program and its output.
> As a pity, if I load the image the elements of
> are non numerical; returns de.dd[2] + de.dd[3]
Kapoor, Bharat wrote:
> Hello
>
> Well I am relatively new so some of these issues may not fall under the
> subject that I have used.
>
> 1. How do I do a Pareto. Following is the approach I took.
>
> My data looks like this
> df2_9
>Reaason.for.failure Frequency
>
> 1 Ph
Hi all
Two questions:
- I would like to learn more on oriented object programming with R. Is
there any tutorial for that?
- Without waking up a troll, I am not very familiar with diffusion
list, I am more use to forum. On a diffusion list, how can I check if
someone already asks a question? I
I have the following code
library(zoo)
miedate <- yearmon((2006)+seq(0,23)/12)
tab <- zoo(cbind(A = c(79.47, 89.13, 84.86, 75.68, 72.82, 78.87, 93.46,
78.18, 82.46, 77.25, 80.95, 84.39, 81.7, 74.76, 65.29, 60.3,
66.59, 73.19, 9
On 09/12/2007 12:21 PM, Christophe Genolini wrote:
> Hi all
>
> Two questions:
> - I would like to learn more on oriented object programming with R. Is
> there any tutorial for that?
The usual reference is the book
_Programming with Data_ (John M. Chambers, Springer, 1998)
I don't know of a
Hi all,
Is there a mailing list for spatstat R package?
Another question: Can I find some ebooks on Spatial Pattern Analysis to be
downloaded?
Finally, how can I use a shapefile as "mask" on spatstat? I read the polygon
using maptools::read.shape function and now I would like to use this shape
On Dec 9, 2007 12:21 PM, Christophe Genolini <[EMAIL PROTECTED]> wrote:
> Hi all
>
> Two questions:
> - I would like to learn more on oriented object programming with R. Is
> there any tutorial for that?
> - Without waking up a troll, I am not very familiar with diffusion
> list, I am more use to
On Sun, 9 Dec 2007, Milton Cezar Ribeiro wrote:
> Hi all,
>
> Is there a mailing list for spatstat R package?
The website is http://www.spatstat.org (seems to be down for me, but
that's probably me); the maintainers also read this list.
> Another question: Can I find some ebooks on Spatial Patt
>From: Duncan Murdoch <[EMAIL PROTECTED]>
>Date: 2007/12/09 Sun AM 11:40:45 CST
>To: Christophe Genolini <[EMAIL PROTECTED]>
>Cc: r-help@r-project.org
>Subject: Re: [R] Oriented object programming
Roger Peng has some objected oriented documentation
at his site as part of a powerpoint presentation.
- You don't need a panel function since there is only one.
- There should not be more than 2 components on lty and lwd
since tab has two columns and that is the cause of the
warning message you are seeing.
- with the recent version of zoo plot.type uses screen so
you can write the shorter screen =
Charilaos Skiadas wrote:
>
> On Dec 8, 2007, at 1:02 AM, Joe W. Byers wrote:
>
>> In a post on R-devel, Prof Ripley add the following comment
>> | > BTW, 1:dim(names)[1] is dangerous: it could be 1:0. That was the
>> | > motivation for seq_len.
>>
>> I use the dim(names)[1] and dim(x)[2] along wit
Thank you. I have downloaded the compressed package "growth-1.tgz". I have
extracted the archive. But all my attempts to install such a package so that
I can call its functions from R command line ... have failed.
It is not clear to me how to operate the "install.packages" command where
the package
Hi Folks,
Say I have 2 continuous variables X,Y.
I can of course plot (X,Y) with
plot(X,Y,pch="+",col="blue")
say, and add the regression line from lm(Y~X)
by extracting the coefficients 'a' of Intercept
and 'b' of X from Y.lm <- lm(Y~X).
Now, however, I want to have not only a general
expla
On Sun, 2007-12-09 at 13:08 -0600, Maura E Monville wrote:
> Thank you. I have downloaded the compressed package "growth-1.tgz". I
> have extracted the archive. But all my attempts to install such a
> package so that I can call its functions from R command line ... have
> failed.
> It is not clear
If x <- round(...whatever...) try:
title(capture.output(print(x)))
On Dec 9, 2007 2:15 PM, Ted Harding <[EMAIL PROTECTED]> wrote:
> Hi Folks,
>
> Say I have 2 continuous variables X,Y.
>
> I can of course plot (X,Y) with
>
> plot(X,Y,pch="+",col="blue")
>
> say, and add the regression line from
On Sun, 2007-12-09 at 19:15 +, [EMAIL PROTECTED] wrote:
> Hi Folks,
>
> Say I have 2 continuous variables X,Y.
>
> I can of course plot (X,Y) with
>
> plot(X,Y,pch="+",col="blue")
>
> say, and add the regression line from lm(Y~X)
> by extracting the coefficients 'a' of Intercept
> and 'b
Hey
do anyone know why this error occurs??
> for(i in 1:n_trait){
+for( j in 1:n_trait){
+rG[i,j] <- (G_o[i,j]/(sqrt(G_o[i,i]%*%G_o[j,j]))
+rP[i,j] <- (P_o[i,j]/(sqrt(P_o[i,i]%*%P_o[j,j]))
Error: unexpected symbol in:
" rG[i,j] <- (G_o[i,j]/(sqrt(G_o[i,i]%*%G_o[j,j]))
rP"
>rE
I would like to know if is it possible implemet a partitioning tree
using a function like rpart, or mvpart, and with formula a glm-object
(as a logistic models) or a robust linear regression (as least sum of
absolute errors).
In this case, the appropriate "method" to use is "mrt"? Or another
On 09/12/2007 3:04 PM, Rina Oldager Miehs wrote:
> Hey
>
> do anyone know why this error occurs??
>
>> for(i in 1:n_trait){
> +for( j in 1:n_trait){
> +rG[i,j] <- (G_o[i,j]/(sqrt(G_o[i,i]%*%G_o[j,j]))
> +rP[i,j] <- (P_o[i,j]/(sqrt(P_o[i,i]%*%P_o[j,j]))
> Error: unexpected symbol in:
>
Dear R-list,
May I ask for help in interpretating the output of 'lmer' (from the lme4
package) when dealing with negative binomial data ?
I'm using the functions glm.nb (from the MASS package) and lmer (from the
lme4) to fit respectively fixed-effects and mixed-effects generalized linear
models
On Sunday 09 December 2007 15:09, Duncan Murdoch wrote
> [...] To get what you
> want, use match.call() to get the call in unevaluated form, and then
> process it. For example,
>
> f2 <- function( ... ) {
>call <- match.call()
>args <- as.list(call)[-1]
>unlist(lapply(args, deparse))
>
Thank you, this seems promissing. How can I address then a certain
element e.g. to add the second and the third component?
Faithfully,
Mag. Ferri Leberl
Am Sonntag, den 09.12.2007, 07:11 -0800 schrieb jim holtman:
> You program contains:
>
> de.dd<-c("Dresden",51,2,13,44)
>
> and since you are
Hi there,
I need to save a series of lattice plots as a PDF,
this is my code so far:
windows(height=8,width=6)
plot.new()
library('grid')
lattice.options(layout.heights=list(top.padding=list(x=0.15,
units="inches")))
print(plot1, split=c(1,1,2,3), more=TRUE)
print(plot2, split=c(1,2,2,3), mor
Hi:
There, I've got a question about how to read in expressions as function
parameters and it really bothered me.
I'm going to use curve() function to plot curves, and I'd like to write a
menu function to let use input math expressions.
say, if I'd like curve(3*x*x-4/x, 10, 100), I use sca
I'm trying to get Cairo and RApache to work together on OS X leopard.
Initializing Cairo stops the whole HTTP response. Apache logs show
this message:
Break on
__THE_PROCESS_HAS_FORKED_AND_YOU_CANNOT_USE_THIS_COREFOUNDATION_FUNCTIONALITY___YOU_MUST_EXEC__
() to debug.
The process has forked
Several different ways. You can look in the Intro to R to learn more:
> my.list <- list(a=matrix(1:9,3), b=1:4, c="this is a string", d=list(1,2,3))
> my.list
$a
[,1] [,2] [,3]
[1,]147
[2,]258
[3,]369
$b
[1] 1 2 3 4
$c
[1] "this is a string"
$d
$d[[1]
On Thu, 06-Dec-2007 at 07:16AM -0800, Svempa wrote:
|>
|> I want to fit a fairly long main title for a plot, supposedly by changing row
|> after a while. As for now it starts way outside the picture margin at the
|> left and continues way out right passed the right margins.
|>
|> >plot(A,main="T
On 12/9/07, Judith Flores <[EMAIL PROTECTED]> wrote:
> Hi there,
>
> I need to save a series of lattice plots as a PDF,
> this is my code so far:
>
> windows(height=8,width=6)
> plot.new()
> library('grid')
> lattice.options(layout.heights=list(top.padding=list(x=0.15,
> units="inches")))
> pri
At 3:27 PM -0500 12/9/07, Duncan Murdoch wrote:
>On 09/12/2007 3:04 PM, Rina Oldager Miehs wrote:
>> Hey
>>
>> do anyone know why this error occurs??
>>
> >> for(i in 1:n_trait){
>> +for( j in 1:n_trait){
>> +rG[i,j] <- (G_o[i,j]/(sqrt(G_o[i,i]%*%G_o[j,j]))
>> +rP[i,j] <- (P_o[i,j
Hello,
I'm using lattice to create a multi-panel figure. I would like
to draw each panel's y-axis ticks and labels on the right-hand
side of the panel. Ordinarily, I would do this by specifying
scales=list(y=list(draw=T, alternating=2)). But in this case, I am
using relation="sliced" t
Hi Deepayan and everyone,
I need to add a common legend to a group of latice
graphs, I have tried different ways using viewport and
grid.draw without success.
Here is what I have:
plot.new()
library(grid)
library('IDPmisc')
print(plot1, split=c(1,1,2,4), more=TRUE)
print(plot4, split=c(2,1,
Hi, There are two unanswered reports of this error (below) in read.spss()
when used with a .por file. I had the same problem and in order to
successfully read the file into R I downloaded spss, saved the data as
.dat tab-delimited, and then used read.table(). The point is that spss
successfully
On Dec 9, 2007 6:58 PM, Judith Flores <[EMAIL PROTECTED]> wrote:
> Hi Deepayan and everyone,
>
>I need to add a common legend to a group of latice
> graphs, I have tried different ways using viewport and
> grid.draw without success.
Try looking at the first example in
http://dsarkar.fhcrc.org
Thanks in advance - am looking for a Tutorial for doing basic stats. I have
already looked/looking at the R-intro.pdf at the R site.
Regards
BK
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://stat.ethz.ch/ma
I'm having no luck getting GARCH parameter estimations. It seems simple
enough, but I don't know what I'm doing. I'm a newbie both at R and GARCH
models, so whatever is going wrong, it's probably very basic. Here's what I
do:
1. I first load the tseries package with:
library("tseries")
2.
Dear all,
I am encountering a cyclic dependency error when running R CMD check on an R
package I wrote (R version 2.6.1, Mac OS X 10.4), see the error message below.
Creating a new generic function for "print" in "clValid"
Creating a new generic function for "summary" in "clValid"
Creating
Charles C. Berry wrote:
> On Sat, 8 Dec 2007, Charles C. Berry wrote:
>
>> Gad Abraham ms.unimelb.edu.au> writes:
>>
>>> Hi,
>>>
>>> I'm using natural cubic splines from splines::ns() in survival
>>> regression (regressing inter-arrival times of patients to a queue on
>>> queue size). The queue s
You are probably using Windows: you didn't say. If so, this is a known
bug in the tseries package: it uses Fortran I/O and that fails in Rgui.
The fix (until the package is updated) is to use Rterm instead.
On Sun, 9 Dec 2007, caffeine wrote:
>
> I'm having no luck getting GARCH parameter esti
Dear helper,
After some cleaning, I found out that I cannot open my R2.6.0 console
with message: "Fatal error: uanble to retore saved data in .Rdata",
and every time I retry it just kick me off. I am wondering if I still
be able to retrieve the data that I've been working for a month.
I have inst
On 09/12/2007, Frank E Harrell Jr <[EMAIL PROTECTED]> wrote:
> Dieter Menne wrote:
>
Thank you both!
Best wishes,
Mark
--
Dr. Mark Wardle
Specialist registrar, Neurology
Cardiff, UK
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mail
Hi all,
I am meetting one problem.I am estimating the variability of a
population using one sample.I must do two-side and one-side hypothesis
test,and estimate the confidence interval,But I don't know which
function I can use !
Thank you !
__
R
58 matches
Mail list logo
