[EMAIL PROTECTED] wrote:
I am working on a Dell Latitude D800 laptop computer with an Intel
processor running Windows XP Professional. The problem I describe below
has occurred while running R versions 1.9.1 and 2.0.1. The error occurs
while using a program to perform categorical regression that
I'm sure others with more experience will answer this, but for what it is worth
my experience suggests that memory issues are more often with the user and not
the machine. I don't use Linux so I can't make specific comments about the
capacity of your machine. However it appears that there is oft
On Thu, 3 Feb 2005, Dan Lipsitt wrote:
I am trying to write a C function that reads lines from a file and
passes them back to R.
Please take a closer look at the posting guide, which clearly indicates
this was the wrong list:
Questions likely to prompt discussion unintelligible to non-programme
> Suppose I have a linear mixed-effects model (from the package nlme) with
> nested random effects (see below); how would I present the results from
the random effects part in a publication?
>
Have you tried anova(lme())?
Your asin(sqrt()) looks a bit like these are percentages of counts.
I'm trying to read in datasets with roughly 150,000 rows and 600
features. I wrote a function using scan() to read it in (I have a 4GB
linux machine) and it works like a charm. Unfortunately, converting the
scanned list into a datafame using as.data.frame() causes the memory
usage to explode (it c
Hi
maybe there are other options but
x[rowSums(x==-99)>=1,]
seems to extract all rows which have at least one -99 in them
Cheers
Petr
On 3 Feb 2005 at 14:57, Tim Howard wrote:
> I am trying to extract rows from a data.frame based on the
> presence/absence of a single value in any column. I'
Hi
I have put up some web pages containing a number of plots (and diagrams)
produced using R (they correspond to the figures for a book that I am
working on about R graphics), with the relevant R code provided for each
plot (or diagram), at
http://www.stat.auckland.ac.nz/~paul/RGraphics/rgraphics.
Please excuse the off topic question, but some of you
certainly have very significant computer savvy beyond
R and statistics and I would like to tap that
knowledge. I do not believe this to be an R problem,
but to date I believe that R is the only application
on my machine that has been impacted. T
On Thu, 3 Feb 2005, Christofer Ecklund wrote:
> Could you tell me how to uninstall R from Mac OS X (10.3)? I am not
> too familiar with terminal or the command line so a GUI removal would
> be better, however I would like to remove every trace of the program
> from my system. Usually with a pack
Statistics Practitioner Fluent in R - San Francisco CA
Loyalty Matrix Inc., in downtown San Francisco, is expanding our team. We
are a young, dynamic and growing team of multidisciplinary marketing and
technical professionals. We deliver value to our clients by discovering
actionable tactical and
Have you considered Monte Carlo?
Invent a test, then Monte Carlo it to get a p-value.
Something similar to what you are asking is provided by Monte
Carlo confidence bounds on a QQ plot, discussed on this list last June.
To find it, go to www.r-project.org -> Search -> "R s
On Mon, Jan 31, 2005 at 07:52:26PM +0100, silvia simoni wrote:
> I have a problem in installing R-2.0.1, downloaded from the R web site,
> on mac o sx version 10.3.7.
> when i launch the command ./configure i get the following error message:
Did you read the OS X FAQ on www.r-project.org? There
I tried to use shapiro.test or ks.test to check the
normality of some data, the problem is, the
distribution function is a mixture of a Gaussian and
some other distributions at the tails. The hypothesis
is that if the tails are excluded, the distribution is
perfect Gaussian, and I want to test that
I am trying to write a C function that reads lines from a file and
passes them back to R.
Here is a simplified version:
--- C code ---
#include
void read_file(char **filename, char **buf, char **buflen) {
FILE *infile;
infile = fopen(*filename, "r");
fgets(*buf, *buflen, infile);
fclos
I did not attempt to follow your code or discussion but you could
try these:
1. try to pin down what part of your code is taking the time
2. try to eliminate the loop, if possible
3. use matrices rather than data frames -- matrices are faster
Daniel E. Bunker columbia.edu> writes:
:
: Hi All
I am working on a Dell Latitude D800 laptop computer with an Intel
processor running Windows XP Professional. The problem I describe below
has occurred while running R versions 1.9.1 and 2.0.1. The error occurs
while using a program to perform categorical regression that was
originally written in
"Tim Howard" <[EMAIL PROTECTED]> writes:
> I am trying to extract rows from a data.frame based on the
> presence/absence of a single value in any column. I've figured out how
> to do get the positive matches, but the remainder (rows without this
> value) eludes me. Mining the help pages and arch
Do the -99 entries really mean NA? In that case, I think it
would be clearer to recode your data frame with NAs and then select
out the complete or incomplete rows:
x[x == -99] <- NA
x[compete.cases(x),] # or na.omit(x)
x[!complete.cases(x),]
Tim Howard gw.dec.state.ny.us> writes:
:
: ap
On Thu, 03 Feb 2005 14:57:58 -0500 Tim Howard wrote:
> I am trying to extract rows from a data.frame based on the
> presence/absence of a single value in any column. I've figured out
> how to do get the positive matches, but the remainder (rows without
> this value) eludes me. Mining the help pa
On Thu, 2005-02-03 at 14:57 -0500, Tim Howard wrote:
> I am trying to extract rows from a data.frame based on the
> presence/absence of a single value in any column. I've figured out how
> to do get the positive matches, but the remainder (rows without this
> value) eludes me. Mining the help pa
Hello,
Could you tell me how to uninstall R from Mac OS X (10.3)? I am not
too familiar with terminal or the command line so a GUI removal would
be better, however I would like to remove every trace of the program
from my system. Usually with a package installer there is some sort of
uninstal
I tried it on Windows XP with R 2.1.0 and could not replicate it either.
Suggest you start up a fresh session and try it again.
By the way, you could consider this:
xtabs(count ~., data.frame(sex = sex, income = income))
Carlos Ortega yahoo.es> writes:
:
: Hi,
:
: That is something strang
Check out ?reshape .
b3e.jussieu.fr> writes:
:
: So in fact :
:
: I have a data wich colnames are :
: "id","VEMS5A","VEMS6A","VEMS7A","VEMS8A","VEMS9A","VEMS10A","VEMS11A","VEMS12
A","VEMS13A","VEMS14A","VEMS15A","VEMS16A","VEMS17A","VEMS18A","cohort")
:
: With "VEMS5A" is the value when t
apply, of course, does the trick exceptionally well. Thank you,
everyone, for the help.
tim
>>> Chuck Cleland <[EMAIL PROTECTED]> 02/03/05 03:10PM >>>
How about this?
#extract data.frame of rows with -99 in them
subset(x, apply(x, 1, function(x){any(x == -99)}))
#extract data.frame of rows not
Vera Hofer asks:
" Dear colleagues,
has anyone an idea how to carry out a nonparametric manova for comparing
3 groups?
Thank you for your help.
Vera"
One approach would be to carry out a permutation test, randomly
reassigning the vector observations among groups under the null
hypothesis that t
On the first issue, there's a recent post you can find in the archives.
It's from Deepayan Sarkar on January 2 this year. It would probably pop up
on a search for "black white lattice" or something similar.
The key part of his answer:
I'd do something like this as part of the initialization:
<<
On Thursday 03 February 2005 13:52, Dean Sonneborn wrote:
> I'm working on a graphic but have run into a road block about two
> issues. I don't have a color printer so I want to produce the
> graphic in black and white. I'm currently using this
> statement trellis.device(bg="white")
Doesn't th
How about this?
#extract data.frame of rows with -99 in them
subset(x, apply(x, 1, function(x){any(x == -99)}))
#extract data.frame of rows not containing -99 in them
subset(x, apply(x, 1, function(x){all(x != -99)}))
hope this helps,
Chuck Cleland
Tim Howard wrote:
I am trying to extract rows from
Tim Howard wrote:
I am trying to extract rows from a data.frame based on the
presence/absence of a single value in any column. I've figured out how
to do get the positive matches, but the remainder (rows without this
value) eludes me. Mining the help pages and archives brought me,
frustratingly,
I am trying to extract rows from a data.frame based on the
presence/absence of a single value in any column. I've figured out how
to do get the positive matches, but the remainder (rows without this
value) eludes me. Mining the help pages and archives brought me,
frustratingly, very close, as yo
I'm working on a graphic but have run into a road block about two
issues. I don't have a color printer so I want to produce the graphic in
black and white. I'm currently using this
statement trellis.device(bg="white") but the body of the graphic contains
color. What is the code to create the
So in fact :
I have a data wich colnames are :
"id","VEMS5A","VEMS6A","VEMS7A","VEMS8A","VEMS9A","VEMS10A","VEMS11A","VEMS12A","VEMS13A","VEMS14A","VEMS15A","VEMS16A","VEMS17A","VEMS18A","cohort")
With "VEMS5A" is the value when the subject (id) was 5 year old,... "VEMS18A" is
the value when the
Hi,
That is something strange, I could not replicate it...
Regards,
Carlos.
+++
> version
_
platform i386-pc-mingw32
arch i386
os mingw32
system i386, mingw32
status
major2
In the future please be a little more specific with your questions. i.e.
give us an example of your data the way it looks now and the way you want
it to look.
I think what you want is something like:
time<-seq(as.Date("2005/1/1"), as.Date("2005/1/10"), "days")
age<-round(runif(10,30,50))
measur
XSolutions Corp (www.xlsolutions-corp.com) is proud to announce
our "Advanced R/Splus programming" course taught by R Development
Core Team Guru!
www.xlsolutions-corp.com/Radv.htm
*San Francisco --> March 31st-April 1st, 2005
*Boston, MA--> TBD
Early-bird d
Dear all,
I wanted to make a two-way-table of two variables with a counting
variable stored in another column of a dataframe. In version 1.9.1, the
behavior is as expected as shown in the simplified example code.
> sex <- rep(c("F", "M"), 5)
> income <- c(rep("low", 5), rep("high", 5))
> count <
Well, in principle there are general stats lists and newsgroups, and
this is not specifically an R question. However, people around here
have been known to contribute information about statistical substance
from time to time. You're still not excused for not using a meaningful
subject line though.
Excuse me.
I don't understant how to use the function "stack".
I'm not sure to know the good syntax.
I've got a data witch column names are the time when the mesures was done and
an
other who indicates the age of the subjetcs.
I would like to do a new matrix with patients older than a value and
Le 03.02.2005 16:42, Romain Francois a écrit :
Hello,
try the following :
l <- list(NULL)
l[[1]] <- 1:2
l[[2]] <- 2:3
l[[3]] <- 4:5
l
[[1]]
[1] 1 2
[[2]]
[1] 2 3
[[3]]
[1] 4 5
matrix(unlist(l),ncol=2)
Sorry, i forgot the byrow argument, the correct comand was :
matrix(unlist(l),ncol=2,byrow=T)
I n
Hi,
I am puzzled by the relationship between the p-values asociated with the
coefficients of a univariate logistic regression involving categorical
variables and the p-value I get from Fisher's exact test of the
associated 2 x 2 contingency table.
(1) The 2-sided p-value for the table is ~ 0.00
Hi All,
I am writing a simulation that examines the effects of species
extinctions on ecological communties by sequentially removing
individuals of a given species (sometimes using weighted probabilities)
and replacing the lost individuals with species identities randomly
sampled from the remai
Charles and Kimberly Maner wrote:
Hi all. I'm reading in a flat, comma-delimited flat file using read.csv.
It works marvelously for the most part. I am using the colClasses argument
to, basically, create numeric, factor and character classes for the columns
I'm reading in. However, a couple of t
Hello,
try the following :
l <- list(NULL)
l[[1]] <- 1:2
l[[2]] <- 2:3
l[[3]] <- 4:5
l
[[1]]
[1] 1 2
[[2]]
[1] 2 3
[[3]]
[1] 4 5
matrix(unlist(l),ncol=2)
[,1] [,2]
[1,]13
[2,]24
[3,]25
Is that what you want ?
Romain
Le 03.02.2005 16:08, michael watson (IAH-C) a écrit :
H
list1 <- list(x=c(1,2), y=c(3,4), z=c(4,5))
matrix(unlist(list1), nrow=length(list1), byrow=T)
HTH
On Thu, 3 Feb 2005, michael watson (IAH-C) wrote:
> Hi
>
> Sorry to ask such a basic question. I have a list, each element of
> which is a vector of two values. What I actually want is a matrix w
Will do.call('rbind',yourlist) do what you want?
Sean
On Feb 3, 2005, at 10:08 AM, michael watson ((IAH-C)) wrote:
Hi
Sorry to ask such a basic question. I have a list, each element of
which is a vector of two values. What I actually want is a matrix with
two columns, and one row per element of t
do.call() is your friend:
mat <- do.call("rbind", myList)
Andy
> From: michael watson (IAH-C)
>
> Hi
>
> Sorry to ask such a basic question. I have a list, each element of
> which is a vector of two values. What I actually want is a
> matrix with
> two columns, and one row per element of th
try this:
mat <- matrix(unlist(lis), ncol=2, byrow=TRUE)
I hope it helps.
Best,
Dimitris
Dimitris Rizopoulos
Ph.D. Student
Biostatistical Centre
School of Public Health
Catholic University of Leuven
Address: Kapucijnenvoer 35, Leuven, Belgium
Tel: +32/16/336899
Fax: +32/16/337015
Web: http://w
Faheem Mitha wrote:
On Sun, 30 Jan 2005, Paul Roebuck wrote:
and if you use the following instead of "Rinternals.h"?
extern "C" {
#include
}
Still the same errors. As Dirk pointed out, putting iostream first makes
the errors go away in either case.
You have to define R_NO_REMAP in your code bef
Hi
Sorry to ask such a basic question. I have a list, each element of
which is a vector of two values. What I actually want is a matrix with
two columns, and one row per element of the list. Obviously I have
tried as.matrix(), and as.vector() but I didn't expect the latter to
work.
I feel so l
Hello!
I want R to call some functions from my DLL, which is written on C++
using Microsoft Visual C++ 6.0. I've looked through "Writing R
extensions" manual and "FAQ for Windows". I achieved, that R can
call void functions without parameters from my C++ DLL.
Then I included R headers files to my
If I'm not mistaken, you are dealing with a total of:
> 414180 * 1515 * 4 / 1024^2
[1] 2393.657
or close to 2.4GB, for a single copy of the (combined) array. Unless you're
on a 64-bit OS running on hardware with sufficient physical RAM, that's not
likely to work so well.
Also, I believe you sho
Chandra H. wrote:
Hello,
I am Hukum Chandra and I'm a statistician. I am using R. I have gone
through R but I could not get the way how to find "Information matrix"
in R. In particular, var-Cov matrix of components of variances in
linear mixed model. Could you please help me in getting the way to
Hi all
I am trying to use R for some data editing using the "array" function to
write binary data to a text file. I realise that R was not designed for
this purpose but I am no C programmer and would prefer to use R (as I
know how to do it and hate c). Basically I get the following error.
Error
Hello,
I am Hukum Chandra and I'm a statistician. I am using R. I have gone
through R but I could not get the way how to find "Information matrix"
in R. In particular, var-Cov matrix of components of variances in
linear mixed model. Could you please help me in getting the way to
produce it using R
I am running R 2.0.1 and I get this error message when I'm trying to install
different packages: No package 'file10181' was found.
Can anyone help?
JJ
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read t
Hi all. I'm reading in a flat, comma-delimited flat file using read.csv.
It works marvelously for the most part. I am using the colClasses argument
to, basically, create numeric, factor and character classes for the columns
I'm reading in. However, a couple of the fields in the file are date
fi
michael watson (IAH-C) wrote:
Hi
I found this post from 2001
(https://stat.ethz.ch/pipermail/r-help/2001-April/011073.html) and was
just wondering if this has been updated and included in an R package
yet?
Cheers
Mick
see mantel() in package vegan on CRAN. The web site for vegan is:
http://cc.oulu.
library(Design)
d <- datadist(mydata); options(datadist='d')
f <- lrm(incidence ~ sun*trees) # lrm is for binary or ordinal response
plot(f, sun=NA, trees=NA)
# add method='image' or 'contour' to get other types of graphs
plot(f, sun=NA, trees=NA, fun='plogis') # probability scale
Correction: fun
Hi
I found this post from 2001
(https://stat.ethz.ch/pipermail/r-help/2001-April/011073.html) and was
just wondering if this has been updated and included in an R package
yet?
Cheers
Mick
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/ma
[EMAIL PROTECTED] wrote:
Dear R-helpers,
I tried to create a 3D surface showing the interaction between two
continuous explanatory variables; the response variable is binary (0/1).
The model is:
model<-glm(incidence~sun*trees,binomial)
then I used "wireframe" to create a 3D plot:
xyz<-expand.grid(s
Dear all,
I have just uploaded to CRAN the version 1.2 of the package
'adehabitat'. This package now contains new functions for
the estimation of home ranges (Nearest neighbour convex
hull, brownian bridge approach for kernel estimation),
estimation of schoener's ratio. This package now allows a
b
> I tried to create a 3D surface showing the interaction between two
> continuous explanatory variables; the response variable is
> binary (0/1).
>
> The model is:
>
> model<-glm(incidence~sun*trees,binomial)
>
> then I used "wireframe" to create a 3D plot:
>
> xyz<-expand.grid(sun=seq(30,180,
On Tuesday 01 February 2005 06:42 am, Tae-Hoon Chung wrote:
> Hi, All;
>
> When I generate a "special" paper postscript image larger than "a4" or
> "letter" using R, I can only see one-page portion of all image, of course.
> What will be the simple solution for this? Is there any way I can set the
Dear R-helpers,
I tried to create a 3D surface showing the interaction between two
continuous explanatory variables; the response variable is binary (0/1).
The model is:
model<-glm(incidence~sun*trees,binomial)
then I used "wireframe" to create a 3D plot:
xyz<-expand.grid(sun=seq(30,180,1),tre
Hi Mike,
Do you have a schematic drawing of how exactly your treatments were
applied? In split-plot experiments, it is generally very important to
clearly define the sequence of plot sizes, because if you don´t do this
properly, then the output will be confusing. Checking if your degrees of
fre
Hi Vera,
See:
http://finzi.psych.upenn.edu/R/Rhelp02a/archive/22273.html
it seems like your problem.
Best regards,
Vito
you wrote:
Dear colleagues,
has anyone an idea how to carry out a nonparametric
manova for comparing
3 groups?
Thank you for your help.
Vera
=
Diventare costruttori d
Hello,
Following on from a paper by Ludwig, Horel and Whiteman in 2004, I have
been using PCA to identify surface wind patterns. In the paper they
approach the problem using the "real vector approach", that is, they
create paired (u,v) vectors for the wind data, in effect artificially
doubling the
Dear colleagues,
has anyone an idea how to carry out a nonparametric manova for comparing
3 groups?
Thank you for your help.
Vera
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! htt
Dear useRs,
version 1.1-0 of the fortunes package is available from CRAN. It now
contains 100 R fortunes. Thanks to all who contributed.
Best wishes,
Z
___
R-packages mailing list
[EMAIL PROTECTED]
https://stat.ethz.ch/mailman/listinfo/r-packages
69 matches
Mail list logo