Linda Lei bccrc.ca> writes:
>
> When dealing with the categorical covariates (for example 3 groups), it
> will come out different results if we add the command "factor" in front
> of the categorical covariate or not
The catch is here: the covariate is not a "categorial", but a number. In that
Hi
I have just started to have a look at R. I have used most stats software
packages and can use perl, visual basic etc. I am interested in how well
it handles lots of output e.g. tables or charts. How would you get lots
of output most easily and quickly into a Word document?
Sharon Snowdon
--
Thanks to Richard, Gabor and Marc for some nice solutions to my request.
I have a new problem:
> xify(30.10)
[1] "X.X"
> xify(30.11)
[1] "XXX.XXX"
The problem originates from:
> as.numeric(unlist(strsplit(as.character(15.10), "\\."))
Hi,
an other way without any libraries and written as a one-line-command
would be something like this:
xify <- function(x)
{
gsub("[0-9]","X",
sprintf(
paste(
"%",
ifelse(
format(x,
On Tue, 18 Jul 2006, Bashir Saghir (Aztek Global) wrote:
> Thanks to Richard, Gabor and Marc for some nice solutions to my request.
>
> I have a new problem:
>
> > xify(30.10)
> [1] "X.X"
> > xify(30.11)
> [1] "XXX.XXX"
>
> The problem
Am Montag, den 17.07.2006, 19:42 -0400 schrieb jim holtman:
> You were down at least 5300 levels in subroutine calls. Here is the
> first couple of lines from 'traceback()':
>
> 5363: vavb(a, b, n, v, x)
...
> 5316: vavb(a, b, n, v, x)
> The other thing is that there is no 'x' defined in the
Does it have to be a stop char ".", or could it be a separate
parameter, i.e. put in a comma, then the 10 becomes an integer...
s/
On 18/07/06, Prof Brian Ripley <[EMAIL PROTECTED]> wrote:
> On Tue, 18 Jul 2006, Bashir Saghir (Aztek Global) wrote:
>
> > Thanks to Richard, Gabor and Marc for some
Hi all,
Is there a simple way to convert an object name to a characters string?
Stéphane.
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.ht
> deparse(substitute(a))
[1] "a"
---
Jacques VESLOT
CNRS UMR 8090
I.B.L (2ème étage)
1 rue du Professeur Calmette
B.P. 245
59019 Lille Cedex
Tel : 33 (0)3.20.87.10.44
Fax : 33 (0)3.20.87.10.31
http://www-good.ibl.fr
--
In my opinion there is nothing wrong with asking the user to enter strings
as strings. I am working on some functions for SAS users (with limited R/S
experience) who are used to working with the SAS macro language. I was
trying to keep things simple for them by not forcing the use of quotes (a
styl
As I am trying to maintain consistency with some SAS macros I prefer not to
change the "standard definitions" in use at the moment. If I change things
too much I don't think many people with use my new functions.
Thanks,
Saghir
-Original Message-
From: [EMAIL PROTECTED] [mailto:[EMAIL PRO
On Tue, 18 Jul 2006, Stéphane Cruveiller wrote:
> Hi all,
>
> Is there a simple way to convert an object name to a characters string?
Yes, as.character, as in
> x <- as.name("foo")
> x
foo
> as.character(x)
[1] "foo"
However, I suspect you are not using the words in their technical sense
(a na
Hello,
I just updated from R 2.1.1 to R 2.3.1 (I also updated from Fedora 3 to
Fedora 4 if that is of any importance). However, several packages (e.g.
fields etc.) refuse to be installed. I get the following error message:
[EMAIL PROTECTED] R-files]# R CMD INSTALL fields_2.3.tar.gz
* Installing *
Il giorno lun, 17/07/2006 alle 15.00 -0400, Deepankar Basu ha scritto:
> I have been trying to read up the posting on the R-archive on this
> topic; but I could not really understand all the discussion, nor could I
> reach the "end". So, I am not aware of the current state of consensus on
> the iss
Thanks for your answer.
Here is what I am trying to do.
I have a list which is called MyList. I would like to get only the "name"
of this object as a simple characters string. i.e. Is there a function
in R which is able
to give:
> name<-fun(Mylist)
> name
"MyList"
thanks in advance,
Stéphan
On Tue, 18 Jul 2006, Stéphane Cruveiller wrote:
> Thanks for your answer.
>
> Here is what I am trying to do.
> I have a list which is called MyList. I would like to get only the "name"
> of this object as a simple characters string. i.e. Is there a function in R
> which is able
> to give:
>
> >
That is exactly what I wanted to do...
Thanks for the hint...
Stéphane.
Prof Brian Ripley a écrit :
> On Tue, 18 Jul 2006, Stéphane Cruveiller wrote:
>
>
>> Thanks for your answer.
>>
>> Here is what I am trying to do.
>> I have a list which is called MyList. I would like to get only the "nam
Hello all,
I am doing some spectral data analysis using R, and I already got some
help on how to do savitsky-golay preprocessing. Now I would like to
perform parameter selection (i.e. I would like to assess which spectral
lines are relevant to keep in the model, and which are not).
There exist
Hi,
only if you allow to input x as a string then you can use, maybe,
simply the following one-line-command:
xify <- function(x)
{
gsub(
"[0-9]",
"X",
sprintf(
paste("%",ifelse(regexpr("\\.",x) > 0, x, x<-paste(x,".
On Tue, 18 Jul 2006, Demmler J. wrote:
> Hello,
>
> I just updated from R 2.1.1 to R 2.3.1 (I also updated from Fedora 3 to
How? Did you install an RPM? Before or after OS update?
I checked the RPM for fc4 on CRAN which has
C compiler:gcc -O2 -g -pipe -Wp,-D_FORTIFY_SOURCE=
Hi,
you might try to use the R2HTML package and then import the html files into
word. please see. http://cran.r-project.org/doc/Rnews/Rnews_2003-3.pdf
Another way: use sweave for tex/latex output and transform it to rtf.
Friedrich Schuster
---
Your post was:
sharon snowdon sharonsnowdon at
> I was checking out your function and thinking of trying to make it more
> efficient when I noticed your comment above.
:-D
> Yes, bioinformatics is full
> of this stuff - have a look at Bioconductor.
I'm about it, but it'll take some time. It's a vast repository...
For the meantime, I've a
Hello
Using R v. 2.3.1 which I installed recently, when I start it as root, there
is nothing
I can do. I get "Error: function xx not found" even when trying to leave R
using
q(). R works fine when started by an ordinary user. I rechecked the root
environment
several times and even have set it simil
[EMAIL PROTECTED] writes:
> Hello
> Using R v. 2.3.1 which I installed recently, when I start it as root, there
> is nothing
> I can do. I get "Error: function xx not found" even when trying to leave R
> using
> q(). R works fine when started by an ordinary user. I rechecked the root
> environment
Hi!
Dear All,
I want to implement an autocorrelation analysis and estimated the Moran's I in
a set of
ecological traits. I seek for an R package that do this analysis, however, I
couldn't
found none that implement it. May be I'm lost in the universe of the Contributed
Packages but I couldn't fou
Problem solved: I had a .Renviron file in root's Home directory pointing to
the
wrong R installation...
Thanks anyway
Rick
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www
If its text that you want included in a larger Word document
you can create a text or HTML file from R and then when inserting it
into Word insert it as a link rather than copying it in.
output a file from R
In Word
Insert | File
browse to the file so its name appears in the File name box
Hi,
I am analysing survival data (diagnosis time until death/cens) with
time-dependent
covariates. I would like to fit a cox model using the (start, stop] variable.
In summary, I have the multiple internal time dependent covariates as follows;
1). LAS score (measured weekly on low r
Hi,
One week ago Sam Walker asked for a way to format a number to the
engineering notation?
After my suggestion I received many mails how to use this function
formatEng within plots dealing, e.g., with physics and SI prefixes.
So, I wrote a humble wiki article about that.
If you are interes
R can handle as many graphs/table as you ask it to.
I don't know what you're trying to do, but I generate lots of graphs using
loops and generating filenames on the spot.
See FAQ 7.34 to see how it's done.
http://cran.r-project.org/doc/FAQ/R-FAQ.html#How-can-I-save-the-result-of-each-iteration-i
Please look at
http://www.turbulence.org/Works/song/
This is a website by Martin Wattenberg that visually displays
the types of patterns you are looking for. He gave a paper
at the Joint Statistics Meetings in Minneapolis 2005.
__
R-help@stat.math.eth
On Tue, 2006-07-18 at 07:10 -0500, Ivan Rubio Perez wrote:
> Hi!
> Dear All,
>
> I want to implement an autocorrelation analysis and estimated the Moran's I
> in a set of
> ecological traits. I seek for an R package that do this analysis, however, I
> couldn't
> found none that implement it. May
Hi,
I have two further comments/questions about large datasets in R.
1. Does R's ability to handle large datasets depend on the operating
system's use of virtual memory? In theory, at least, VM should make the
difference between installed RAM and virtual memory on a hard drive
primarily a determ
Hi,
In a recent reply I sent to the list, my initial response was to
suggest, politely that by reading the posting guide and learning to use
the supplied tools one can often help oneself. I then went on to say
that reading the posting guide would have led you to RSiteSearch("insert
query"). Luckil
I think the posting guide may have been written prior to the
existence of RSiteSearch as a builtin command.
On 7/18/06, Gavin Simpson <[EMAIL PROTECTED]> wrote:
> Hi,
>
> In a recent reply I sent to the list, my initial response was to
> suggest, politely that by reading the posting guide and lear
In my experience, the OS's use of virtual memory is only relevant in the rough
sense that the OS can store *other* running applications in virtual memory so
that R can use as much of the physical memory as possible. Once R itself
overflows into virtual memory it quickly becomes unusable.
I'm n
Hello Vincent,
I think there is not much else than tcltk2. It is a complete
implementation of DDE (client/server), so it should fit your needs.
There are lots of examples on the man page. Please, if you think it
could be helpful for other users, do submit examples. I think real-time
acquisitio
On Tue, 18 Jul 2006, Marshall Feldman wrote:
> Hi,
>
> I have two further comments/questions about large datasets in R.
>
> 1. Does R's ability to handle large datasets depend on the operating
> system's use of virtual memory? In theory, at least, VM should make the
> difference between install
Hello. I am having trouble understanding the use of split.screen. I want
to divide the device surface first into 4 equal screens:
split.screen(figs=c(2,2))
This works.
I next want to subdivide each of these 4 screens into 10 subscreens. I do,
for the first of these 4 screens:
screen(1,new=
Hi R-users!
I am trying to create a what I call inflated array (maybe there is
already some other name for that). It is an array that changes
dinamically its dimensions, e.g. the higher the number of third
dimensions, the more rows in the array. So for example the array could
look like the followi
You can access DDE via COM as in this example
which uses DDE to open an Excel file. Note that
Excel also supports COM directly and normally
one would use COM with Excel, not DDE, so you
might check if your application also supports COM.
# opens an excel spreadsheet c:\test.xls using dde
library(R
On Tue, 18 Jul 2006, Bill Shipley wrote:
> Hello. I am having trouble understanding the use of split.screen. I want
> to divide the device surface first into 4 equal screens:
> split.screen(figs=c(2,2))
>
> This works.
>
> I next want to subdivide each of these 4 screens into 10 subscreens.
This means that the margins for 10 screens would take up more room than you
have - essentially the plot area is being squeezed to nothing. You can try
reducing your margins using par.
Also, it looks like you're trying to split into 20 screens there.
Hope this helps,
Matt
-Original Message
Hi,
I have a related question. How differently do other statistical
softwares handle large data?
The original post claims that 350 MB is fine on Stata. Some one
suggested S-Plus. I have heard people say that SAS can handle large
data sets. Why can others do it and R seem to have a problem? Don't
On Mon, 17 Jul 2006, John Wiedenhoeft wrote:
> Hi there,
>
> I'm having myself a hard time writing an algorithm for finding patterns
> within a given melody. In a vector I'd like to find ALL sequences that
> occur at least twice, without having to check all possible patterns via
> pattern matching
S-Plus stores objects as files whereas R stores them in memory.
SAS was developed many years ago when optimizing computer
resources was more important than it is now.
On 7/18/06, Ritwik Sinha <[EMAIL PROTECTED]> wrote:
> Hi,
>
> I have a related question. How differently do other statistical
> sof
> "Gabor" == Gabor Grothendieck <[EMAIL PROTECTED]>
> on Tue, 18 Jul 2006 09:29:07 -0400 writes:
Gabor> I think the posting guide may have been written prior
Gabor> to the existence of RSiteSearch as a builtin command.
yes, that's correct.
Thank you, Gavin, for the suggestion
Here is an example of the data.frame that I have,
df<-data.frame("case"=rep(1:5,each=9),"id"=rep(1:9,times=5),"x"=round(runif(length(rep(1:5,each=9)
"case" represents the cases,
"id" the persons, and
"x" is the binary state.
I would like to know in how many cases any two persons
a. both hav
On Tue, 18 Jul 2006, Ritwik Sinha wrote:
> Hi,
>
> I have a related question. How differently do other statistical
> softwares handle large data?
>
> The original post claims that 350 MB is fine on Stata. Some one
> suggested S-Plus. I have heard people say that SAS can handle large
> data sets. W
Am Dienstag, den 18.07.2006, 22:09 -0400 schrieb Jim Lemon:
> Hi John,
>
> Minor bug - I zeroed the hit counter in the wrong place.
>
> find.replay<-function(tunestring,maxlen) {
> return(matchlist)
> }
Dear Jim,
many, many thanks for your effords :-D!!! Your program is great and very
el
library(gtools)
apply(combinations(9,2), 1, function(x) with(df[df$id %in% x, ], table(x, id)))
---
Jacques VESLOT
CNRS UMR 8090
I.B.L (2ème étage)
1 rue du Professeur Calmette
B.P. 245
59019 Lille Cedex
Tel : 33 (0)3.20.87.10.44
Fax
Hi,
I am trying to fit a model with an ordered response variable (3 levels) and
13 predictor variables. The sample has complex survey design and I've used
'svydesign' command from the survey package to specify the sampling design.
After reading the manual of 'svyglm' command, I've found that you c
Hello!
I work with:
R : Copyright 2006, The R Foundation for
Statistical Computing
Version 2.3.1 (2006-06-01)
On Windows XP Professional (Version 2002) SP2
I think there is a bug in the conditional
execution if (expr1) {expr2} else {expr3}
If I try:
"if (expr1) expr2 else expr3"
it works w
Hello!
I work with :
R : Copyright 2006, The R Foundation for
Statistical Computing
Version 2.3.1 (2006-06-01)
On Windows XP Professional (Version 2002) SP2.
At this moment I use the function "nls" combined
with a selfStar model (SSmicmen, related to
Michaelis-Menten equation, and provided by the
Hadassa Brunschwig gmail.com> writes:
>
> Hi R-users!
>
> I am trying to create a what I call inflated array (maybe there is
> already some other name for that). It is an array that changes
> dinamically its dimensions, e.g. the higher the number of third
> dimensions, the more rows in the arra
Thanks Gabor and Philippe.
Special thanks to Philippe for his tcltk2 nice job.
I'm testing differents approaches for my problem.
I'll return info if I use DDE + tcltk2.
Vincent
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listin
Hi
I would like to know if there is a package that allows me to implement
endogenous Tobit/Probit model.
Example:
res1 <- rnorm(N);
res2 <- res1*0.5 + rnorm(N)
x <- z[,1]*2 + res1;
ys <- x*b + res2;
d <- (ys>0); #dummy variable
y <- d*ys;
y is censored and x is correlated
On Tue, 18 Jul 2006, Xavier Barron wrote:
Hello!
I work with:
R : Copyright 2006, The R Foundation for
Statistical Computing
Version 2.3.1 (2006-06-01)
On Windows XP Professional (Version 2002) SP2
I think there is a bug in the conditional
execution if (expr1) {expr2} else {expr3}
If I try:
Hi
your fit is an object (list) and you could use some functions like
summary or coef to extract usefull information from it or you can
call its components on your own.
> DNase1 <- subset(DNase, Run == 1)
>
> ## using a selfStart model
> fm1DNase1 <- nls(density ~ SSlogis(log(conc), Asym, xmi
Hi
not at all
works for me
> a<-1
> b<-2
> if (a>b) {print("Hallo")} else {print("OK")}
[1] "OK"
> if (a
you probably started your else on new line e.g.
> if (a if (a
HTH
Petr
On 18 Jul 2006 at 17:43, Xavier Barron wrote:
Date sent: Tue, 18 Jul 2006 17:43:42 +0200 (CEST)
Fr
Others have suggested using R2HTML (which is a good option).
Another option is to use Sweave and specifically the new odfWeave
package for R. This works on OpenOffice files rather than word files
(but OpenOffice http://www.openoffice.org/ can inport and export word
documents).
The basic idea is
Dear R-Listeners,
as the Sweave faq says:
http://www.ci.tuwien.ac.at/~leisch/Sweave/FAQ.html
creating several figures from one figure chunk does not work, and for
standard graphics, a workaround is given. Now I have a multipage trellis
plot with an a-priori unknown number of pages, and I don't s
On Tue, 18 Jul 2006, Debarchana Ghosh wrote:
> Hi,
>
> I am trying to fit a model with an ordered response variable (3 levels) and
> 13 predictor variables. The sample has complex survey design and I've used
> 'svydesign' command from the survey package to specify the sampling design.
> After read
Greetings, fellow R'ers.
How can I get this frame in R:
ID meas ID.1 meas.1
1 1.13 1.2
2 2.14 2.2
to look like this (stacking):
ID meas
1 1.1
2 2.1
3 1.2
4 2.2
It's not really the reshape function (or is it?) because we can consider
the additional colum
Thanks a lot for all the responses. The general drift of all the messages was
the suggestion to use some database management package that has a nice
interface with R; and most of the suggestions pointed in the direction of SQL.
I will look into the SQL package and start learning to use it along
Try this:
# set up test data
Lines <- "ID meas ID.1 meas.1
1 1.13 1.2
2 2.14 2.2
"
DF <- read.table(textConnection(Lines), header = TRUE)
# reshape
matrix(t(DF), nc = 2, byrow = TRUE, dimnames = list(NULL, colnames(DF)[1:2]))
On 7/18/06, Jesse Albert Canchola <
Hi,
What is the best general strategy to prevent returned errors from
interrupting whatever it is I am running? Only by using options()?
This is a problem for me in 2 particular cases:
(i) An error breaking my loops.
(ii) I would like to run some regressi
Sorry, in looking at this again my previous code did not give the
same ordering you indicated. Instead using the same DF try
this:
rbind(as.matrix(DF[,1:2]), as.matrix(DF[,3:4]))
Both this and the last piece of code produce matrices so
use as.data.frame if you want a data frame.
On 7/18/06, Gab
I am thrilled to learn tcltk2 has DDE capability.
It is the piece I have been needing to make ESS work directly
with the RGUI on Windows. GNU emacs on Windows has a ddeclient,
but no access to COM. So if R, or tcltk2 talking in both directions to R,
has a ddeserver, all should be possible. I wil
I think this is an easy question, but I would be grateful for any advice
on how to implement this in R.
I simply have a response variable (y) that I am trying to predict with
one explanatory variable (x) but the shape of the scatter plot is
distinctly bilinear. It would be best described by two st
Many thanks, Gabor. That worked great! I'm ecstatic.
Best regards,
Jesse
"Gabor Grothendieck" <[EMAIL PROTECTED]>
07/18/2006 10:12 AM
To
"Jesse Albert Canchola" <[EMAIL PROTECTED]>
cc
r-help@stat.math.ethz.ch
Subject
Re: [R] Reconfiguring wide frame to long frame
Try this:
# set up
Thanks, Gabor. Since the data stacking components are independent, that
didn't matter much but I am grateful for your follow-up code to match the
desired output specifically.
Regards,
Jesse
"Gabor Grothendieck" <[EMAIL PROTECTED]>
07/18/2006 10:36 AM
To
"Jesse Albert Canchola" <[EMAIL P
It appears that you might have a latent (hidden) explanatory variable that
causes the two-population appearance. If you have some ideas on what that
other factor might be, you could try two separate linear regressions for
each value of the latent factor and compare the slopes and intercepts. You
Hello,
suppose I have a multivariate multiple regression model such as the
following:
> DF<-data.frame(x1=rep(c(0,1),each=50),x2=rep(c(0,1),50))
> tmp<-rnorm(100)
> DF$y1<-tmp+DF$x1*.5+DF$x2*.3+rnorm(100,0,.5)
> DF$y2<-tmp+DF$x1*.5+DF$x2*.7+rnorm(100,0,.5)
> x.mlm<-lm(cbind(y1,y2)~x1+x2,dat
Dear R-helpers,
In JGR, how to I get the help() to update when I install a new package?
> sessionInfo()
Version 2.3.1 (2006-06-01)
powerpc-apple-darwin8.6.0
attached base packages:
[1] "datasets" "methods" "stats" "graphics" "grDevices"
"utils" "base"
other attached packages:
I am having trouble fitting correlation structures within nlme. I would like to
fit corCAR1, corGaus and corExp correlation structures to my data. I either
get the error "step halving reduced below minimum in pnls step" or
alternatively R crashes.
My dataset is similar to the CO2 example in th
Hi,
I find a lot of the R-help email traffic overloads my inbox. My IT
managers are not really happy for me to be subscribed to several
high-traffic email lists. I don't want to lose my contact with the
R-help emails, so I'm having to consider various ways of handling the
traffic. Anyhow, I'm w
All,
Recently I ran across a URL documenting published research using R:
http://www.cgd.ucar.edu/ccr/ammann/millennium/CODES_MBH.html
A note on the site indicates that the code is being revised. The code and
data are provided, so that one could reproduce the results without having to
buy a prop
On Tue, 2006-07-18 at 11:30 -0700, Darren Weber wrote:
> Hi,
>
> I find a lot of the R-help email traffic overloads my inbox. My IT
> managers are not really happy for me to be subscribed to several
> high-traffic email lists. I don't want to lose my contact with the
> R-help emails, so I'm havi
On 7/18/2006 2:30 PM, Darren Weber wrote:
> Hi,
>
> I find a lot of the R-help email traffic overloads my inbox. My IT
> managers are not really happy for me to be subscribed to several
> high-traffic email lists. I don't want to lose my contact with the
> R-help emails, so I'm having to conside
Hi,
When I used
x <- c(10,10,10,10,5,5,5,5)
sum(x)
x <- x - .Last.value/8
the result was zero, as I expected.
However, when I used
x <- rnorm(101,0,10)
sum(x)
x <- x - .Last.value/101
sum(x)
I did not get zero, but -2.664535e-15. OK, that's fairly close to
zero, but I tried to improve the figur
Dear Ulrich,
I'll look into generalizing linear.hypothesis() so that it handles
multivariate linear models.
Meanwhile, vcov(x.mlm) will give you the covariance matrix of the
coefficients, so you could construct your own test by ravelling
coef(x.mlm) into a vector.
I hope that this helps,
John
On Tue, 18 Jul 2006, Jim Frederick wrote:
> Hi,
>
> When I used
> x <- c(10,10,10,10,5,5,5,5)
> sum(x)
> x <- x - .Last.value/8
> the result was zero, as I expected.
> However, when I used
> x <- rnorm(101,0,10)
> sum(x)
> x <- x - .Last.value/101
> sum(x)
> I did not get zero, but -2.664535e-15.
Dear R-helpers,
Can someone point me to an extended example of use of R2HTML?
_
Professor Michael Kubovy
University of Virginia
Department of Psychology
USPS: P.O.Box 400400Charlottesville, VA 22904-4400
Parcels:Room 102Gilmer Hall
McCormick
Well, SPSS used to claim that all its algorithms dealt with only one case at
a time and therefore that it could handle very large files. I suppose a
large correlation matrix could cause it problems.
Marsh Feldmman
-Original Message-
From: Ritwik Sinha [mailto:[EMAIL PROTECTED]
Se
Dear R-ers,
I conducted a regression analysis, and then intended to add the regression
function (y=4.33+1.07x) in a graph. But the following code can only give me a
text like y=a+bx. Who can help me out? Thank you very much in advance.
CODE
# Read data
x <- c(1, 2, 3, 4, 5
On Tue, 18 Jul 2006, junguo liu wrote:
> Dear R-ers,
>
> I conducted a regression analysis, and then intended to add the
> regression function (y=4.33+1.07x) in a graph. But the following code
> can only give me a text like y=a+bx. Who can help me out? Thank you very
> much in advance.
This i
[Thomas Lumley]
>People have used R in this way, storing data in a database and reading it
>as required. There are also some efforts to provide facilities to support
>this sort of programming (such as the current project funded by Google
>Summer of Code: http://tolstoy.newcastle.edu.au/R/devel
Hello!
Say I have a function, which creates a design matrix i.e.
myFunc <- function(x)
{
ret <- cbind(x, x*x, x*x*x)
colnames(ret) <- 1:ncol(ret)
return(ret)
}
n <- 200
x <- runif(n=n, min=0, max=100)
y <- myFunc(x) %*% c(1, 0.2, -0.0002) + rnorm(n=n, sd=100)
then I can use this in formu
Or, more succinctly, "Pinard's Law":
The demands of ever more data always exceed the capabilities of ever better
hardware.
;-D
-- Bert Gunter
Genentech Non-Clinical Statistics
South San Francisco, CA
> -Original Message-
> From: [EMAIL PROTECTED]
> [mailto:[EMAIL PROTECTED] On Behal
You could plot y vs. fitted(y.lm) where y.lm is the output of lm or
plot both y and fitted(y.lm) against x on the same chart.
On 7/18/06, Gregor Gorjanc <[EMAIL PROTECTED]> wrote:
> Hello!
>
> Say I have a function, which creates a design matrix i.e.
>
> myFunc <- function(x)
> {
> ret <- cbind(
Can someone please help im trying to use JGR' version 1.4-2 on a IBM
Thinkcentre P4 with windows xp with the R v2.3.1, but any time I try and
plot or use JavaGD I get this error
> plot(1:10)
Error in JavaGD() : unable to start device JavaGD
In addition: Warning message:
Another VM is runnin
Hi Uli,
I suggest that you try to rewrite the model system into a single
mixed-effects model, which would allow direct parameterization of the
tests that you're interested in. A useful article, which may be
overkill for your needs, is:
Hall, D.B. and Clutter, M. (2004). Multivariate multilevel n
I've been meaning send an announcement for this package, but Greg Snow
beat me to the punch today.
Max
The odfWeave package is now available on CRAN at
http://lib.stat.cmu.edu/R/CRAN/src/contrib/Descriptions/odfWeave.html
and your local mirror.
The package extends Sweave to Open Document For
*
* "methods" *
*
You have asked an excellent question. I can provide a partial answer
below. First, however, I wish to pose a question of my own, which could
help answer your question:
How can one obtain a simple list of the available generics for a
class?
Hello, looking to draw a voronoi tessellations in R - can anyone
recommend a package that has tackled this?
some background:
i have a economic data set and created a sammons projection, like to now
overlay a voronoi tessellation over the sammons 2-D solution for a slick
visual, and potentially
Hi,
I'm trying to fit a ordered logistic regression. The response variable
(y) has three levels (0,1,2).
The command I've used is:
/ordlog<-polr(y~x1+x2+x3+x4, data=finalbase, subset=heard, weight=wt,
na.action=na.omit)
/
(There are no NA's in y but there are NA's in X's)
The error I'm getting
Hi,
I'm trying to fit a ordered logistic regression. The response variable
(y) has three levels (0,1,2).
The command I've used is:
ordlog<-polr(y~x1+x2+x3+x4, data=finalbase, subset=heard, weight=wt,
na.action=na.omit)
(There are no NA's in y but there are NA's in X's)
The error I'm getting is
I'll suggest going to the CRAN packages page and doing a search for "voronoi".
Also, search for 'triangulation', since that is one of the uses of them.
-Don
At 11:46 PM -0400 7/18/06, zubin wrote:
>Hello, looking to draw a voronoi tessellations in R - can anyone
>recommend a package that has tack
Hi,
Can anyone pls help me in plotting the following data?
The data-set contains company name, specification-1, specification-2.
The graph would basically plot company name with specification-1
on x-axis, & specification-2 on y-axis.
Simply put - company name should get class
1 - 100 of 103 matches
Mail list logo