Edna Bell gmail.com> writes:
> I have a tiny data set:
> > xs
> [1] 0.7 2.8 0.1 1.9 0.0 1.4 0.2 2.3 0.3 0.2
> >
>
> and generate the density function for it.
>
> I would like to replicate this. Is there a straightforward way to do
> this, please?
See chapter 5 in MASS and the examples in lib
Hi,
i have a data with values against time and date. for plotting against date
i used axis.Date() option. but I need a 24hr time plot as well. i have
extracted and plotted, but not being able to put the time ticks ( eg. 6am, 7
am etc) on the x-axis. Im new to R, any help would be appreciated!!
-
Wow, i had no idea R was capable of producing graphics like that,
smegging awesome!
On 20 Feb, 19:28, Duncan Murdoch wrote:
> On 2/20/2009 1:55 PM, Duncan Murdoch wrote:
>
>
>
> > On 2/20/2009 1:46 PM, kapo coulibaly wrote:
> >> Ideally I would want it to look like a rubik cube with each little c
Thanks yet another time, Gabor -
I think I am slowly understanding - particularly I was confused by persistence
of connections.
So starting with some parts of your example 12,
##
gc()
rm(list=ls())
unlink("mydb")
sqldf("attach 'mydb' as new")
irishead <- file("irishead.dat")
iristail <- file("
Hi listers...
I am trying to do a matrix multiplication...
Suppose my data has 100 obs and 2 variables...
I calculated the estimateurs of location and dispersion...
So now I would like to do the following...
t(Xi-location)%*%solve(dispersion)%*%(Xi-location)
My code is...
Z<-rep(2,200)
Z<-matrix(Z
> wt.filter("d8")
Error in validObject(.Object) :
invalid class "wt.filter" object: invalid object for slot "transform" in
class "wt.filter": got class "function", should be or extend class "character"
How come ? The above function works fine with the other Daubechies filters.
Thank you.
Best
Hi R Gurus:
I have a tiny data set:
> xs
[1] 0.7 2.8 0.1 1.9 0.0 1.4 0.2 2.3 0.3 0.2
>
and generate the density function for it.
I would like to replicate this. Is there a straightforward way to do
this, please?
According to the help file, the FFT is used.
How does this compare to Silverman'
On 2/20/09, Fredrik Karlsson wrote:
> Dear list,
>
> I am sorry for asking you this, but I am trying to do again what I
> thought I have done before, although this time it does not work.
>
> So, given the data set:
>
>> testdf <- data.frame(grfak=sample(c("One","Two"),size=100,replace=TRUE),
>> pa
The dev.new function starts the default graphics device for that platform, so
is the device independent method for doing this. Individual users can change
the default using options().
Hope this helps,
--
Gregory (Greg) L. Snow Ph.D.
Statistical Data Center
Intermountain Healthcare
greg.s...@i
I have the libx11-dev package installed through the aptitude
application. This is on 5.0 power pc. What can I provide to help
with this? I don't know where to begin.
--
Stephen Sefick
Let's not spend our time and resources thinking about things that are
so little or so large that all they rea
thanks you both for the reply.
yes, the second question is two-in-one, a bit confusing: I was thinking
about Matlab's hold
command, but also waiting for user input to proceed. So you have provided
the complete answer.
The first question though - yes, windows() will give you a different window
on
At this point, the function would only return the names because unless
specified otherwise they return the result of the last evaluation. Add
another line with just "a" and assign the results of the function to
something:
> testold=function(){
+ assign("a",c(1,2,3),env=.GlobalEnv)
+ names(a
You might take a look at this competitor to conv2:
http://www.advancedsourcecode.com/ffw.asp
... which is Open Source, not too long, and appears to be built out
of tools which have analogs in the signal package:
http://cran.r-project.org/web/packages/signal/signal.pdf
--
David winsemius
On Fri, 20 Feb 2009, GAF wrote:
>
> Hi there,
> I am trying to measure orthogonal/perpendicular distances from regression
> lines (i.e. the shortest distance from a given point to my regression line).
> As it sounds rather easy (basically just an orthogonal/perpendicular
> residual) I hoped that t
Of interest to R Users, here are notes from a panel Jim Porzak and I
co-hosted this week -- how two companies in the Bay Area are using R.
http://dataspora.com/blog/predictive-analytics-using-r/
"Facebook's Data Team used R in 2007 to answer two questions about new
users: (i) which data points p
On Feb 20, 2009, at 4:53 AM, Simon Pickett wrote:
Hi Oliver,
1) thats easy, if you want to display several graphs at once type
windows(), once for each new graph.
2) Not sure exactly what you mean here but you have complete control
over graphs in R. e.g. plot a blank graph then add axex
Hi David,
This was useful, thanks.
The example was just that, there are no "a"s, "b",s etc I was
just trying to show what I was trying to do.
Using m[v,]
where m was the matrix in my example and v the vector of
index values works great as you suggested Thanks.
R is quite a powerful language
> "KK" == Ken Knoblauch
> on Fri, 20 Feb 2009 14:11:55 + (UTC) writes:
KK> Thomas Lumley u.washington.edu> writes:
>> I want to construct a symmetric band matrix in the Matrix
>> package from a matrix where the first column contains
>> data for the main diagonal,
Since you do not offer reproducible code, I must guess that you were
thinking of this as a matrix with 6 instances of "a" across the top
row. See if this helps:
> mtx<-matrix( rep(letters[1:4], 6), nrow=4)
> mtx
[,1] [,2] [,3] [,4] [,5] [,6]
[1,] "a" "a" "a" "a" "a" "a"
[2,] "b" "b
Dear Peter,
Yes, you're absolutely right -- of course (though I didn't notice it!), the
precedence of / and * in a formula is like division and multiplication in an
arithmetic expression. With this correction, the SSs and error SSs are
properly matched by aov(), producing the same F-tests as Anova
If I assign a variable in a function, as in:
test=function(){
assign("a",c(1,2,3),env=.GlobalEnv)
}
How can I do the equivalent of:
names(a)=c("one","two","three")
within the function?
Merely adding the call to names does not work within the function as
it only affects a local variable:
t
Hello
I have a matrix of size rows x cols.
I also have a vector of size rows.
The vector contains index values that corresponds to rows in the matrix.
I would like to re-arrange/sort the contents of the matrix according
to the entries in the vector. Can this be done efficiently in R and
possib
Hello John, thanks for your reply and correction.
I apologies for my crude mistake in applying the aov (now I have
learned better). I hope to get a hold of "Statistical Models in S",
but I don't predict it could easily happen in the near future.
Also, I would be very happy if you could supply me w
John Fox wrote:
Dear Tal,
I didn't have time to look at all this yesterday.
Since aov() doesn't do what I typically want to do, I guess I've not paid
much attention to it recently. I can see, however, that you appear to have
specified the error strata incorrectly, since (given your desire to co
Dear Tal,
I didn't have time to look at all this yesterday.
Since aov() doesn't do what I typically want to do, I guess I've not paid
much attention to it recently. I can see, however, that you appear to have
specified the error strata incorrectly, since (given your desire to compare
to Anova) th
Rolf Turner wrote:
>
> On 21/02/2009, at 12:54 AM, lauramorg...@bluewin.ch wrote:
>
>> Thank you for your advice, but I didn't manage to make it work...
>> I tried
>>
>> carichi.annui <-
>> data.frame(anno,loadPTG,loadPO4,loadNT,loadNH4,loadNO3,loadBOD5,loadSiO2)
>>
>>
>> And I got this error messa
Dear list,
I am sorry for asking you this, but I am trying to do again what I
thought I have done before, although this time it does not work.
So, given the data set:
> testdf <- data.frame(grfak=sample(c("One","Two"),size=100,replace=TRUE),
> panfak= sample(c("Yes","No"),size=100,replace=TRUE)
On 21/02/2009, at 12:54 AM, lauramorg...@bluewin.ch wrote:
Thank you for your advice, but I didn't manage to make it work...
I tried
carichi.annui <- data.frame
(anno,loadPTG,loadPO4,loadNT,loadNH4,loadNO3,loadBOD5,loadSiO2)
And I got this error message:
Error in data.frame(anno, loadPTG,
Hi,
I am trying to test the random labeling hypothesis with my data. I have a
point pattern of bird nests and each nest has a mark that it is either predated
or not-predated. I want to test the null hypothesis that the K distribution of
predation events is equal to the K distribution of all
On 2/20/2009 1:55 PM, Duncan Murdoch wrote:
On 2/20/2009 1:46 PM, kapo coulibaly wrote:
Ideally I would want it to look like a rubik cube with each little cube
color coded based on the fourth column (data column). Your suggestion might
work if I could color code based on data in the fourth colum
For the bivariate case:
g <- function(b,x,y)
(abs(y - b[1] - b[2] * x))/sqrt(1 + crossprod(b))
url:www.econ.uiuc.edu/~rogerRoger Koenker
emailrkoen...@uiuc.eduDepartment of Economics
vox: 217-333-4558University of Illinois
fax: 2
Hi all,
I am using the function COV.WT to estimate the estimators (location and
scale) of a bivariate cauchy distribution.
My doubt is about the option WT (weight), cause at the R-help shows that the
weight is uniform according to the number of observations. But, checking the
theory, for example,
Hi,
I have a large data set in binary code, no covariates. Say, positions along a
genomic sequence where reference sequence is represented by 0 and changes
represented by 1. I have 99 positions and 2000 sequences to analyze. I want to
run a univariate analysis to isolate positions where majorit
Hi,
Storer and Crowley (JASA 1985) presented an approach for approximating the
changes in maximum partial-likelihood parameter estimates for the Cox model
when a single observation is deleted. Is there an R implementation of this
approach?
Any help is greatly appreciated. Thanks.
Best,
Ravi
Hi there,
I am trying to measure orthogonal/perpendicular distances from regression
lines (i.e. the shortest distance from a given point to my regression line).
As it sounds rather easy (basically just an orthogonal/perpendicular
residual) I hoped that there was some function in R that can do that
On 2/20/2009 1:46 PM, kapo coulibaly wrote:
Ideally I would want it to look like a rubik cube with each little cube
color coded based on the fourth column (data column). Your suggestion might
work if I could color code based on data in the fourth column.
Thanks
There's no primitive "cube" symbo
Ideally I would want it to look like a rubik cube with each little cube
color coded based on the fourth column (data column). Your suggestion might
work if I could color code based on data in the fourth column.
Thanks
On Fri, Feb 20, 2009 at 10:56 AM, Duncan Murdoch wrote:
> On 2/20/2009 10:50 AM
baptiste auguie wrote:
> Hi,
>
> something like this perhaps,
>
>
> create_string <- function(.s){
> result <- read.table(textConnection(.s))
> sapply(result, as.character)
> }
>
> (test <- create_string("ab cd ef"))
or:
unlist(strsplit('ab cd ef', split=' '))
vQ
>
> hope this helps
>
Imagine the following data of noises that people make, and you want to know
who makes the noise "blub."
data=data.frame(rbind(c("blah, boo, blub"),c("blabber, baeh,
blib"),c("quack, cheer, blub")))
names(data)=c("noises")
data
searchTerm="blub"
#The first gives a list of all noises that people m
Dear R-users,
I would have a question regarding the multiplication between a
matrix (A, mXn) and a block-diagonal matrix, B, (m*n)X(m*n).
The easy solution would be to, first, vectorize A and use simple matrix
multiplication, but my aim is to avoid such vectorization. Is it
possible?
Pleas
Still searching what's going on In fact, I can download manually any
package and install it. The problem occurs when updating the html files. I
can reproduce it with this:
> link.html.help()
Error in .readRDS(pfile) : unknown input format
> traceback()
5: .readRDS(pfile)
4: .packages(all.avail
par(ask=TRUE)
Uwe Ligges
Will Stone wrote:
Hi,
How do you set up the graphics window so that when clicked by the mouse; it
cycles through your different plots?
ie, run a function with multiple graphs in it: one graphics window pops up
and clicking it goes to the next graph.
I have seen this
Hi Sergii,
>> I have trouble obtaining the same results for nested Anova with two fixed
>> factors when using
>> lm and aov functions.
There is no difference between the two if you treat them equally, i.e. if
you summarize them in the same way.
## Try:
anova(e2)
summary(e1)
## Or:
summary.lm(
On Fri, 20 Feb 2009, Philipp Pagel wrote:
On Fri, Feb 20, 2009 at 04:12:53PM +, Jarrod Hadfield wrote:
I was wondering whether there was an easy way to change the attributes
of all elements in a data.frame (rather than looping through elements)?
Specifically, I would like to set the "dim"
Hi,
How do you set up the graphics window so that when clicked by the mouse; it
cycles through your different plots?
ie, run a function with multiple graphs in it: one graphics window pops up
and clicking it goes to the next graph.
I have seen this done before but cannot remember how to replicate
Hi,
something like this perhaps,
create_string <- function(.s){
result <- read.table(textConnection(.s))
sapply(result, as.character)
}
(test <- create_string("ab cd ef"))
hope this helps
baptiste
On 20 Feb 2009, at 16:38, Sean Zhang wrote:
My dear R-helpers:
I am a novic
Dear all,
I got a code for e1071 package in R for SVM regression. I
have used *m$coefs* for extracting the coefficients but I am getting only
72 . How can I extract coefficients of the predictors set? Does it mean
that I will get only 72 as *Number of Support Vectors: 72. *
**
On Fri, Feb 20, 2009 at 04:12:53PM +, Jarrod Hadfield wrote:
> I was wondering whether there was an easy way to change the attributes
> of all elements in a data.frame (rather than looping through elements)?
>
> Specifically, I would like to set the "dim" attributes to NULL
Maybe I just don'
My dear R-helpers:
I am a novice in R and have the following text string manipulation question.
Is there a function that performs the job described below?
Say,
wanted_output <- c("ab", "cd", "ef")
#the function_wanted can generate c("ab", "cd", "ef") using ab cd ef as the
single input argument
w
I take the liberty of pasting the on-line documentation for package "wavelets",
functions "dwt" and "wt.filter", pertinent to the wavelet/filter choice.
While I apologize for my verbose messages with huge attachments, I encourage
whoever loads packages to CRAN to take some minutes to reread th
Right, thanks. I feel a bit sheepish about asking the question now.
I'm using the latex function in the Hmisc package to generate the
tables, which is why it wasn't obvious to me that I just needed to
insert the new line myself.
Thanks,
Ista
On Fri, Feb 20, 2009 at 11:24 AM, Duncan Murdoch wrote:
On 2/20/2009 11:08 AM, Ista Zahn wrote:
Dear Professor Leisch + R helpers,
For the past few days I have been trying to figure out why the LaTeX
endfloat package is not working with Sweave. I figured it out, and
it's trivial: The endfloat package requires that \begin{table} and
\end{table} command
Try this:
c4 <- cutree(cluster, k=4)
by(data, c4, mean)
HTH
Marcelino
On 2009-02-20 Jgaspard wrote:
Hi all!
I'm new to R and don't know many about it. Because it is free, I managed to
learn it a little bit.
Here is my problem: I did a cluster analysis on 30 observations and 16
variab
jgaspard wrote:
Hi all!
I'm new to R and don't know many about it. Because it is free, I managed to
learn it a little bit.
Here is my problem: I did a cluster analysis on 30 observations and 16
variables (monde, figaro, liberation, etc.). Here is the .txt data file:
"monde","figaro","liberat
Hi,
I was wondering whether there was an easy way to change the attributes
of all elements in a data.frame (rather than looping through elements)?
Specifically, I would like to set the "dim" attributes to NULL
Thanks for any help,
Jarrod
--
The University of Edinburgh is a charitable bod
Dear Professor Leisch + R helpers,
For the past few days I have been trying to figure out why the LaTeX
endfloat package is not working with Sweave. I figured it out, and
it's trivial: The endfloat package requires that \begin{table} and
\end{table} commands are on their own line. But Sweave is wri
You can try rgl package.
It utilizes OpenGL library.
Or, check the CRAN to find other gtk based package
- Original Message -
From: "kapo coulibaly"
To:
Sent: Friday, February 20, 2009 11:50 PM
Subject: [R] 3D or 4D plot
I have data on a regular grid in the format: x y z data. I wou
On 2/20/2009 10:50 AM, kapo coulibaly wrote:
I have data on a regular grid in the format: x y z data. I would like to
be able to plot them in 3d kind of like a volume or a mesh with colors. Is
that possible in R and how?
There are lots of ways, depending on what you want the plot to look
lik
Marc Schwarz's solution solved my problem.
Thank you Marc!
John
John David Sorkin M.D., Ph.D.
Chief, Biostatistics and Informatics
University of Maryland School of Medicine Division of Gerontology
Baltimore VA Medical Center
10 North Greene Street
GRECC (BT/18/GR)
Baltimore, MD 21201-1524
(Phone)
I have data on a regular grid in the format: x y z data. I would like to
be able to plot them in 3d kind of like a volume or a mesh with colors. Is
that possible in R and how?
Thanks
[[alternative HTML version deleted]]
__
R-help@r-project.or
Hi all!
I'm new to R and don't know many about it. Because it is free, I managed to
learn it a little bit.
Here is my problem: I did a cluster analysis on 30 observations and 16
variables (monde, figaro, liberation, etc.). Here is the .txt data file:
"monde","figaro","liberation","yespeople","n
I met today a computer crash and our maintenance officer had to reinstall
some components of the OS (MS Windows XP Pro) as well as the Internet
browser (among other things). Now, I cannot install packages from a distant
repository:
> utils:::menuInstallPkgs()
Error in .readRDS(pfile) : unknown inp
You need unixODBC or iodbc. Specifically on F10
gannet% rpm -q --whatprovides /usr/include/sql.h
unixODBC-devel-2.2.12-9.fc10.x86_64
This is in the RODBC/README (well, not the exact RPM for F10).
On Fri, 20 Feb 2009, John Sorkin wrote:
Fedora 10
R 2.8.1
I hope someone can tell me the meanin
on 02/20/2009 09:27 AM John Sorkin wrote:
> Fedora 10
> R 2.8.1
>
> I hope someone can tell me the meaning of error I received trying to install
> RODBC and how I can get around the problem.
> n.b. The error originally occurred then I was installing Rcmdr. I then tried
> to install RODBC separa
On Fri, 20 Feb 2009, Michael Friendly wrote:
In an Rweave tutorial written for possibly naive R users, I felt it necessary
to explain why table() had to be used inside with(), whereas other tools
like xtabs() had a data= argument.
I think the simple rule is: If there is a "formula", there is
Fedora 10
R 2.8.1
I hope someone can tell me the meaning of error I received trying to install
RODBC and how I can get around the problem.
n.b. The error originally occurred then I was installing Rcmdr. I then tried to
install RODBC separately and received the same error.
checking sql.h usabil
on 02/20/2009 09:11 AM Michael Friendly wrote:
> In an Rweave tutorial written for possibly naive R users, I felt it
> necessary to explain why
> table() had to be used inside with(), whereas other tools like xtabs()
> had a data= argument.
>
> with() is quite nice for such cases, but it seems an
In an Rweave tutorial written for possibly naive R users, I felt it
necessary to explain why
table() had to be used inside with(), whereas other tools like xtabs()
had a data= argument.
with() is quite nice for such cases, but it seems an unnecessary thing
to learn right off.
Before I turn thi
Here's how I do this:
The last lines of my .First() function are these:
## Run command line program if there is one
if(length(.cmd <- commandArgs(trailingOnly = TRUE)) > 0)
try(source(textConnection(.cmd), echo = T, prompt.echo = "> "))
and on my Linux path I have this csh script, called "runR"
> "TL" == Thomas Lumley
> on Fri, 20 Feb 2009 05:49:10 -0800 (PST) writes:
TL> On Fri, 20 Feb 2009, David Winsemius wrote:
>> If you can tolerate the subdiagonal and superdiagonal zero elements
being
>> populated, then perhaps this is useful. If there is a subset functio
Then? :
xxx <- data.frame(x1 =rnorm(5), offd= c(1:4,0) )
M <-Matrix(0,5,5)
rrr <- dim(M)[1]
for (rr in 1:rrr){ M[rr,rr] <- xxx$x1[rr] }
for (rr in 1:(rrr-1)){ M[rr+1,rr] <- xxx$offd[rr] ; M[rr,rr+1] <- xxx
$offd[rr]}
> M
5 x 5 sparse Matrix of class "dgCMatrix"
[1,] -1.494930 1.000 .
Thank you!!!
unlist() worked perfectly!!!
Have a nice weekend
Laura
Messaggio originale
Da: petr.pi...@precheza.cz
Data: 20.02.2009 14.58
A:
Copia:
Oggetto: Re: [R] write.table
Hi
r-help-boun...@r-project.org napsal dne 20.02.2009 12:54:41:
> Thank you for your advice, but I didn't man
I also plan to create a R website in Taiwan (locale: Traditional Chinese).
My reason is the same as you.
Welcome people in taiwan to give me suggestion.
- Original Message -
From: "UsuarioR España"
To: ;
Cc:
Sent: Friday, February 20, 2009 5:59 PM
Subject: Re: [R] everybody loves R..
Thanks for this, Frank. Quick comment (decided to put on list rather
than just send to you directly, hope that's OK). Under mixed models, you
might consider creating two smaller columns, one for the fixed effects
and another for the random effects. Under the random effects, you might
consider check
on 02/20/2009 07:23 AM Frank E Harrell Jr wrote:
> C.H. wrote:
>> Dear experts,
>>
>> I would like to know how to plot the log-minus-log plot for survival
>> analysis (to check the proportional assumption) in R.
>>
>> Using the AML example.
>>
>> fit <- survfit(Surv(time, status) ~ x, data=aml)
>>
Thomas Lumley u.washington.edu> writes:
> I want to construct a symmetric band matrix in the Matrix
> package from a matrix where the first column
> contains data for the main diagonal, the second column
> has data for the first subdiagonal/superdiagonal
> and so on.
>
> Since the Matrix will
Hello,
I do not know any such "community of R users in Spain and Latin America"
but it sounds like a great idea.
A number of R official documents have already been translated into
Spanish, but enhancing local language support on basic documentation
would facilitate adoption of the language by un
Hi
r-help-boun...@r-project.org napsal dne 20.02.2009 02:22:35:
> I have to admit I don't fully understand how MatLab's step() function
> works, but I'm learning that part. Can someone point me to
> equivalent&related functions in R?
Maybe ?step. Howewer without knowing what Matlab's step fu
Hi
r-help-boun...@r-project.org napsal dne 20.02.2009 12:54:41:
> Thank you for your advice, but I didn't manage to make it work...
> I tried
>
> carichi.annui <- data.frame
> (anno,loadPTG,loadPO4,loadNT,loadNH4,loadNO3,loadBOD5,loadSiO2)
>
> And I got this error message:
>
> Error in data.f
Hi
r-help-boun...@r-project.org napsal dne 20.02.2009 11:21:26:
>
> Hi
>
> I have a data set with two variables "q" and "depth" as follows:
>
> q<-c(tapply(weight[Soil=="Jy"], Depth[Soil=="Jy"], mean)). This commns
> returns 7 "q" values:
> 0.68790 0.84555 0.405416667 0.15277 0.033
On Fri, 20 Feb 2009, David Winsemius wrote:
If you can tolerate the subdiagonal and superdiagonal zero elements being
populated, then perhaps this is useful. If there is a subset function as you
suggest, then perhaps further reduction would be feasible from this as a
starting point. At least i
Here is one way of doing it. It does not check for errors in the case
that there are not enough values in the first dataframe:
> x1 <- read.table(textConnection("CHR_NR diffdatoperiode
+ 1137729
+ 1137759
+ 1137778"), header=TRUE)
>
> x2 <- rea
May you put your code and data in the internet?
It will be easy for others to find your real problem.
- Original Message -
From:
To:
Sent: Friday, February 20, 2009 7:54 PM
Subject: Re: [R] write.table
Thank you for your advice, but I didn't manage to make it work...
I tried
carichi.
On 2/20/2009 7:07 AM, Cetinyürek Aysun wrote:
Dear all,
I have a problem in calling my C code from R. If I do not need to use
#include, the codes work properly. But when I need to use it, my
codes produce an error as:
logpostCLM.cpp:7:15: error: R.h: No such file or directory
logpostCLM.cpp:8:24
Hi
r-help-boun...@r-project.org napsal dne 20.02.2009 10:34:02:
>
>
> Hi,
>
> I want to add a dashed vertical line to a number of xyplots.
>
> Here is a simple script of the type of plot I have but then I want to
add a
> reference line to 1995 on each of the panels. I have tried panel.ablin
If you can tolerate the subdiagonal and superdiagonal zero elements
being populated, then perhaps this is useful. If there is a subset
function as you suggest, then perhaps further reduction would be
feasible from this as a starting point. At least it would reduce the
size from 10^5 x 10^5
C.H. wrote:
Dear experts,
I would like to know how to plot the log-minus-log plot for survival
analysis (to check the proportional assumption) in R.
Using the AML example.
fit <- survfit(Surv(time, status) ~ x, data=aml)
length(fit$surv) #20
as the length of fit$surv is shorter than aml$x and
Have just added an example 12 on the home page:
http://sqldf.googlecode.com
that shows an example. Note use of notation
main.mytable to refer to an existing table in the
main database (as opposed to a data frame in R).
On Thu, Feb 19, 2009 at 11:55 PM, Stephen Tucker wrote:
> Hi all,
>
> I am
At 10:18 20/02/2009, mau...@alice.it wrote:
Your email has gone to hundreds of people all over the world many of
whom read their mail over quite slow connections especially in less
developed countries. Sending all of them half a megabyte of unwanted
mail is not very polite.
[mail content del
Hello,
I have a poblem with the Rpad server version. I've created a
three-step-calculation with Rpad. It consists of three .Rpad pages
linked to each other. It works perfect on he local version, because on
every page Rpad is communicating with the same R process. So all the
variables are sti
Dear all,
I tried help.search("karber") and RSiteSearch("karber") and
RSiteSearch("*karber*") to find whether the (trimmed) Spearman-Karber method
for LD50 evaluation in toxicity data (e.g. according to Hamilton 1977) has been
implemented in R. Or does this method feature under a different n
Dear all, I would like to hear about the best way between an dynamic web page
an R. (rpad, rapache, rphp son on ???).
I am producing dynamic web pages with PHP and need some data analysis in these
pages.
Best regards.
JL
_
[[e
Dear all,
I have a problem in calling my C code from R. If I do not need to use
#include, the codes work properly. But when I need to use it, my
codes produce an error as:
logpostCLM.cpp:7:15: error: R.h: No such file or directory
logpostCLM.cpp:8:24: error: Rinternals.h: No such file or directory
Note that using solve can be numerically unstable for certain problems.
On Fri, Feb 20, 2009 at 6:50 AM, Kenn Konstabel wrote:
> Decyphering formulas seems to be the most time consuming part of lm:
>
> mylm1 <- function(formula, data) {
> # not perfect but works
> F <- model.frame(formula
Thank you for your advice, but I didn't manage to make it work...
I tried
carichi.annui <-
data.frame(anno,loadPTG,loadPO4,loadNT,loadNH4,loadNO3,loadBOD5,loadSiO2)
And I got this error message:
Error in data.frame(anno, loadPTG, loadPO4, loadNT, loadNH4, loadNO3, :
the arguments have a di
Decyphering formulas seems to be the most time consuming part of lm:
mylm1 <- function(formula, data) {
# not perfect but works
F <- model.frame(formula,data)
y <- model.response(F)
mt <- attr(F, "terms")
x <- model.matrix(mt,F)
coefs <- solve(crossprod(x), crossprod(x,y))
Hello
I have one DF (detheleny1periode), with some variables that mathes, in some
way, variables in another DF (y2).
The DF named detheleny1periode look like this (i have not included alle
variables):
CHR_NR diffdatoperiode
1137729
1137759
11377
On Fri, 20 Feb 2009, lauramorg...@bluewin.ch wrote:
Sorry, I'm using R 2.8.1 on Microsoft Windows XP professional 2002 Service Pack
2.
The error I get is
Error in write.table(x, file, nrow(x), p, rnames, sep, eol,
na, dec, as.integer(quote), :
'list' type not implemented in 'EncodeElement
Using the zoo package:
Lines <- "SST_dateSST
2008-01-01 22.2
2008-01-02 21.8
2008-01-03 22.8
2008-01-04 22.9
2008-01-05 23.1
2008-01-06 23.2"
library(zoo)
# z <- read.zoo("myfile.dat", header = TRUE)
z <- read.zoo(textConnection(Lines), header = TRUE)
merge(z, avg3 = rollmean(z, 3, align = "
Dear R users,
I have trouble obtaining the same results for nested Anova with two fixed
factors when using lm and aov functions.
The formulas are:
> e1=aov(y~x/z)
> e2=lm(y~x/z)
summary(e1)
Df Sum Sq Mean Sq F valuePr(>F)
x 47 260.0 5.5 18.0088 < 2.
1 - 100 of 128 matches
Mail list logo
