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 one of
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 is 0.
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: Reached
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
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,
Mark Wardle mark at 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
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 reverse
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 = special,encoding =
Dieter Menne wrote:
Mark Wardle mark at 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
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
f1(
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 best
Dear Everybody,
Enclosed are the tiny beginning of a program and its output.
As a pity, if I load the image load(kanal.RData) the elements of
de.dd are non numerical; de.dd[2]+de.dd[3] returns Fehler in
de.dd[2] + de.dd[3] : nicht-numerisches Argument für binären Operator.
How can I keep the
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 51
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 load(kanal.RData) the elements of
de.dd are non numerical;
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 Phy Conn
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?
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 an
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 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 Pattern
- 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 with length(x)
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
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 to me
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 lm(Y~X)
by
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' of X from
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,j]
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
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
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 mixing
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),
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
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
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]]
[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=This is
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)))
print(plot1,
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,
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
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
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
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 a new
Charles C. Berry wrote:
On Sat, 8 Dec 2007, Charles C. Berry wrote:
Gad Abraham g.abraham at 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 size
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
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
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
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 !
__
51 matches
Mail list logo