Frank E Harrell Jr wrote:
Cody,
How amazing that SAS is still used to produce reports that reviewers
hate and that requires tedious low-level programming. R + LaTeX has it
all over that approach IMHO. We have used that combination very
successfully for several data and safety monitoring rep
Hi:
Another solution is found in package BHH2, associated with the second
edition of the classic Box, Hunter and Hunter text. Illustration of both
BHH2 and plotrix:
x <- c(2, 4, 5, 5, 6, 4, 5, 2, 1)
library(BHH2)
dotPlot(x, pch = 16, xlab = "")
library(plotrix)# thanks to Peter
dotp
aegea wrote:
>
>
> Can R make an usual dotplot just like Minitab and other softwares?
>
> I have the following data, and can use dotchart to graph a dotplot:
>
>
>
> But I don't want the dots on different level of y-axis since y-axis has no
> meaning here.
>
y = rnorm(10)
plot(y,rep(
rkevinburton wrote:
>
> I am testing 'qr' with an admittedly contrived matrix and I am getting
> different results than I am from another package. The matrix that I am
> using is:
>
> x <- matrix(seq(.1, by=.1, length.out=12), 4)
>
> So the whole test is:
>
> x <- matrix(seq(.1, by=.1, length
I agree with Bill's advice, but if you want the easy way out,
try dotplot.mtb in package plotrix. Jim Lemon's done the job
for us.
-Peter Ehlers
On 2010-02-17 23:23, bill.venab...@csiro.au wrote:
R can't do anything.
The question is whether you can do it with R. R puts you in the driver seat
R can't do anything.
The question is whether you can do it with R. R puts you in the driver seat
and if necessary gives you the mechanic's tookkit so that you can tinker
underneath the bonnet if need be.
If the function you like isn't quite there, have a look at the code and adjust
it. It's
Exact to the degree that real numbers can be exactly expressed. Look into
?uniroot
You want to specify the function in a one-sided way, i.e., x/(1-exp(-x))-2.2
HTH,
Daniel
-
cuncta stricte discussurus
-
-Original Message-
From: r-help-boun
Dear R experts,
Can R make an usual dotplot just like Minitab and other softwares?
I have the following data, and can use dotchart to graph a dotplot:
y=c(2.873438152e-01, -8.732895642e-01,
4.579001889e-01, 1.047395204e+00,
8.491182299e-02 , -1.938007105e+00,
-1.273708343e+00, 9.8
Good catch, Dennis. Thanks for permission copy this to the list to set
the record straight.
-- David.
On Feb 17, 2010, at 9:34 PM, Dennis Murphy wrote:
Hi David:
I don't know if anyone corrected you on your use of bquote in the
title below,
but
plot(1,1)
title(main= bquote(ss~"is sum of
Rhett Harrison wrote:
Hi,
I have been having problems installing the newest version on Linux
(Ubuntu 9.10) (tried on two machines).
The ./configure appears to work but I get the following error on the
'make' command.
Don't know about this, but if there's no particular reason you need to
com
Frank E Harrell Jr wrote:
William Dunlap wrote:
Both plyr and rms contain an object called ".".
In plyr it is a "closure" (the common kind of
function) and in rms is is NA. If plyr is attached
in front of rms then you get your problem with
Predict().
Bill Dunlap
Spotfire, TIBCO Software
wdunlap
Hi,
I have been having problems installing the newest version on Linux
(Ubuntu 9.10) (tried on two machines).
The ./configure appears to work but I get the following error on the
'make' command.
$ make
make: *** No targets specified and no makefile found. Stop.
Below I'll post the contents of
Hi:
Alternatively, you could use rollmean in the zoo package:
library(zoo)
x <- c(2, 4, 5, 5, 6, 4, 5, 2, 1)
rollmean(x, 4)
[1] 4.00 5.00 5.00 5.00 4.25 3.00
which avoids the NAs that get produced from filter:
filter(x, rep(1/4, 4))
Time Series:
Start = 1
End = 9
Frequency = 1
[1] NA 4.00 5.00
Erik Iverson wrote:
Frank E Harrell Jr wrote:
Cody,
How amazing that SAS is still used to produce reports that reviewers
hate and that requires tedious low-level programming. R + LaTeX has
it all over that approach IMHO. We have used that combination very
successfully for several data and
Frank E Harrell Jr wrote:
Cody,
How amazing that SAS is still used to produce reports that reviewers
hate and that requires tedious low-level programming. R + LaTeX has
it all over that approach IMHO. We have used that combination very
successfully for several data and safety monitoring rep
Hello!
I want to use R to calculate the variable x which is in some eqation, give
an example:
2.2=x/(1-exp(-x)),
how to solve equation to get the exact x in R?
Thank you.
Have a nice day.
--
View this message in context:
http://n4.nabble.com/Solve-equation-using-R-tp1559576p1559576.html
Dear all,
I have some questions for count data.
First, I want to calculate the expectation value and variance from truncated
data when I assume that the truncated data follow Poisson or negative
binomial distribution.
For example,
N freq
1 30
2 35
3 26
4 8
5
Dear R users,
Overdispersion is often a problem in binomial data. I attempt to model a
binary response (sex-ratio) with three categorical explanatory variables,
using GLM, which could assume the form:
y<-cbind(sexf, sample-sexf)
model<-glm(y ~ age+month+year, binomial)
summary(model)
Output:
(D
Sorry Just a generic list
Is<-list()
forgot to add that from my actual code
From: jim holtman [mailto:jholt...@gmail.com]
Sent: Wednesday, February 17, 2010 3:58 PM
To: ROLL Josh F
Cc: r-help@r-project.org
Subject: Re: [R] Procedure not working for actual data
Yo
Jim,
Thanks for the input so far. Perhaps a solution would be to do what im doing
below, differently, and trying it with my actual data set and seeing if i have
the same issues. To be honest, i am kind of a novice R programmer and the
below took me a substantial amount of time to construct.
Agreed. The tables in the pdf the poster at
http://www.stata.com/statalist/archive/2008-01/msg00098.html links to
look terrible compared to the standard I am used to from
Hmisc::latex(). Just saying.
-Ista
On Wed, Feb 17, 2010 at 9:33 PM, Frank E Harrell Jr
wrote:
> Cody,
>
> How amazing that SA
Cody,
How amazing that SAS is still used to produce reports that reviewers
hate and that requires tedious low-level programming. R + LaTeX has it
all over that approach IMHO. We have used that combination very
successfully for several data and safety monitoring reporting tasks for
clinical
Try this on your real data:
> #Sample data
> Bldgid<-c(1000,1000,1001,1002,1003,1003)
> Maplot<-c(2,20001,3,30001,4,40001)
> Area<-c(40,170,50,100,100,4.9)
> #Construct Sample dataframe
> MultiLotBldgs..<-data.frame(Bldgid,Maplot,Area)
> #Get Building Area Proportions
> MultiLotBldgs..
This implies that you have indeed a time series. You cannot run a survival
model on a single unit of obvservation (i.e., one population) unless you can do
the things (or similar things) that David suggested to create a larger dataset
by disaggregating information. However, your initial approach
You might want to put some debugging statements to see what the values are.
As you say the example seems to work, so look at the data:
> MultiLotBldgs2.. <- MultiLotBldgs..[Is$NotTooSmall.X, ]
> # some debugging statements
> str(MultiLotBldgs2..)
'data.frame': 4 obs. of 3 variables:
$ Bldgid:
Yes, this can be easily computed analytically (even though my result is a bit
different).
--- On Fri, 12/2/10, dav...@rhotrading.com wrote:
> From: dav...@rhotrading.com
> Subject: Re: [R] Integral of function of dnorm
> To: "Greg Snow" , "Trafim Vanishek"
> , "Peter Dalgaard"
> Cc: r-help@r
Your example does not work since "Is" is not defined. What is it supposed
to be?
On Wed, Feb 17, 2010 at 6:34 PM, LCOG1 wrote:
>
> Hello all,
> I have what i feel is a unique situation which may not be resolved with
> this inquiry. I have constructed the below data set so that i may give an
Hello all,
I have what i feel is a unique situation which may not be resolved with
this inquiry. I have constructed the below data set so that i may give an
example of what im doing. The example works perfectly and i have no issues
with it. My problem arises with my actual data, which inclu
Hello,
Are there any packages/functions available for testing the assumptions
underlying assumptions for a good GLM model? Like linktest in STATA
and smilar. If not, could somebody please describe their work process
when they check the validity of a logit/probit model?
Regards,
Jay
_
On Feb 17, 2010, at 5:32 PM, FishR wrote:
We are looking the extinction of a species of freshwater fish. The
logistic
regression was derived by scoring the anecdotal descriptions of the
species'
former population size (1 for a positive description of the
population e.g.
abundant, and 0 f
Yes, you were right. Although I have browsed the FAQ I did not connect the
dots ...
Many thanks
Ed.
-Original Message-
From: David Winsemius [mailto:dwinsem...@comcast.net]
Sent: Wednesday, February 17, 2010 8:22 PM
To: Eduardo M. A. M.Mendes
Cc: r-help@r-project.org
Subject: Re: [
We are looking the extinction of a species of freshwater fish. The logistic
regression was derived by scoring the anecdotal descriptions of the species'
former population size (1 for a positive description of the population e.g.
abundant, and 0 for a negative description e.g. scarce) and plotting
On Feb 17, 2010, at 5:07 PM, Eduardo M. A. M.Mendes wrote:
Hello
I wonder whether someone can tell me what I am doing wrong. Here is
the
code (from Bayesian Computation with R - Chapter 2.3.R) that I am
trying to
run
library(lattice)
windows
Tena koe Erik
Like this boxplot(rnorm(100), las=2)?
las is an argument of par() which sets the graphics parameters. It is
well worth the effort to become familiar with par and its arguments.
HTH
Peter Alspach
> -Original Message-
> From: r-help-boun...@r-project.org
> [mailto:r
Hello
I wonder whether someone can tell me what I am doing wrong. Here is the
code (from Bayesian Computation with R - Chapter 2.3.R) that I am trying to
run
# Section 2.3 Using a Discrete Prior
graphics.off(
==> CRASH INSTALLING PACKAGES CHRON IN:
Intel(R) Core(TM) i7 CPU 860 @ 2800GHz
==> IT'S COMPILED DIRECTIVE WRONG? REGARD YOUR ATENTION
==> I RUN:
SO:Linux 2.6.31.12-174.2.19.fc12.x86_64 x86_64 Fedora release 12
(Constantine)
==> AND RECEIVE NEXT MESSAGENS IN INSTALLING:
> install.packages()
L
In the function boxplot: Is it possible to have the numbers on the vertical
axix written horizontally (so that tilting of the head is unneccessary)?
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://stat.ethz.ch/m
Dear all,
There have been a variety of discussions on the R list regarding the use of R
in clinical trials. The following post from the STATA list provides an
interesting opinion regarding why SAS remains so popular in this arena:
http://www.stata.com/statalist/archive/2008-01/msg00098.html
Re
I ran the follow code for an ordered logit, but don't know why two levels of my
dependent variable are at the topic of my list of variables.
I don't know why this appears, and what I'm supposed to take from them
y>=0. Haven't thought much about this
y>=1. Favor
library(Design)
two <- lrm(
On Wed, 17 Feb 2010, Randall Wrong wrote:
Dear R users,
I have multi-argument returns in a function and I am warned by the program
they are deprecated.
Defunct as from the next R release.
I have found this in the R-help archives :
in 2001!
http://tolstoy.newcastle.edu.au/R/help/01c/0319
Hi,
On Wed, Feb 17, 2010 at 4:09 PM, Randall Wrong wrote:
> Dear R users,
>
> I have multi-argument returns in a function and I am warned by the program
> they are deprecated.
>
> I have found this in the R-help archives :
>
> http://tolstoy.newcastle.edu.au/R/help/01c/0319.html
> http://tolstoy.
Hi all,
I'm trying to make a little script to determine an "unknown" rate for a
number of known exponential trials.
My Code:
#Set Trials and generate number
trials=100
rand<-runif(1,0,1)
vector=0
#Generate vector of 100 random exponentials and sum them
for (i in 1:100) {
vector<-rexp(trials
Dear R users,
I have multi-argument returns in a function and I am warned by the program
they are deprecated.
I have found this in the R-help archives :
http://tolstoy.newcastle.edu.au/R/help/01c/0319.html
http://tolstoy.newcastle.edu.au/R/help/01c/0356.html
Since I am not too good at programmi
R 2.10.0
Windows XP
Please forgive a question that is not strictly speaking R centric
I have 24 raters, each of who characterizes 100 pictures into one of three
categories, A, B, or C. Each rater rates each picture twice in random order
with sufficient time between ratings so that the rater is
Hi All,
I have a function that is used with data frames having multiple id's
per row and it aggregates the data down to 1 id per row. It also
randomly selects one of the within-id values of a variable (mod),
which often differ within-id. Assume this data frame (below) is much
larger and I want to
There a numerous issues, some of which David has pointed out. I will add some
and address some:
1. As far as I understand, you look at only one population. For a survival
model, you would need an indicator when the species was extinguished (rather
than a probability). However, with only one ext
Everyone:
If you work in Higher Education, in either an office of Institutional
Research or an office Institutional Effectiveness, I'm keen to hear from you
on how you use R, individually or with colleagues.
Why do you use R? Is the cost savings an issue? Do you also maintain a
license to
Try this:
> tmp <- as.function(list(x=NULL, D( expression( x^3 ), 'x' ) ) )
> curve(tmp)
--
Gregory (Greg) L. Snow Ph.D.
Statistical Data Center
Intermountain Healthcare
greg.s...@imail.org
801.408.8111
> -Original Message-
> From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-
On Feb 17, 2010, at 1:54 PM, FishR wrote:
Dear all
I have a dataset examining the probability of a population surviving
(calculated from a logistic regression) of a species over a 200yr
period.
The predictor variables are either continuous but non-normal (e.g.
temperature, oxygen) or catego
Hello,
I'm interested in doing Bayesian kriging using R. I see that the
package geoR has a function that will allow one to do this
(krige.bayes). However, my data are not in the form of points, but
rather they are blocks that represent spatial averages (i.e., the number
of fishing hooks
Dear all
I have a dataset examining the probability of a population surviving
(calculated from a logistic regression) of a species over a 200yr period.
The predictor variables are either continuous but non-normal (e.g.
temperature, oxygen) or categorical (e.g. channelisation), unfortunately I
also
rle is extremely cool!!
This is really an impressive method to reduce vectors to point where you
need to do calculations. This really helped! Thanks a lot!
jorgusch
--
View this message in context:
http://n4.nabble.com/delete-repeated-values-not-unique-tp1557625p1559171.html
Sent from the R he
Hi Lu,
Just store the id number as a variable at the top of your script, and
modify the title() call as shown below.
score.id <- "003"
ID <-
c("001","001","002","002","002","003","003","003","004","004","005","005","006")
test <-
c("A","B","A","B","C","A","B","C","A","B","A","B","A")
score <- c(
I have data that comes into R already ordered. When I use ggplot, it orders
them which I don't want. How do I fix this without changing
options("contrast")?
The data I have is number of days:
30
29
...
20
19
...
10
9
...
1
When I plot with ggplot, it orders them by the first number only. So 3 en
Yes,
I paste the wrong ERROR MESSAGE... Cause I was checking many options it the
problem was about files with separate names...
But the problem still persists...
I am using the following code:
file<-read.table("/Users/Márcio/UdeM/Travail Dirigé/Data/MU284
Population.txt",header=T,skip=24)
And ge
I'm trying to run an ordered logistic regression model. I've run the following
code, but the output does not provide the p-values. Is there some command to
include the p-values in the output.
reg2 <- polr(trade1 ~ age2 + education2 + personal2 + economy2 + partisan2 +
employment2 + union2 +
Thanks a lot, Brian!
--
View this message in context:
http://n4.nabble.com/Problems-with-boxplot-in-ggplot2-qplot-tp1555338p1558810.html
Sent from the R help mailing list archive at Nabble.com.
__
R-help@r-project.org mailing list
https://stat.ethz.ch
Hi,
I use lapply with a function from a package, and every 20th to 50th function
execution (out of 500 to 10 000 times) fails for some unknown reason. RAM
consumption ever increases (out of limit) and execution time is endless. I am
not going to debug the package. (I rather feel like an end-use
On Feb 17, 2010, at 2:13 PM, Terrence Ireland wrote:
TestHistogram_A is a function that plots a histogram,
It creates a plot-description but does not "plot" it.
then returns 3. For
me it fails to display the plot, at least in the order given below
TestHistogram_B is a function that plots
TestHistogram_A is a function that plots a histogram, then returns 3. For
me it fails
to display the plot, at least in the order given below.
TestHistogram_B is a function that plots a histogram, then exits. The only
difference between
the 2 functions is that B does not return 3. It d
Hello,
I have trouble with the following: I am trying do display two groups in a
wireframe using two different color schemes. I have been able to display the
two groups, but the color scheme for both is identical. I would like to
display one group in a color scheme red-yellow (or whatever) and the
On 18/02/2010, at 3:25 AM, Christine SINOQUET wrote:
> I want to unsuscribe from the R list.
> Thanks.
Go right ahead!!! Don't let us stop you!
cheers,
Rolf Turner
##
Attention:\ This e-mail message
On Feb 17, 2010, at 1:30 PM, David Winsemius wrote:
On Feb 17, 2010, at 1:06 PM, Lu Wang wrote:
Hi,
I want to produce a pie chart with titles. I used the following
code. So it created a pie chart for the percentage of the score for
ID="002".
ID <-
c
("001
","001
","002","002","002
On Feb 17, 2010, at 1:48 PM, Gabor Grothendieck wrote:
Try this:
do.call(curve, list(D(expression(x^3), "x")))
This also worked, after noting that str() told me the D() object had a
"language" attribute, but given that it was not the chosen method by
Gabor, it may have some deficiency ab
Your results show that you have a rank 2 matrix.
> $rank
> [1] 2
Therefore the third column of Q for your x is arbitrary.
A better example would be
x <- matrix(rnorm(12), 4, 3)
qr(x)
Rich
__
R-help@r-project.org mailing list
https://stat.ethz.ch/m
Now available online is our March R / S-PLUS courses schedule in New
York City and other USA cities
http://www.xlsolutions-corp.com/Rcourses
*** S-PLUS / R : Programming Essentials
*** R Fundamentals and Programming Techniques
*** R/S-PLUS Functions by Example.
*** S/R-PLUS Programming 3: Adv
Try this:
do.call(curve, list(D(expression(x^3), "x")))
On Wed, Feb 17, 2010 at 1:31 PM, Walmes Marques Zeviani
wrote:
>
> Hi all,
>
> I want to provide the result from D() to curve(), because I want to plot the
> k-th derivative of some functions. Actually, I copy from console the result
> g
On Feb 17, 2010, at 1:06 PM, Lu Wang wrote:
Hi,
I want to produce a pie chart with titles. I used the following
code. So it created a pie chart for the percentage of the score for
ID="002".
ID <-
c
("001
","001
","002","002","002","003","003","003","004","004","005","005","006")
tes
Hi all,
I want to provide the result from D() to curve(), because I want to plot the
k-th derivative of some functions. Actually, I copy from console the result
given by D() and paste inside curve(). With a lot of functions and high degree
differentiation this process is tedious. Can I provide
Hi everyone,
I'm trying to perform the following regressions in order to compare linear
vs. sigmoidal fit of the relationship between my dependent variable (y) and
one explaining parameter (x2), both including the confounding effects of a
third variable (x1):
quasi-pois-lin <- glm(y ~ x1 + x2, fa
Hi,
I want to produce a pie chart with titles. I used the following code. So it
created a pie chart for the percentage of the score for ID="002".
ID <-
c("001","001","002","002","002","003","003","003","004","004","005","005","006")
test <-
c("A","B","A","B","C","A","B","C","A","B","A","B","A"
Hello,
I am implementing the BoxCox transformation for multiple regression
using the function boxcox() from library MASS (this seems to be the
best one). To avoid negative values in the variable to be transformed
however, one must specify a shift parameter. However different
choices for the valu
Dear R users,
Can someone explain me the interpretation of empirical p-value (
max(t)) used to test the genetic association with PLINK?
Any help will be appreciated
Regards
M
--
Mohamed Lajnef,IE
INSERM U955 eq 15
Pôle de Psychiatrie
Hôpital CHENEVIER
40, rue Mesly
94010 CRETEIL Cedex
Flana gmail.com> writes:
>
> Hi,
>
> First, thank you all for your help.
>
> Here is my problem (simplified):
>
> Say I have a list:
> a=list(matrix(50,nrow=5,ncol=5),
> matrix(25,nrow=5,ncol=5),
> matrix(10,nrow=5,ncol=5))
>
> I'd like to use rbinom with a different probability for each m
It looks like I had to reconfigure with Readline=yes.
When I installed R I had to set Readline to something and I picked "no"
because I had no idea what Readline was.
Thanks
- Original Message -
From: "Erik Iverson"
To: "Steven Martin"
Cc:
Sent: Tuesday, February 16, 2010 8:59 PM
you mean something along the lines of
filter(x, rep(1/4, 4))
(which you can combine with na.omit)
?
On Wed, Feb 17, 2010 at 5:18 PM, Mohsen Jafarikia wrote:
> Hello All:
>
> If I do have:
>
> x = (2, 4, 5, 5, 6, 4, 5, 2, 1)
> y = (9, 11.5, 12.5, 13, 14, 19, 20, 21, 22)
>
> I wanted to find a s
William Dunlap wrote:
Both plyr and rms contain an object called ".".
In plyr it is a "closure" (the common kind of
function) and in rms is is NA. If plyr is attached
in front of rms then you get your problem with
Predict().
Bill Dunlap
Spotfire, TIBCO Software
wdunlap tibco.com
I think I
Hi,
I am trying to test whether a series is return series stationary, but
before proceeding I wanted to make sure I understand correctly how to
use the adf.test function and interpret its output... Could you please
let me know whether I am correct in my interpretations?
ex: I take x such as I kno
Hi everyone,
I've tried your various suggestions and through a combination of them
I think I've managed to make it work. This latest version appears to
be much more finicky than earlier ones. Thank you so much for all your
help!
Best to all,
Jessica
> On Wed, Feb 17, 2010 at 6:32 AM, Ista Zahn w
Hello All:
If I do have:
x = (2, 4, 5, 5, 6, 4, 5, 2, 1)
y = (9, 11.5, 12.5, 13, 14, 19, 20, 21, 22)
I wanted to find a simple function in R which calculates the averages of X
and Y for every 4 unit increase of Y. The results should look like:
x = (4, 6, 3) #where (2+ 4+ 5+ 5)/ 4
Christine SINOQUET wrote:
I want to unsuscribe from the R list.
then visit this website: https://stat.ethz.ch/mailman/listinfo/r-help
Best,
Dimitris
Thanks.
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE
I am testing 'qr' with an admittedly contrived matrix and I am getting
different results than I am from another package. The matrix that I am using is:
x <- matrix(seq(.1, by=.1, length.out=12), 4)
So the whole test is:
x <- matrix(seq(.1, by=.1, length.out=12), 4)
qr(x)
And the output from 'R
Does anyone know how to dock graphics windows in Eclipse when using the
StatET plug-in for R?
Right now every time I make a graph it pops up as a separate window, which
takes up too much real estate.
By the way, if you have not tried the StatET thing, you should. It's really
nice, and I'm sure I
David, this is very nice. Simple and straightforward. Thank you!
-Ista
On Wed, Feb 17, 2010 at 10:09 AM, David Winsemius
wrote:
>
> On Feb 17, 2010, at 8:55 AM, David Winsemius wrote:
>
>>
>> On Feb 17, 2010, at 8:26 AM, Ista Zahn wrote:
>>
>>> Hi all,
>>> I am designing a study in which I am sel
On Wed, 17 Feb 2010, S. Few wrote:
Currently using R 2.92
Win XP
Goal: scale up computing capacity for large datasets (1-5 million records)
I realize under 32 bit versions of R the memory.limit is maxed at 4GB.
Q:
1. What are the limits under 64 bit versions of R? Are those limits OS
depende
Hi,
First, thank you all for your help.
Here is my problem (simplified):
Say I have a list:
a=list(matrix(50,nrow=5,ncol=5),matrix(25,nrow=5,ncol=5),matrix(10,nrow=5,ncol=5))
I'd like to use rbinom with a different probability for each matrix. I
tried:
b=c(.8,.1,.9)
brep=rep(b,each=25)
lapply(
All,
I have searched the previous help boards and discovered the problem
with Rdonlp2 - Specifically, its non-availability. I thought that this
was my solution, but perhaps there is a better way that you all could
help me with. I imagine that this problem is trivial to people such
as the
Currently using R 2.92
Win XP
Goal: scale up computing capacity for large datasets (1-5 million records)
I realize under 32 bit versions of R the memory.limit is maxed at 4GB.
Q:
1. What are the limits under 64 bit versions of R? Are those limits OS
dependent?
2. Are there limits to the size of
Hi Bill, thanks for your response. There are many cases in which I find direct
use of `names<-`, `mode<-`, `body<-`, `[<-`, `[[<-`, `$<-`, etc. to be
helpful... so wrapping each function individually seems a bit verbose. Thanks
for pointing out the compatibility with S-Plus - that was originally
Hi,
Usually species richness data have very heterogeneous variance, so
they are better described by a Poisson distribution than by a
Gaussian. you should probably analyse your data with Generalized
Linear Models using function lmer (from package lme4) and specify
family=poisson.
hope this helps
Please look at the maiz example, the last example in ?MMC
in the HH package, to see how to display pairwiuse contrasts in
split-plot experiments.
Look at the entire example, all the way to the end of the help file.
Rich
__
R-help@r-project.org mailing l
On Feb 17, 2010, at 8:55 AM, David Winsemius wrote:
On Feb 17, 2010, at 8:26 AM, Ista Zahn wrote:
Hi all,
I am designing a study in which I am selecting a subset of college
courses to be randomly assigned to one of two conditions. I would
like
to create matched pairs of courses, and then
Thank you for your recommendations and the submitted code.
That is excactly what I was searching for.
I am working on it to adapt it to my data set.
Best regards and thanks a lot
Johannes
--
View this message in context:
http://n4.nabble.com/Triangular-filled-contour-plot-tp1557386p1558760.html
There is also the OptimizeR project on R-forge
http://r-forge.r-project.org/R/?group_id=395. Other related projects are there also, but
I'll let their authors speak for themselves. Stefan Theussl did a good job in the Task
View, but it is from last June, and it would be a monumental effort to ke
Thank you for the tip. I was used to inserting write statements and was
surpised when it didn't work and reading this section I see that I shouldn't
have been doing this anyway.
One more question. Is there another call that I can use to print out a
2-dimensional array? Since FORTRAN stores as c
I want to unsuscribe from the R list.
Thanks.
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained,
Hi all,
I'd like to compare species richness patterns of two groups with paired
observations with an extra continuous fixed effect as a covariate. Normally
I would use a standard general mixed effect model for such data, with one
group variable and a continuous variable as fixed effects, plus a r
Rob,
There have been a few reports of issues with the default value of 100 for
rows_at_time and I have had issues on OSX connecting to an Oracle server. I
have set the value to 1 in the wrapper functions that I have written to obtain
data from our server where I call odbcConnect() and sqlQuery(
No no...
This is not in R code. I mean that when in the R, I use arules
library, and try to create a association rules for this transaction
file:
a,c,f,3,4,5
b,e,1,2,4
a,c,e,f,1,3,4,5
d,5
b,c,e,f,1,2,3,4
a,c,e,f,1,3,4,5
b,c,e,f,1,3,4
b,e,1,2,4
a,c,e,f,1,3,4,5
a,b,c,e,f,1,3,4
a,c,d,f,3,4,5
I
On Feb 17, 2010, at 8:26 AM, Ista Zahn wrote:
Hi all,
I am designing a study in which I am selecting a subset of college
courses to be randomly assigned to one of two conditions. I would like
to create matched pairs of courses, and then randomly assign them to
condition within each pair. I woul
1 - 100 of 125 matches
Mail list logo
