Dear Sebastian
plot(1:10, 1:10)
text(4, 9, expression(paste("<", k, ">")))
should work here.
Best regards,
Christoph
--
Credit and Surety PML study: visit our web page www.cs-pml.org
-
Of course that should have been differences in the log-likelihoods in my
previous post. Aaargh.
Simon.
--
Simon Blomberg, B.Sc.(Hons.), Ph.D, M.App.Stat.
Centre for Resource and Environmental Studies
The Australian National University
Canberra ACT 0200
Australia
T: +61 2 6125 7800 email: Simon.Bl
Remember that -2 * the difference in the likelihoods between the two models is
asymptotically chi-squared distributed, with degrees of freedom equal to the
difference in number of parameters between the models. So you can just
calculate that for your preferred and null models, then use the pchisq f
You can use read.zoo in the zoo package to read in the data
and then see:
https://www.stat.math.ethz.ch/pipermail/r-help/2006-December/122742.html
See ?axis for creating additional axes with classic graphics and
library(lattice)
?panel.axis
in lattice graphics. Search the archives for examples
Dear all R users,
Suppose I have a data set like this:
date price
1-Jan-02 4.8803747
2-Jan-02 4.8798430
3-Jan-02 4.8840133
4-Jan-02 4.8803747
5-Jan-02 4.8749683
6-Jan-02 4.8754263
7-Jan-02 4.8746628
8-Jan-02 4.8753500
9-Jan-02 4.8882416
10-Jan-02
Hi all,
I met a problem while using the debug package, I have the following
program:
mainfun<- function(){
beta<-1
result<-subfun(beta+x)
}
subfun<-function(expr){
y <- eval(expr, envir=list(x=c(1,2)),enclos = parent.frame())
packages gstat and geoR both have kriging functions. There are probably others.
Have a look at the spatial task view on CRAN.
HTH,
Simon.
[EMAIL PROTECTED] wrote:
> Dear R-list members,
>
> I wish everyone a happy and successful 2007!
>
> Does anyone know of R-based software for
> optimal spatia
Dear R-list members,
I wish everyone a happy and successful 2007!
Does anyone know of R-based software for
optimal spatial prediction (kriging)?
We are working on a seismic event characterisation
technique and need to do some kriging.
Any help would be greatly appreciated.
Augusto
I've just bought a couple of "iButton Thermochrons" (logging
thermometers), and I'd like to access them through their Java interface
from R. But I've never really used Java, so I'm running into a problem,
and I hope there's a very simple solution.
I've managed to use rJava to create an object
Hi all,
I am hoping someone can help with a problem I have. I want to do a
zero-inflated negative binomial model on some count data. I have found
how to get the model (using zicounts), and the test of each predictor on
both the negative binomial and zero-inflated parts of the distribution.
C
On 1/4/07, Seth Falcon <[EMAIL PROTECTED]> wrote:
Prof Brian Ripley <[EMAIL PROTECTED]> writes:
> [I am not sure who is actually maintaining RSQLite, so am Cc: both the
> stated maintainer and the person who prepared the package for
> distribution. The posting guide asked you to contact the maint
To Jeff: a tip of the hat.
I have another question after Jeff's solution:
On Jan 3, 2007, at 5:09 PM, Jeffrey Horner wrote:
> Michael Kubovy virginia.edu> writes:
>> I tried and it gave a strange result. See
>> http://people.virginia.edu/~mk9y/mySite/twoGaussian.R
>> and
>> http://people.virgin
On 1/3/2007 5:30 PM, [EMAIL PROTECTED] wrote:
> Hello list,
>
> I have encountered a problem trying to save graphs using the R-graphics
> menu: File|Save as. The menu suggests that files may be saved as either
> Metafile, Postscript, pdf, png, bmp, jpeg.
> When I specify any of those file format
Hello list,
I have encountered a problem trying to save graphs using the R-graphics
menu: File|Save as. The menu suggests that files may be saved as either
Metafile, Postscript, pdf, png, bmp, jpeg.
When I specify any of those file formats a menu comes up requesting a file
name. After providing
Sorry, Bert. I didn't notice your use of "apply", which will indeed give
you separate regression results using all available data. But I was
wondering, if there was a way to modify "lm" to be able to accomplish this,
since it is doing separate regressions anyway.
Ravi.
-
Hi Bert.
Thank you so much, your solution with "apply" works perfectly. Sorry, I
know this was an elementary question, and I saw the statement you referred
to on the Help page. I just wasn't sure why, considering that there is a
facility for na options, the option of treating the dependent va
No, Bert, lm doesn't produce a list each of whose components is a separate
fit using "all" the nonmissing data in the column. It is true that the
regressions are independently performed, but when the response matrix is
passed from "lm" on to "lm.fit", only the complete rows are passed, i.e.
rows w
Ravi:
You misinterpreted my reply -- perhaps I was unclear. I did **not** say that
lm() with a matrix response would do it, but that the apply construction or
an explicit loop would. As you and the poster noted, lm() produces a
separate fit to each column of only the rowwise complete data.
Bert
Dieter Menne wrote:
> Michael Kubovy virginia.edu> writes:
>
>> I tried and it gave a strange result. See
>> http://people.virginia.edu/~mk9y/mySite/twoGaussian.R
>> and
>> http://people.virginia.edu/~mk9y/mySite/twoGaussian.pdf
>>
>> *
>> Session Info
>> *
Thanks all for your hints and extensive codes. With list it seemed to work.
lars
Gabor Grothendieck wrote:
>
> Or you could define it as your own class and define your own print and
> other methods, e.g.
>
>> X <- structure(array(1:8, c(2,2,2)), class = "twomats")
>> attr(X, "DIMNAMES") <- li
As the Help page says:
If response is a matrix a linear model is fitted separately by least-squares
to each column of the matrix
So there's nothing hidden going on "behind the scenes," and
apply(cbind(y1,y2),2,function(z)lm(z~x)) (or an explicit loop, of course)
will produce a list each of whose
On 1/3/07, Derek Eder <[EMAIL PROTECTED]> wrote:
> My question is so basic that I am (almost too) embarrassed to admit that
> I could not find an answer after an hour's worth of homework.
>
> What is the Trellis / Lattice analog for the axis(graphics) function
> that enables the creation of axes in
Hi.
I am running regressions of several dependent variables using the same set
of independent variables. The independent variable values are complete, but
each dependent variable has some missing values for some observations; by
default, lm(y1~x) will carry out the regressions using only the o
Hello,
thanks for your help.
I tried with index in database (primary keys and index) but nothing
changed.
My hard disk has two partition FAT 32: in the first there is ubuntu 6.10, in
the second (19 GB with 6GB of free space, never defragmented) there is
Windows XP HE. I have the problem under Wind
I wrote my own ewma function to deal with the somewhat odd way that
filter handles missing values.
The function I wrote works as long as the NA isn't first but when it is
first I still get a zero in the output.
I'm not expert enough to look at filter and undeerstand what it is
doing.
# 1) THE FIR
Hi folks,
I have assumed that ratios of variance components (Fst and Qst in
population genetics) could be estimated using the output of mcmcsamp
(the series on mcmc sample estimates of variance components).
What I have started to do is to use the matrix output that included
the log(variances
Hi,
I'm trying to compare models, one of which has all parameters fixed
using offsets. The log-likelihoods seem reasonble in all cases except
the model in which there are no free parameters (model3 in the toy
example below). Any help would be appreciated.
Cheers,
Jarrod
x<-rnorm(100)
y<-rn
Prof Brian Ripley <[EMAIL PROTECTED]> writes:
> [I am not sure who is actually maintaining RSQLite, so am Cc: both the
> stated maintainer and the person who prepared the package for
> distribution. The posting guide asked you to contact the maintainer:
> what response did _you_ get?]
For the reco
On Wed, 3 Jan 2007, Thomas Lumley wrote:
> On Wed, 3 Jan 2007, ONKELINX, Thierry wrote:
>
>> This is due to the internal representation of 0.1, which is not exactly
>> 0.1 but very close to it. If you want to do an integer divide, you
>> should only use integers to divide with.
>
> This must be mo
Dear All,
My name is José Cortiñas Abrahantes, I am statistician and work at the
university in Belgium. I started working recently with machine learning
techniques and I finding a fascinating field. The reason of my email is to
ask you a question related to regression forest. I am interested to co
Hi Lars,
in a 3-dim array you have 3 axes, day x, y and z (which I named
xNo/xYes, yNo/yYes and zSmall/zBig)
to assign directly the labels use :
> tab <- array(1:8, c(2, 2, 2),dimnames= list(c("xNo","xYes"),
c("yNo","yYes"),c("zBig","zSmall")))
> tab
, , zBig
yNo yYes
xNo13
xYes
Or you could define it as your own class and define your own print and
other methods, e.g.
> X <- structure(array(1:8, c(2,2,2)), class = "twomats")
> attr(X, "DIMNAMES") <- list(list(c("No", "Yes"), c("No", "Yes")),
+dimnames = list(c("No", "Yes"), c("big", "small")))
>
> print.twomats <- fun
You can't do that. If you want to have different labels
on the first two dimensions, then a 3-dimensional array
doesn't seem to be the natural data structure.
I would suggest two matrices held in a list.
Patrick Burns
[EMAIL PROTECTED]
+44 (0)20 8525 0696
http://www.burns-stat.com
(home of S Poe
Hi Wolfgang,
thanks for your hint. But I am desperating. I have an 3 dim array of say 10
matrices where every matrix has to stick an other label. I found no way to
direct assign the labels.
for example
tab <- array(1:8, c(2, 2, 2))
dimnames(tab[,,1]) <- list(c("No","Yes"), c("No","Yes"))
dimna
Matthew,
You don't seem to say what linux release you are using, They can't very well
help you without that information. Not all releases are equal.
JWD
On Wednesday 03 January 2007 03:49, Matthew Dowle wrote:
> Hi All,
>
> I'm getting the following error, could anyone help please?
>
> $ R C
try this :
x <-
array(1:24,dim=c(2,3,4),dimnames=list(letters[1:2],LETTERS[1:3],letters[23:26]))
Cheers,
Wolfgang
downunder03 a écrit :
> hi all. how can i adress a array directly. for example i wanna give array 1
> other labels than array 2. How can I overcome this problem?
>
> ...this doesn't
Dear all,
I'm trying to manage with user defined split function in rpart
(file rpart\tests\usersplits.R in
http://cran.r-project.org/src/contrib/rpart_3.1-34.tar.gz - see bottom of
the email).
Suppose to have the following data.frame (note that x's values are already
sorted)
> D
y x
1 7 0.428
2
RSQLite can import data from a large file directly (via
"dbWriteTable"). This future is quite appealing.
On 1/3/07, Prof Brian Ripley <[EMAIL PROTECTED]> wrote:
> I guess you are using package RSQLite without telling us (or telling us
> the version), and that your example is incomplete?
>
> Using
If you just want to label and identify outliers after creating a plot
then look at the identify function. In your case you could just run the
following command after creating your plot:
> out.index <- identify( c$lb, c$index )
Then click on (or near) the outliers or other interesting points. Th
On Sat, 30 Dec 2006, Gabor Grothendieck wrote:
> That has two disadvantages:
>
> (1) it only works if the user is defining ll himself; however, if the
> user is getting
> ll from somewhere else then its not applicable since the user no
> longer controls its
> scope whereas resetting the environmen
On Sat, 30 Dec 2006, Sebastian P. Luque wrote:
> On Sat, 30 Dec 2006 15:46:01 -0600 (CST),
> Luke Tierney <[EMAIL PROTECTED]> wrote:
>
>> It is much cleaner to do this sort of thing with lexical scope. For
>> example,
>
>> mkll <- function(x, y) {
>>function(ymax=15, xhalf=6) {
>>
On Wed, 3 Jan 2007, ONKELINX, Thierry wrote:
> This is due to the internal representation of 0.1, which is not exactly
> 0.1 but very close to it. If you want to do an integer divide, you
> should only use integers to divide with.
This must be more-or-less correct, but it is worth noting that
> 0
Hello everyone!
I'm trying to plot some mathematical expression along my axis, but
demo(plotmath) did not have the symbol I was looking for. In particular,
I would like to denote the mean of an observable by writing
which I tried to enter with
expression(group("<", k, ">"))
However, my naive
Try the following:
act <- c('good','good','bad','bad','good','good','bad','bad')
pred <- c('good','bad','bad','bad','good','good','good','bad')
table(pred,act)
table(pred,act)/apply(table(pred,act),1,sum)
Cheers,
Andreas
On 1/3/07, Feng Qiu <[EMAIL PROTECTED]> wrote:
> Hi everybody, I'm trying
Try this:
> actual <- factor(c("good", "good", "bad", "bad", "good", "good", "bad",
> "bad"))
> pred <- factor(c("good", "bad", "bad", "bad", "good", "good", "good", "bad"))
> table(actual, pred)
pred
actual bad good
bad31
good 13
> prop.table(table(actual, pred), 1)
Hi everybody, I'm trying to do a statistic on the error rate of a prediction
algorithm.
suppose this is the real category
[good, good, bad, bad, good, good, bad, bad]
this is the predicted category
[good, bad, bad, bad, good, good, good, bad]
I'm trying to do a statistic on the error rate for ea
Thierry Onkelinx wrote:
> If you want to do an integer divide, you should only use integers to
> divide with.
I think this should go into ``fortunes''.
cheers,
Rolf Turner
[E
hi all. how can i adress a array directly. for example i wanna give array 1
other labels than array 2. How can I overcome this problem?
...this doesn't work
tab <- array(1:8, c(2, 2, 2))
dimnames(tab[,,1]) <- list(c("No","Yes"), c("No","Yes"),c("ARRAY1"))
dimnames(tab[,,2]) <- list(c("big","smal
Without any actual example Ias requested in the footer of this message) I
can only guess, but the most common cause of slow queries is the lack of
indices in the database, so did you create any?
You haven't told us your actual OS (beyond 'Windows'), but a guess is that
your processes are I/O bo
Dear List,
when reading MS Excel files in R using package RODBC I encountered the
problem of having the first line of data getting omitted.
I read the data as :
> library(RODBC)
> channel1 <- odbcConnectExcel("myFile.xls")
> sheet1 <- sqlQuery(channel1, "SELECT * FROM [Cell measures (1)$]")
Nevermind the CPU usage, the likely problem is that your queries are
inefficient in one or more ways (i.e., you don't use indexes when you
really should - it's impossible to guess without knowing how the data
and the queries look like, which somehow you've decided are not
important enough to descri
This is due to the internal representation of 0.1, which is not exactly
0.1 but very close to it. If you want to do an integer divide, you
should only use integers to divide with.
Cheers,
Thierry
ir. Thierry Onkelinx
I am confused about why the following occurs:
> version
_
platform i386-pc-mingw32
arch i386
os mingw32
system i386, mingw32
status
Hello everyone,
I have a problem when I execute queries using R 2.3.1 and MySql server 5.0.
What I do: I load data in different csv files (every file represents a
particular temporal step of a simulation) using Mysql query "load data" with
RMySQL command DbSendQuery (but the same problem there is
This should do it. You can play around with the alignment of the labels.
input <- 'ID index lb
100 FLINDYTHNIPLI 1.84770221 9.087463
101 none 0.06657547 8.927778
102 GDDKVYSANGFTT -0.22922544 8.599913
103 GDFTQGPQSAKTR 0.01203925 8.483816
104 GDKEFSDALGYL
Hi
On 3 Jan 2007 at 12:54, Jenny persson wrote:
Date sent: Wed, 3 Jan 2007 12:54:50 +0100 (CET)
From: Jenny persson <[EMAIL PROTECTED]>
To: r-help@stat.math.ethz.ch
Subject:[R] How to add characters on graph ?
> Dear R-users,
>
This is not a 'compile error' but an error from the configure script,
probably a run error. As it says
> See `config.log' for more details.
and then ask your local Linux guru what is broken locally.
This is neither an R nor an RODBC problem, and for what it is worth
RODBC_1.1-7.tar.gz configu
On Wed, 3 Jan 2007, Atte Tenkanen wrote:
[...]
> By the way, right now the demo(Hershey) seems not to work in OSX version R
> 2.4.1. ...
> I get a message
>
>> i <- i + 1
> Error in deparse(ei, control = c("showAttributes", "useSource")) :
> invalid multibyte string
But it does should wo
Dear R-users,
I have following data
# Plot coloured scatter plot
c<-dat[100:110,c(5,7,8)]
par(mfrow=c(3,2))
plot(c$lb,c$index, pch=1, col=5,cex=1, lwd=2,
xlab="LB", ylab="Index",cex.main =1,font.main= 1,
main="scatterplot")
ID index
The function needs to have a single parameter. Then extract each
parameter.
For example, or see first example in help for optim.
errorFunction <- function(params) {
localShifts <- params[1]
etc
>>> "Gerster Sarah" <[EMAIL PROTECTED]> 01/03/07 9:22 PM >>>
Hi!
I'm trying to figure out how to us
Hi,
I am interested in integrating R with a Java front end. Before deciding to
use R I am concerned about multi-threading.
I have been investigating R's capability of dealing with multiple requests
simultaneously (multi-threading)
and have looked at Luke Tierney's 2001 notes for ideas for futur
Hi Derek,
see ?xyplot and ?panel.axis
Hint: RSiteSearch("panel.axis") will point you to examples.
Thomas
Derek Eder wrote:
> My question is so basic that I am (almost too) embarrassed to admit that
> I could not find an answer after an hour's worth of homework.
>
> What is the Trellis / Latti
Hi All,
I'm getting the following error, could anyone help please?
$ R CMD INSTALL RODBC_1.1-7.tar.gz
* Installing *source* package 'RODBC' ...
checking for gcc... gcc -std=gnu99
checking for C compiler default output file name... a.out
checking whether the C compiler works... yes
checking whet
Hi,
I am considering using R to integrate with a Java application. However,
before deciding upon R I need to understand if R is capable of dealing with
multiple requests simulataneously.
Is a single instance of R capable of dealing with multiple simulataneous
requests or does a new instance of
My question is so basic that I am (almost too) embarrassed to admit that
I could not find an answer after an hour's worth of homework.
What is the Trellis / Lattice analog for the axis(graphics) function
that enables the creation of axes in locations other than the default
(i.e., bottom for X a
Hi!
I'm trying to figure out how to use optim... I get some really strange results,
so I guess I got something wrong.
I defined the following function which should be minimized:
errorFunction <-
function(localShifts,globalShift,fileName,experimentalPI,lambda)
{
lambda <- 1/sqrt(147)
# error
Hello Christophe,
Thanks a lot! This is what I need. My purpose is to generate chords and I need
interactive responses straight in R. I can output midi event lists as csv-files
and convert them with a nice midicsv-program in linux- or OSX-console. Then it
is easy to convert midi files to notes
>>> Prof Brian Ripley<[EMAIL PROTECTED]> 3/1/2007 09:22 >>>
>In this particular case 'An Introduction to R' has a comprehensive
>description of graphical parameters with figures (as do all good books on
>S/R e.g. MASS4 - since it has the same first author).
Thanks, Brian,
I've reached both "
Hi, Atte
> De : [EMAIL PROTECTED]
> [mailto:[EMAIL PROTECTED] De la part de Atte Tenkanen
> Envoyé : mercredi 3 janvier 2007 09:17
> Hi,
>
> I'd like to know if it is possible to use Hershey vector
> fonts to create very primitive musical notation.
> [...]
There is an example of a music scor
In this particular case 'An Introduction to R' has a comprehensive
description of graphical parameters with figures (as do all good books on
S/R e.g. MASS4 - since it has the same first author).
On Tue, 2 Jan 2007, Ricardo Rodríguez - Your XEN ICT Team wrote:
--
Ricardo Rodríguez
Your XEN IC
Hi,
I'd like to know if it is possible to use Hershey vector fonts to create very
primitive musical notation.
If I can hang some whole notes on these lines
X11()
plot(0,0, xlim=c(0,10), ylim=c(0,10))
# Staves:
for (i in c(seq(from=2,to=2.8,by=0.2),seq(from=4,to=4.8,by=0.2)))
{
abline(h=
On Tue, 2 Jan 2007, [EMAIL PROTECTED] wrote:
> Hello,
>
> I am building R-2.4.1 on an SGI IRIX platform, using gcc 3.3.
> gmake check failed, and the arith-true.Rout.fail file indicated:
>> is.na(mean(c(1,NA,NA)[-1], trim = .1, na.rm = TRUE))
> [1] FALSE
>>
>
> I tried the mean() command in R and
I guess you are using package RSQLite without telling us (or telling us
the version), and that your example is incomplete?
Using RSiteSearch("RSQLite Windows") quickly shows that this is a
previously reported problem with the package, e.g.:
http://finzi.psych.upenn.edu/R/Rhelp02a/archive/7251
73 matches
Mail list logo