K. Fleischer wrote:
Dear all,
I have built glm models based on presences/absences and a number of
predictor maps and would like to compute habitat suitability based on
the modelled coefficients.
I thought this is pretty straight forward and wanted to use predict()
and supply the new data
On Wed, 15 Oct 2008, Göran Broström wrote:
In the R for Windows FAQ (2.7.2) we have under 2.2:
There is no specific version for x64 Windows, but the standard 32-bit
version works well enough.
Does that include the handling of huge data sets?
Yes, but 'huge' is not the largest size in the
Dear R developers,
i'm developing a Java application with a very efficient image transfer api
between ImageJ and R.
In my next release i can transfer bytes to R very fast with the help of the
Rserve application.
Furthermore i added an easy to use interface to the clustering algorithm
clara in
Gad Abraham wrote:
This approach leaves much to be desired. I hope that its
practitioners start gauging it by the mean squared error of predicted
probabilities.
Is the logic here is that low MSE of predicted probabilities equals a
better calibrated model? What about discrimination? Perfect
An alternative approach would be to store 0 x 0 matrices instead of
NULLs. This way every object in your list is a consistent type.
Hadley
On Wed, Oct 15, 2008 at 5:23 AM, Muhammad Azam [EMAIL PROTECTED] wrote:
Dear friends
There is a list of arrays comprising different no of rows and columns
Benoit Boulinguiez benoit.boulinguiez at ensc-rennes.fr writes:
I'd like to know whether R is capable to assess parameters in a model
describing the kinetic of a pollutant adsorption onto activated carbon.
Start with the deSolve package (especially demo(CCL4model))
Ben Bolker
Dears,
I'm trying to find the parameters (a,b, ... l) that optimize the function
(Model)
described below.
1) How can I set some constraints with MLE2 function? I want to set p10,
p20,
p30, p1p3.
2) The code is giving the following warning.
Warning: optimization did not converge (code 1)
Meesters, Erik Erik.Meesters at wur.nl writes:
is there a way to define the color of the title for the legend in
lattice?
Getting the right par to set in lattice can be intimidating. I keep the result
of trellis.par.get() in a text file and search for the closest match.
Dieter
Hi everyone,
I have a dataset(named Mydata) which includes 4 different variables named;
s1,s2,s3,s4 .Each variable(symptom) has 14 patients.
I need to use random sampling to make, 5 different samples from my data with
5 patients in each sample. i.e. using all 4 variables I need to make 5
You did not read Peters response close enough.
Do:
names(areadata)
And you will see that the names of your areadata object do not match with what
the predict function is expecting. When you create areadata, you need to name
the columns (or name them afterwards), it does not automatically
Hello,
Writing 19000 lines of R-script in Excel sounds terrible.
Could you provide us some code (please NOT 19000 lines), and the way you
generate this with Excel,
maybe we can figure out a much shorter/faster way to accomplish the same
result?
Bart
Duncan Murdoch-2 wrote:
On 10/14/2008
Dear all,
thanx for the quick reply but your ideas have not gotten me further
unfortuantely..
after implementing the idea from Thierry my code looks like that now:
#b00all is the dataframe with only element used for model fitting
#variables .._zui contain the whole are over which to predict
maybe:
result - list()
for (i in 1:10) result[[i]] - rnorm(i)
Bart
Megh Dal wrote:
Can anyone please tell me how to define a list. Suppose I want to define
a list object result with length n then want to fill each place of
result with different objects. For e.g.
i=1
result[1] =
On Wed, Oct 15, 2008 at 5:59 AM, Michael Pearmain [EMAIL PROTECTED] wrote:
Hi All,
I've been trying to compare if the previous value in a variable is equal to
a binary value..(i.e i want to check if the last event was a yes or no)
i've been trying to write some code for this, but it seems
Alex99 wrote:
Hi everyone,
I have a dataset(named Mydata) which includes 4 different variables named;
s1,s2,s3,s4 .Each variable(symptom) has 14 patients.
I need to use random sampling to make, 5 different samples from my data with
5 patients in each sample. i.e. using all 4 variables I need
Dear Alex,
Is this what you want?
# Data set
my=read.table(textConnection(
X8 X9 X10 X102 X110 X177 X283 X284 X286 X292 X297 X306 X308 X314
0 1 000100000000
0 0 001000000010
0 1 00000
Try
x[ !sapply(x, is.null) ]
hadley wickham wrote:
An alternative approach would be to store 0 x 0 matrices instead of
NULLs. This way every object in your list is a consistent type.
Hadley
On Wed, Oct 15, 2008 at 5:23 AM, Muhammad Azam [EMAIL PROTECTED] wrote:
Dear friends
There
Hello,
I would like to extract the standard error for the coefficients returned
when using survreg. Does anyone know how to do this?
Thank you.
Joan
--
View this message in context:
http://www.nabble.com/Extracting-standard-error-from-survreg--tp19994742p19994742.html
Sent from the R help
Thanks for the reply Peter,
that works BUT the reason I didn't use it, is because I need to calculate
the mean and Standard deviation for S1,S2,S3 and S4 in each sample. which
means I'll have 5 means for S1, and 5 means for S2 and 5 means for S3 and 5
means for S4 (and at the end I have to get
Run Rprof on your script that is updating the dataframe. A dataframe
is a list and everytime you access something in the list it can be
expensive. Rprof will probably show that a lot of time is spent in
the function [[ which is accessing portions of the dataframe.
Vectors are much faster because
Thanks for the reply Jorge,
that works BUT the reason I didn't use it, is because I need to calculate
the mean and Standard deviation for S1,S2,S3 and S4 in each sample. which
means I'll have 5 means for S1, and 5 means for S2 and 5 means for S3 and 5
means for S4 (and at the end I have to get
Another way to look through the possible par settings that hopefully makes it a
little less intimidating is:
library(TeachingDemos)
TkListView(trellis.par.get())
This requires the TeachingDemos package (obviously) and the tcltk package (most
have it, but sometimes you need to run R a certain
Greg Snow Greg.Snow at imail.org writes:
Another way to look through the possible par settings that hopefully makes it
a little less intimidating is (completed by DM):
library(lattice)
library(TeachingDemos)
TkListView(trellis.par.get())
Nice; should come handy in other contexts. But with
jpl wrote:
Hello,
I would like to extract the standard error for the coefficients returned
when using survreg. Does anyone know how to do this?
Thank you.
Joan
Apparently, it is
summary(mymodel)$table[,2]
(Why, oh why, can't Terry allow coef(summary(...)) like we have in
lm/glm???
Dear Thomas Lumley, and R-help list members,
I have read your article Network meta-analysis for indirect treatment
comparisons (Statist Med, 2002) with great interest. I found it very
helpful that you included the R code to replicate your analysis;
however, I have had a problem replicating your
Hey,
My dataset has 1 dependent variable(Logloss) and 7 independent dummy
variables(AS,AM,CB,CF,RB,RBR,TS) , it's attached in this email. The problem
is I cant finish Khmaladze test because there's an error Error in
switch(mode(x), NULL = structure(NULL, class = formula), : invalid
formula
Yes, I agree that for searching for specific terms, the file approach is great.
I suggested the TkListView more for browsing through the possible parameters
to learn what is there and can be changed (and to see what the current values
are of those of interest).
--
Gregory (Greg) L. Snow Ph.D.
On Wed, 15 Oct 2008, Peter Dalgaard wrote:
jpl wrote:
Hello,
I would like to extract the standard error for the coefficients
returned when using survreg. Does anyone know how to do this?
Apparently, it is
summary(mymodel)$table[,2]
In order not to compute on the internal structure
?mathplot
but this should work:
qplot(1:10,1:10)+scale_y_continuous(expression(Respiration, pmol *O[2]*h^-1))
On Wed, Oct 15, 2008 at 7:45 AM, Adam Marsh [EMAIL PROTECTED] wrote:
How are sub/superscripts designated for text in axis labels? As in
this y-axis label:
Also:
str(trellis.par.get())
show it compactly.
On Wed, Oct 15, 2008 at 11:18 AM, Greg Snow [EMAIL PROTECTED] wrote:
Another way to look through the possible par settings that hopefully makes it
a little less intimidating is:
library(TeachingDemos)
TkListView(trellis.par.get())
This
Hi:
After a couple of days trying to synchonize or connect MiKTex and R I'm getting
a little frustrated. My knowledge about creating and running batch
files is kind of limited so, I was wondering if someone could explain in plain
english how to do this task. I'm able to run MiKtex by itself and
Dear Alex,
Is this what you want?
my=read.table(textConnection(
X8 X9 X10 X102 X110 X177 X283 X284 X286 X292 X297 X306 X308 X314
0 1 000100000000
0 0 001000000010
0 1 0000000
Please forgive this repost, it's been a week without a squeak. No
comments?
Original post:
https://stat.ethz.ch/pipermail/r-help/2008-October/176340.html
Hello,
Two lapply questions (system info and sample data below):
1) Why does the first form of command1 add the name of y _after_ the
str()
Ruben,
Thankyou for the advice. I'll do what I can with it.
Ruben Wrote:
If I understand your problem correctly, you have that the magnitude of
deviations from the mean/median/mode in the volume of your requests for
background checks in month m predicts a multivariate response that
represents
Can anyone please tell me how to define a list. Suppose I want to define a
list object result with length n then want to fill each place of result
with different objects. For e.g.
i=1
result[1] = rnorm(1)
i=2
result[2] = rnorm(2)
...
i=n
result[n] = rnorm(n)
What would
Thanks! This is exactly what I needed.
Peter Dalgaard wrote:
jpl wrote:
Hello,
I would like to extract the standard error for the coefficients returned
when using survreg. Does anyone know how to do this?
Thank you.
Joan
Apparently, it is
summary(mymodel)$table[,2]
I wondered was people on this list felt about this article:
http://www.voxeu.org/index.php?q=node/2363
which talks about the problems of obtaining sound answers in numerical
optimisation in settings such as MLE or NLS.
--
Ajay Shah
Suppose foo.Rnw is in \a\b and you have placed
MiKTeX in the default place when you installed it
and suppose that when you installed R you did
not prevent it from writing to the reigstry. Then:
1. Go to the link I gave already, read everything
there, go to the download page and download
Hey back,
If you had RTFM you would know that the first argument of the
function KhmaladzeTest was supposed to be a formula, so try:
T - KhmaladzeTest(Logloss~AS+AM+CB+CF+RB+RBR, data=CPBP,
nullH=location)
url:www.econ.uiuc.edu/~rogerRoger Koenker
email[EMAIL
Gabor Grothendieck ggrothendieck at gmail.com writes:
str(trellis.par.get())
Never thought of that. Really nice.
Dieter
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide
Hi, I apologise in advance for the naïve question. I have large matrices that I
want to plot. I currently use color2D.matplot. However, these matrices contain
many values of no interest (i.e. where there is no data, the figure -999 is
automatically displayed). Is there any way of removing these
Achim Zeileis wrote:
On Wed, 15 Oct 2008, Peter Dalgaard wrote:
jpl wrote:
Hello,
I would like to extract the standard error for the coefficients
returned when using survreg. Does anyone know how to do this?
Apparently, it is
summary(mymodel)$table[,2]
In order not to compute on
What do you mean by removing them? Do you want the whole row/column
deleted? Have you tried setting them to NA so that the points are not
plotted? You have to be a little more specific on what type of action
you want to have happen when the criteria is met. Of course the
answer is anything is
Hi Christian,
I believe that the argument var has changed name to weights.
The following lines work for me:
...
sigma - rep(sqrt(.5), nrow(lumley1)) # not nrow=
lme1 - lme(Y1 ~ trt.B + trt.C + trt.D + trt.E, random = ~ 1 | trtpair,
data=lumley1, weights = varConstPower(form=~sigma,
On 10/15/08, Dieter Menne [EMAIL PROTECTED] wrote:
Gabor Grothendieck ggrothendieck at gmail.com writes:
str(trellis.par.get())
Never thought of that. Really nice.
There's also
http://lmdvr.r-forge.r-project.org/figures/figures.html?chapter=07;figure=07_03
All of which doesn't really
I am trying to investigate the interaction term in the below. The
paradigm in aquatic systems is that algal production is either
nitrogen (TIN) or Phosphorus limited, and I am trying to investigate
this- what is the best way to go about investigating the interaction
term. I have some thoughts
In the R for Windows FAQ (2.7.2) we have under 2.2:
There is no specific version for x64 Windows, but the standard 32-bit
version works well enough.
Does that include the handling of huge data sets? What I have read in
other places indicates that the answer is No. If so, will there be
a 64-bit
Dear R users,
This is to announce new versions of our packages granova and PSAgraphics, both
now v.1.2.
granova derives from ‘graphical analysis of variance.’ The package consists of
four functions that facilitate seeing basic data as well as standard summary
statistics for various ANOVA
Gad Abraham explained :
Max wrote:
Hi everyone,
This is not so much of an R question as a statistics question. I currently
work for the largest pre employment screening company in Canada. Upper
management has noticed that noticed that usually a month or so before any
big kind of economic
Dear R users,
is there a way to define the color of the title for the legend in
lattice?
The help page on xyplot has a lot of details on key options just as the
new book, but no mentioning of a color attribute for the title.
Should I use ltext or is there any other way?
Best wishes,
Erik
Michael Just wrote:
Hello,
I write most of my R code in excel and then paste it into R. I am wondering
if there is a limit to how much I can paste? I want to paste about 19,000
lines of code should this work? I am doing this because when I did it chunks
it took about an hour and half. I thought
Try this:
lapply(1:n, rnorm)
On Wed, Oct 15, 2008 at 8:19 AM, Megh Dal [EMAIL PROTECTED] wrote:
Can anyone please tell me how to define a list. Suppose I want to define
a list object result with length n then want to fill each place of
result with different objects. For e.g.
i=1
result[1]
Hi,
I'd like to know whether R is capable to assess parameters in a model
describing the kinetic of a pollutant adsorption onto activated carbon.
A common relation is for instance the Adam-Bohart-Thomas' one:
dx/dt = K1 * (qm-x)*C - K2x
where {K1,K2} are the unknown paramters and {qm,C} are
L.S.,
I would like to add a downsized image to a pdf file and got stuck on the
code to use.
My goal is to open a pdf device, add a plot, add the downsized image, and
then close the pdf device.
Making the downsized image is easy using png(),
png(file=x.png, width = 600, height = 600,
Hello R-list members!
I tried to do the following with my dataset that contains factor and
numerics, (80columns,about 600 rows)
Dataset.afdm-AFDM(Dataset[282:595,], type=TypeVector, ncp=3)
Fehler in svd(X) : infinite or missing values in 'x'
TypeVector
[1] n n n n n n n n n n n n n n n n n n
Hello all,
I'm working with a large data-set, and upgraded my RAM to 4GB to help
with the mem use.
I've got a 32bit kernel with 64GB memory support compiled in.
gnome-system-monitor and free both show the full 4GB as being available.
In R I was doing some processing and I got the following
Dear useRs,
I am struggling to use gamboost function form the 'mboost' package. More
precisely, I am trying to extract the *partial fit* for each of the
covariates estimated in a model and I usually end up with this annoying: Error
in newdata[[xname]] : subscript out of bounds . I hope that the
Thanks again Jorge,
but when I type in lapply(res,function(x) rbind(rowMeans(x)) it's like it's
waiting for something else, it doesnt run. I tried help(lapply) to learn
more about lappy but didnt give me any result. could you tell me why it
doesnt run and what is lapply for?
Thanks a lot for
For the first question, ?str tells us that the str function does not return
anything, but has a side effect of printing information to the console. So in
version 1 when you call the cat function it evaluates its arguments and as it
evaluates str(...) the information is printed, then nothing is
Hi everyone,
I have just started using R, and I have a simple question.
How can I get the Standard deviation for rows. basically I am looking for
something like rowMeans()
but for Standard deviation (I tried rowSds() didn't exist)
Thanks,
--
View this message in context:
I want to thank everyone for the help. I ended up having to use a loop to
assign values from the table to NinYear. However, as I have played with the
full datasets I have noticed that R is MUCH faster if I use vectors in the
loop rather than columns of a dataframe. In the specific case of 43,000
Em Ter, 2008-10-14 às 17:13 +0200, Arnau Mir Torres escreveu:
Hello.
I need to know how can R compute AIC when I study a regression model?
For example, if I use these data:
growth tannin
1 12 0
2 10 1
3 8 2
4 11 3
5 6 4
6 7 5
See ?Memory-size
On Wed, 15 Oct 2008, B. Bogart wrote:
Hello all,
I'm working with a large data-set, and upgraded my RAM to 4GB to help
with the mem use.
I've got a 32bit kernel with 64GB memory support compiled in.
gnome-system-monitor and free both show the full 4GB as being available.
?apply
e.g. apply(matrix,1,sd)
-Original Message-
From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED]
On Behalf Of Alex99
Sent: Wednesday, October 15, 2008 1:17 PM
To: r-help@r-project.org
Subject: [R] Standard deviation for rows
Hi everyone,
I have just started using R, and I have a
How are sub/superscripts designated for text in axis labels? As in
this y-axis label:
scale_y_continuous(Respiration, pmol O2 h-1)
where the 2 needs to be subscripted and the -1 needs to be
superscripted.
Thanks.
Adam
---
Adam G.
Deepayan Sarkar deepayan.sarkar at gmail.com writes:
All of which doesn't really answer the original question.
...
So, the title color cannot be currently specified, either directly or
through the settings system.
Krrr.. it was key$title
Dieter
Dear Mahbub
Thanks a lot for the help.
- Original Message
From: Mahbub Latif [EMAIL PROTECTED]
To: Muhammad Azam [EMAIL PROTECTED]
Sent: Wednesday, October 15, 2008 12:30:13 PM
Subject: Re: [R] request: How can we ignore a component of list having no
element
try this,
junk -
Em Ter, 2008-10-14 às 18:44 -0700, Lavan escreveu:
Hi R users,
I'm trying to have the symbols for sigma[1] in my legend. the code is given
below, please have a look.
Thanks,
lavan
try legend(bottom,legend=expression(sigma[1]),...)
legend(bottom,legend=paste(Sigma1=,
Dear Stephen,
I'm not sure exactly what you have in mind, but you might take a look
at the effects package.
I hope this helps,
John
On Wed, 15 Oct 2008 12:59:33 -0400
stephen sefick [EMAIL PROTECTED] wrote:
I am trying to investigate the interaction term in the below. The
paradigm in
This approach leaves much to be desired. I hope that its practitioners
start gauging it by the mean squared error of predicted probabilities.
Is the logic here is that low MSE of predicted probabilities equals a
better calibrated model? What about discrimination? Perfect calibration
implies
Thank you very much Greg.
Apparently, my learning curve is still well below it's horizontal
asymptote.
And no wonder in the vast universe of the R-project, eh?
Thanx, DaveT.
-Original Message-
From: Greg Snow [mailto:[EMAIL PROTECTED]
Sent: October 15, 2008 01:52 PM
To: Thompson, David
Dear all
I am intrested in making iterative estimation (thro' loop statements) of, say,
linear regression model. For this purpose, I have written the following
programme and that I have made use of a sample data (viz., exp.txt):
Programme:
# Linear regression modelling with sample data
On 10/10/2008 8:55 AM, Oliver Bandel wrote:
Zitat von Duncan Murdoch [EMAIL PROTECTED]:
On 10/10/2008 8:13 AM, Oliver Bandel wrote:
Hello,
I tried to use rgl.snapshot and it failed.
The error message was not very verbose:
==
plot3d( motion[[idx+2]],
Doen't work.
\misiek
Prof Brian Ripley wrote:
See ?Memory-size
On Wed, 15 Oct 2008, B. Bogart wrote:
[...]
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do
apply( mymatrix, 1, sd )
Will give the row standard deviations of a matrix or matrix like data frame.
There are also some functions in add-on packages that do row sd's or var's.
--
Gregory (Greg) L. Snow Ph.D.
Statistical Data Center
Intermountain Healthcare
[EMAIL PROTECTED]
801.408.8111
Thanks for the reply,
that's exactly what I wanted. I also figured that the function rowSds() does
exist. I just needed to load the generfilter package.
Thanks again.
Nutter, Benjamin wrote:
?apply
e.g. apply(matrix,1,sd)
-Original Message-
From: [EMAIL PROTECTED]
Hi there,
I am looking for R/parallel package or some other package that would speed
up the analysis.I am working on computatioanly intensive data so any
suggestions would be really helpful.
Kindly let me know if any
--
View this message in context:
i'm in the process of switch from Maple to R and am trying to code the
following function:
w(d)=\int -Inf_Inf A(x) |int 0_(t(d)-x) Fy dydx
can some one point me in the right direction? i don't seem to be able to
figure it out on my own.
Dr. Wade Winterhalter
University of Central Florida
Dear all,
I have been fitting models to do meta-analysis in R. This includes a
mixed effect model with a weighting function. This weighting function
is the sample size and the random factor is study or species for
instance (Nakagawa 2007). I have tried this method and I find some
strange things.
Hi,
have a look to Dirks tutorial at the UseR2008. This should be a good
starting point:
http://www.statistik.uni-dortmund.de/useR-2008/tutorials/eddelbuettel.html
Markus
Rajasekaramya wrote:
Hi there,
I am looking for R/parallel package or some other package that would speed
up the
Or ?Memory-limits (and the posting guide of course).
On Wed, 15 Oct 2008, Prof Brian Ripley wrote:
See ?Memory-size
On Wed, 15 Oct 2008, B. Bogart wrote:
Hello all,
I'm working with a large data-set, and upgraded my RAM to 4GB to help
with the mem use.
I've got a 32bit kernel with 64GB
In the example in the documentation, I see:
rs - dbSendQuery(con,
select Agent, ip_addr, DATA from pseudo_data order by Agent)
out - dbApply(rs, INDEX = Agent,
FUN = function(x, grp) quantile(x$DATA, names=FALSE))
Maybe I am a bit thick, but it took me a while, and a kind
I have examined the documentation for batch mode use of R:
R CMD BATCH [options] infile [outfile]
The documentation for this seems rather spartan.
Running R CMD BATCH --help gives me info on only two options: one for
getting help and the other to get the version. I see, further on, that
Hi,
Is there a way to simulate a population with R and pull out m samples,
each with n values
for calculating m means?
I need that kind of data to plot a graphic, demonstrating the central
limit theorem
and I don't know how to begin.
So, perhaps someone can give me some tips and hints
Look at the clt.examp function in the TeachingDemos package, it generates
samples from normal, uniform, gamma (default exponential), and beta (default
U-shaped) distributions and plots histograms of the means along with a
reference line of a normal distribution with the same mean and sd. The
Hello R-users,
I have code that gives me the important variables from an analysis. I
need to input these variables into a different analysis. To do this, I
need to modify them slightly... 1) remove all numbers at the end of the
variables, 2) remove all periods.
I tried to do it with the
On 10/15/2008 3:48 PM, Jörg Groß wrote:
Hi,
Is there a way to simulate a population with R and pull out m samples,
each with n values
for calculating m means?
I need that kind of data to plot a graphic, demonstrating the central
limit theorem
and I don't know how to begin.
The easiest
Hi all,
I am running OLS models with a power parameter: y=a+bx^p, and I want to iterate
through p with a set increment. The function I constructed is something like
below:
test-function(data, model,...){
f-formula(model) #model=y~I(x^p)
p-0.1
n-10
result-list()
while (p=1) {
for (i in 1:n) {
Hi Joerg,
Is there a way to simulate a population with R and pull out m samples,
each with n values
for calculating m means?
I need that kind of data to plot a graphic, demonstrating the central
limit theorem
and I don't know how to begin.
So, perhaps someone can give me some tips and
A case for [g]sub:
gsub(., , sub([0-9]*$, , pick), fixed = TRUE)
On Wed, 15 Oct 2008, John Poulsen wrote:
Hello R-users,
I have code that gives me the important variables from an analysis. I need
to input these variables into a different analysis. To do this, I need to
modify them
Try this;
gsub(\\.|[0-9], , pick)
On Wed, Oct 15, 2008 at 4:42 PM, John Poulsen [EMAIL PROTECTED] wrote:
Hello R-users,
I have code that gives me the important variables from an analysis. I need
to input these variables into a different analysis. To do this, I need to
modify them
To remove all numbers and periods at the end do:
vars - sub([0-9.]+$, , pick )
If there are other periods that this does not remove (not at the end) then you
can do a second pass:
vars - sub(\\., , vars)
Part of the reason that your code for the period does not work is that an
unescaped
I can't help with your inquiry generally, but I can address the
issue of ``optional braces''.
On 16/10/2008, at 8:34 AM, Ted Byers wrote:
snip
I take it FUN = function(x, grp) quantile(x$DATA, names=FALSE) is the
function definition for a function called FUN. I would guess,
Jörg Groß wrote:
Hi,
Is there a way to simulate a population with R and pull out m samples,
each with n values
for calculating m means?
I need that kind of data to plot a graphic, demonstrating the central
limit theorem
and I don't know how to begin.
So, perhaps someone can give me some
Galton's 19th century mechanical version of this is the quincunx. I
have a
(very primitive) version of this for R at:
http://www.econ.uiuc.edu/~roger/courses/476/routines/quincunx.R
url:www.econ.uiuc.edu/~rogerRoger Koenker
email[EMAIL PROTECTED]
Hello,
I've been learning R functions recently and I've come across a problem that
perhaps someone could help me with.
# I have a a chron() object of times
hours=chron(time=c(01:00:00,18:00:00,13:00:00,10:00:00))
# I would like to subtract 12 hours from each time element, so I created a
The if statement is for program flow control (do this bunch of code if this
condition is true, else do this), the vectorized version is the ifelse
function, that is the one that you want to use.
Also note (this is for times in general, not sure about chron specifically)
subtracting noon from
Joe Kaser wrote:
Hello,
I've been learning R functions recently and I've come across a problem that
perhaps someone could help me with.
# I have a a chron() object of times
hours=chron(time=c(01:00:00,18:00:00,13:00:00,10:00:00))
# I would like to subtract 12 hours from each time element,
Hi Siang Li. It would help if you explained what packages you are
using. auto.arima() is in the forecast package and arimax appears to
be from the TSA package.
Using auto.arima() to select the orders is inappropriate because you
are ignoring the regressors. auto.arima() does not currently handle
Another request for help implementing the 'apply' functions to avoid a
loop structure...
I am working with a data set that includes lab measurements taken at
different dates for the subjects, with some subjects having more
results than others. I would like to average lab results for each
subject
1 - 100 of 163 matches
Mail list logo