Hi All,
I have questions about MANOVA which I am still not sure if appropriately I
should use it.
For example I have a data set like this:
BloodPressure (BP) Weight Height
120115165
125145198
15699 176
I know that BloodPressure is correlated with both
You only have one response variable, so MANOVA is not appropriate. One
option would be to compare BP ~ Weight + Height with BP ~ 1. That would
give you a joint test of weight and height together. Since they are
collinear, that should tell you the overall effect of size. There are
other options,
ONKELINX, Thierry Thierry.ONKELINX at inbo.be writes:
I have managed to create a figure on the screen with multiple plots in
it. Something like the example below. When I save that with ggsave(),
only the last plot gets saved (pPath in the example) instead of the
entire figure. Any suggestions
Hello all,
I am using the e1071 SVM with the tune options for classification, which work
pretty well, given the examples of using tune.svm function for classification.
But I have not found any example to tune the SVM novelty detection
(one-classification) parameters (gamma, cost, nu), for
Hi All,
act_2
DateDtime Hour Min Second Rep
51 2006-02-22 14:52:18 14 52 18 useractivity_act
52 2006-02-22 14:52:18 14 52 18 4
55 2006-02-22 14:52:49 14 52 49 4
57 2006-02-22 14:52:51 14 52 51
Dear all,
Thanks for all your responses. It looks like I was misinterpreting the
capabilities of ggsave(). That can only handle the basic ggplot2 plots
(which can be very elaborate). This behaviour is decribed at the end of
chapter 8 in the ggplot2 book, but I missed that.
Opening and closing a
Hi All,
act_2
DateDtime Hour Min Second Rep
51 2006-02-22 14:52:18 14 52 18 useractivity_act
52 2006-02-22 14:52:18 14 52 18 4
55 2006-02-22 14:52:49 14 52 49 4
57 2006-02-22 14:52:51 14 52 51
Dajiang Liu wrote:
Dear All,
I am writing an R library myself. The package build ok, and can be installed. But when I load the library into R using library('xyz'), then I type ?xyz for help in one of the functions named xyz as well, the following error occured,
No documentation for 'xyz' in
Hi All,
act_2
DateDtime Hour Min Second Rep
51 2006-02-22 14:52:18 14 52 18 useractivity_act
52 2006-02-22 14:52:18 14 52 18 4
55 2006-02-22 14:52:49 14 52 49 4
57 2006-02-22 14:52:51 14 52 51
Hi All,
act_2
DateDtime Hour Min Second Rep
51 2006-02-22 14:52:18 14 52 18 useractivity_act
52 2006-02-22 14:52:18 14 52 18 4
55 2006-02-22 14:52:49 14 52 49 4
57 2006-02-22 14:52:51 14 52 51
Hi All,
act_2
DateDtime Hour Min Second Rep
51 2006-02-22 14:52:18 14 52 18 useractivity_act
52 2006-02-22 14:52:18 14 52 18 4
55 2006-02-22 14:52:49 14 52 49 4
57 2006-02-22 14:52:51 14 52 51
On 11/03/2009 11:35 PM, Oliver wrote:
hi, all
I am looking at R package RGL to draw a colored mesh/surface plot like
this one (from matlab).
http://www.mathworks.com/access/helpdesk/help/techdoc/visualize/cbar.gif
The key features I am looking for is surfaced with grid and color, but
not the
Hi everybody.
I'm just learning to use R. I've got a strange problem and I couldn't find any
answer to my problem. So here am I.
I want to visualize a model of regression on repeated measures :
5 measures (y) by x-value (lx for log10(x) ; lxf colonn is for as.factor) ; 5
x-values.
I fail to
Hi,
I'm working with custom slides(Cy5) and working in the normalization of the
arrays.
I have three arrays (technical replicates).
I have sucesfully normalized the data using vsn, however i would like to
compare it to the normalization using spike in controls.
My controls are annotated as
Hi Tammy,
You can do:
act_2temp - act_2[which(!duplicated(act_2[,c(1:5)])),]
act_2temp[1:5,]
act_2 - act_2temp
rm(act_2temp)
Atenciosamente,
Leandro Lins Marino
Centro de Avaliação
Fundação CESGRANRIO
Rua Santa Alexandrina, 1011 - 2º andar
Rio de Janeiro, RJ - CEP: 20261-903
R (21) 2103-9600
I assume the problem is to only keep those rows whose
Rep value is not equal to the Rep value in the prior row.
In that case just compare c(, Rep[-nr]) to Rep and
use the resulting vector, ix, to select out rows.
Rep - as.character(act_2$Rep) could
be simplified to Rep - act_2$Rep if Rep
is
Thanks a lot for your reply.
I expect a more flexible solution. Because I couldn't check that column one by
one.
If another Rep colum is:
useractivity_act
2
2
3
3
3
4
5
5
6
useractivity_idle
i should be able to write a program to make it as:
useractivity_act
2
3
4
5
6
useractivity_idle
Hello dear R-help members,
I recently became interested in using biglm with leaps, and found myself
somewhat confused as to how to use the two together, in different settings.
I couldn't find any example codes for the leaps() package (except for in the
help file, and the examples there are not
Sorry,
act_2temp - act_2[which(!duplicated(act_2[,c(1: 6 )])),] # not 5
Atenciosamente,
Leandro Lins Marino
Centro de Avaliação
Fundação CESGRANRIO
Rua Santa Alexandrina, 1011 - 2º andar
Rio de Janeiro, RJ - CEP: 20261-903
R (21) 2103-9600 R.:236
( (21) 8777-7907
( lean...@cesgranrio.org.br
Its not clear what more flexible refers to. The posted solution
already produces the sample output you showed. What else
should it do?
On Thu, Mar 12, 2009 at 8:15 AM, Tammy Ma metal_lical...@live.com wrote:
Thanks a lot for your reply.
I expect a more flexible solution. Because I couldn't
Dear useRs,
I'm trying to write a loop to sum my data in the following way:
(the second - the first) + (the third - the second) + (the fourth - the third)
+ ...
for each column.
So, I wrote something like this:
c - list()
for(i in 1:ncol(mydata)) {
for(j in 2:nrow(mydata)) {
c[[i]] -
Why use a loop? Try using diff()
x - c(4, 19, 21, 45, 50, 73, 78, 83, 87, 94)
sum(diff(x))
-Original Message-
From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On
Behalf Of Rafael Moral
Sent: Thursday, March 12, 2009 9:04 AM
To: r-help@r-project.org
Subject: [R]
is your data a data frame or a matrix? do you want to compute the
differences columnwise, i.e., for each column independently? consider
this example:
# generate and display dummy data
(d = as.data.frame(replicate(3, sample(5
# compute successive differences columnwise
Hi,
Try this;
lapply( mydata, function(x){
sum( diff( x ) )
} )
Romain
Rafael Moral wrote:
Dear useRs,
I'm trying to write a loop to sum my data in the following way:
(the second - the first) + (the third - the second) + (the fourth - the third)
+ ...
for each column.
So, I wrote
I'm trying to write a loop to sum my data in the following way:
(the second - the first) + (the third - the second) + (the fourth -
the third) + ...
for each column.
This is just sum(diff(x)), or even x[length(x)] - x[1].
Regards,
Richie.
Mathematical Sciences Unit
HSL
Well actually, what about that (Assuming mydata is a data frame)
tail( mydata, 1 ) - head( mydata, 1)
since:
(the second - the first) + (the third - the second) + (the fourth - the third)
= the last - the first
Romain
Rafael Moral wrote:
Dear useRs,
I'm trying to write a loop to sum my
This is a telescoping sum that can be calculated analytically as:
(a[2] - a[1]) + ... + (a[n] - a[n-1]) = a[n] - a[1]
On Thu, Mar 12, 2009 at 9:04 AM, Rafael Moral
rafa_moral2...@yahoo.com.br wrote:
Dear useRs,
I'm trying to write a loop to sum my data in the following way:
(the second -
Dear Rafael,
Perhaps:
sum(diff(x))
where x is your vector. To apply above to your data set (by rows), you could
use
apply(mydata,1,function(x) sum(diff(x)))
See ?diff, ?sum and ?apply for more information.
HTH,
Jorge
On Thu, Mar 12, 2009 at 9:04 AM, Rafael Moral
On 12 Mar 2009, at 13:22, richard.cot...@hsl.gov.uk wrote:
I'm trying to write a loop to sum my data in the following way:
(the second - the first) + (the third - the second) + (the fourth -
the third) + ...
for each column.
This is just sum(diff(x)), or even x[length(x)] - x[1].
I think
Dear R experts,
I have a folder names 200209 and in this folder there are many data files,
such as:
BA020902.txt,BA020903.txt,BA020904.txt,
BA020905.txt,BA020906.txt,BA020909.txt,
BA020910.txt,BA020911.txt,BA020912.txt,
BA020913.txt,BA020916.txt,BA020917.txt,
Hi all,
I would like to know how to create n equidistant random points along a
multisegment line/object of class psp and to retrieve the coordinates of
those generated random points.
thanks
P.L.
--
View this message in context:
Hi:
I can't install the R-2.8.1 on the machine IBM AIX according the
instruction configuratioin(OBJECT_MODE=64). I can successfully
./configure ***, and get a Makefile. But during compiling the source code,
there is some problem I can't fix it. Please tell me some tips about how to
Hi all and thanks for your time in advance,
I can't figure out why summary.lmrob complains when lmrob is used on a
zooreg object. If the zooreg object is converted to vector before
calling lmrob, no problems appear.
Let me clarify this with an example:
library(robustbase)
library(zoo)
Sorry for not answering long. After working a lot with new version of bind
function I find it actually really useful. Just to avoid any possible
hard-to-debug errors I try to stick to some good practices, like always
specifying default values for formal parameters and calling bind with
named
I need a package that can compute missing values of n-dimensional vectors for n
2. This is a kind of interpolation, complicated in dimensions higher than 2.
The idea is that I have a set of fully specified vectors (i.e., with no missing
values) and I get a new vector that has one or more
Dear List,
I have just started to familiarize with artificial neural networks and I find
two names of packages by chance using R search. But unfortunately I can not
find these packages on the CRAN. Is it possible that they were renamed? The
names of the packages in question are:
netalg
One can't assume an arbitrary function necessarily works with zoo
objects although many do; however, for functions that are sufficiently
similar to lm, like lmrob, dyn in the dyn package can transform it to
work with zoo (as well as with ts and a number of other time series
classes). Just preface
Hi,
I have the following R-code:
pdf(file=test.pdf)
hist(rcauchy(100))
dev.off()
How can I set the size of the histogram in my pdf?
I tried to do with par(...) but I didn't find the right parameter yet.
Thanks, Max
--
View this message in context:
I am estimating daily electricity prices using GARCH (1,1). What I would like
to see is whether there is some kind of daily or seasonal effect in variance of
the price series. For instance, variance of electricity prices might be
different (higher) during weekdays as opposed to during weekend.
This is a consequence of changes made to support multiple imputation, where
extracting the standard errors was needed. I will try to fix it. A work-around
is to use lapply() and subset()
lapply(c(E,M,H), function(s) svyquantile(~api99,subset(dclus1, stype==s),
ci=TRUE,quantiles=0.5))
?pdf reads: width, height the width and height of the graphics region in
inches. The default values are 7.
ir. Thierry Onkelinx
Instituut voor natuur- en bosonderzoek / Research Institute for Nature
and Forest
Cel
If you run the example from ?biglm
data(trees)
ff-log(Volume)~log(Girth)+log(Height)
chunk1-trees[1:10,]
chunk2-trees[11:20,]
chunk3-trees[21:31,]
a - biglm(ff,chunk1)
a - update(a,chunk2)
a - update(a,chunk3)
summary(a)
you can then do
b -regsubsets(a, method=forward)
summary(b)
to get the
Set which size? Of the page? - then you need to specify in pdf().
Of the plot margins? - then par$mai
Of something else? Then you need to be more specific about your
intention.
Sarah
On Thu, Mar 12, 2009 at 10:05 AM, Maxl18 maximilianwoe...@hotmail.com wrote:
Hi,
I have the following R-code:
Another method would be to use a summary that incorporates both as a
measure of obesity, In medical investigations it is common to use the
BMI which is the ratio of (weight in Kg) to (height in meters squared).
Yet a third method would be to investigate for nonlinearity on the
response
Much simpler would be:
act2[!duplicated(act2$Rep),] #use the negation of the duplicated
function on $Rep and indexing
DateDtime Hour Min Second Rep
51 2006-02-22 14:52:18 14 52 18 useractivity_act
52 2006-02-22 14:52:18 14 52 18 4
See ?pdf and its width and height arguments.
Also if you want to have the graph centered on a standard page, there are
additional arguments to help you achieve that effect:
pdf(test.pdf, height=5, width=5, paper=letter, pagecentre=TRUE)
hist(rcauchy(100))
dev.off()
-Christos
-Original
Hi,
see
?dev.off
but I think (guessing you use a windows system), you would be better of
in using win.metafile() instead of x11() in initiating the graph, see
?win.metafile
There is no need to show all 100+ graphs on your display if you actually
want them in files.
hth.
BaKaLeGuM schrieb:
Model simplification in lme
Pinheiro and Bates suggest that model simplification of fixed effects in lme
with anova (model,model2)
is anticonservative. What is the best way to simplify a model in lme?
Testing the p value for each fixed
factor with anova (model, type=marginal) and removing it if it
Hi,
I would like to some excel files with some worksheets. I tried this with
the following R script:
library(gdata)
i-1
rc-0
while(rc != try-error) {
wksh-try(read.xls(cluster-microarray-FW.xls,sheet=i,verbose=TRUE,perl=perl))
rc-class(wksh)
print(sprintf(--- i=%2d rc=%s
Here is a program that will get a list of the sheet names and then
read in the indicated name. Modify to suit.
library(rcom)
xls - C:\\test.xls
oxl - comCreateObject('Excel.Application')
comSetProperty(oxl, Visible, TRUE)
owb - comGetProperty(oxl, Workbooks)
ob - comInvoke(owb, Open, xls)
If I understand your question and if I understood Dr Harrell's
description of the Hmisc/Design packages capabilities, this should all
be available. I have been blessed with a dataset so large and with so
relatively few missing values that I have not yet needed to put his
methods into
I am trying to print graphs as pdf's or eps in a for loop, but I can't seem
to get it right
Either it prints only a single eps graph (overwrites) or when I use
#pdf(paste(i,.pdf, sep=)) .. it prints all pdf's but they are empty
Can someone suggest which method shud I use for such.
-Original Message-
From: Rolf Turner [mailto:r.tur...@auckland.ac.nz]
Sent: Wednesday, 11 March, 2009 5:09 PM
To: Szumiloski, John
Cc: r-help@r-project.org
Subject: Re: [R] Easy Recall to get ls(..., all.names=TRUE)?
On 12/03/2009, at 3:16 AM, Szumiloski, John wrote:
Dear useRs,
Dear R help,
This seems to be a commonly asked question and I am able to run examples that
have been proposed, but I can't seems to get this to work with my own data.
Reproducible code is below. Thank you in advance for any help you can provide.
The main problem is that I can not get the
On Wed, 11 Mar 2009, Camarda, Carlo Giovanni wrote:
Dear R-users,
I am searching to the best way to compute a series of n matrix
multiplications between each matrix (mXm) in an array (mXmXn), and each
column of a matrix (mXn).
Please find below an example with four possible solutions.
The
Hi,
your example is quite messy (neither reproducible or minimal). I think
you could try the following,
mdf - data.frame(1:3)
names(mdf) - 147
i - 147
mdf[ as.character(i) ]
Hope this helps,
baptiste
On 11 Mar 2009, at 22:34, Mohan Singh wrote:
Hi everyone
I am trying to use
Sean,
would like to receive expert opinion to avoid potential trouble
[..]
i think the following is the most secure way if one really
really has to do assignment in a function call
f({a=3})
and if one keeps this convention, - can be dropped altogether.
secure is relative, since due to
On Mar 12, 2009, at 11:14 AM, Michael Denslow wrote:
Dear R help,
This seems to be a commonly asked question and I am able to run
examples that have been proposed, but I can't seems to get this to
work with my own data. Reproducible code is below. Thank you in
advance for any help you
On Thu, Mar 12, 2009 at 03:06:30PM +, Mohan Singh wrote:
I am trying to print graphs as pdf's or eps in a for loop, but I can't seem
to get it right
Either it prints only a single eps graph (overwrites) or when I use
#pdf(paste(i,.pdf, sep=)) .. it prints all pdf's but they are
Hi everyone,
I have data from an experiment in which human participants were
instructed to generate a random sequence of yes/no answers under 4
different conditions. I want to test how successful they were in doing
this. More specifically, I want to test the null hypothesis that the 4
conditions
On Mar 12, 2009, at 11:45 AM, David Winsemius wrote:
On Mar 12, 2009, at 11:14 AM, Michael Denslow wrote:
# I am not sure how to get the lines to draw correctly here
matlines(pred.frame$a,pp, lty=c(1,2,2),col=black)
The x values are your sequence whereas the y values are in the
Hello everyone,
Do we have package in R which will estimates the parameter in
Seemingly Unrelated Spatial Regression Model? Any help would be appreciated.
Thanks
Arun
--
View this message in context:
If these are the only files in the directory, then you might try...
(File paths will need to change if the folder 200209 isn't in the
working directory)
fpath - ./200209
a - list.files(fpath)
for(i in 1:length(a)){
assign(paste(y,i,sep=),read.table(paste(fpath,a[i],sep=/)))
}
Another option
Dear all,
I did this replacement below and its work for the first case but not in the
second case! Please tell me what's wrong with the second case.
Case 1:
m
[,1] [,2]
[1,]14
[2,]25
[3,]36
l
[,1] [,2]
[1,]9 10
matric1-replace(m,m[3,],l)
Hi, All
I have so many files in my directory. I got the following error. At this time,
I couldn't use debug() to check error one by one because so many files getting
invloved. How can I get some idea about my error??
.
Skipping book_log-20060222T143612.txt
Skipping
Thank you all for the suggestions and answers.
Oliver
On Mar 12, 6:51 am, Duncan Murdoch murd...@stats.uwo.ca wrote:
On 11/03/2009 11:35 PM, Oliver wrote:
hi, all
I am looking at R package RGL to draw a colored mesh/surface plot like
this one (from matlab).
Tammy Ma wrote:
Hi, All
I have so many files in my directory. I got the following error. At this time,
I couldn't use debug() to check error one by one because so many files getting
invloved. How can I get some idea about my error??
The usual advice: simplify the problem until it's
Dear all,
Is there a function to fit nested logit available somewhere ?
Sincerly
Justin BEM
BP 1917 Yaoundé
Tél (237) 99597295
(237) 22040246
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
Hello. I have a script in which I repeatedly fit a nonlinear regression to
a series of data sets using nls and the port algorithm from within a loop.
The general structure of the loop is:
for(i in 1:n){
extract relevant vectors of dependent and independent variables
estimate starting
Hi Bill,
You should have a look at the try() function that lets you enclose an
expression and handle failure in a way that prevents the entire execution
from halting.
Josh
On Thu, Mar 12, 2009 at 1:13 PM, Bill Shipley
bill.ship...@usherbrooke.cawrote:
Hello. I have a script in which I
If the goal here is to repeatedly run a particular model from
different starting values
then the nls2 package will do that for you.
On Thu, Mar 12, 2009 at 2:13 PM, Bill Shipley
bill.ship...@usherbrooke.ca wrote:
Hello. I have a script in which I repeatedly fit a nonlinear regression to
a
?try
For example,
for (i in 1:n) {
try (fit - nls(...), silent=TRUE)
if (class(fit) != try-error) dowhateverthatneedstobedonewiththeresults
else fit - NA
}
Ravi.
---
Ravi Varadhan, Ph.D.
Assistant
Hi all,
I'd like to do cross-validation on lm and get the resulting lift curve/table
(or, alternatively, the estimates on 100% of my data with which I can get
lift).
If such a thing doesn't exist, could it be derived using cv.lm, or would we
need to start from scratch?
Thanks!
--
Eric Siegel,
I think I answered a very similar question from you yesterday
but perhaps the mail went astray. The subject line is not
informative.
It may make it easier to think about if you use a function like
isFirstInRun - function(x) c(TRUE, x[-1]!=x[-length(x)]
Given a vector x (without NA's in it) it
Simpler, but maybe wrong. Not duplicated
was my first response as well, but then I began
wondering if the question implied globally
duplicated or duplicated within subgroups.
Patrick Burns
patr...@burns-stat.com
+44 (0)20 8525 0696
http://www.burns-stat.com
(home of The R Inferno and A Guide
Thank you for you help Dr. Winsemius.
The problem seems to stem from the fact that I have used the incorrect name in
the prediction dataframe.
The following code seems to work correctly.
Thank you again,
Michael
wt.data - data.frame(code = factor(LETTERS[1:24]),
area =
I would like to use grep and gsub to manipulate a vector to make the names used
consistent, i.e. reduce a level or two.
However, here is what I found when I attempted to use grep and gsub:
tmp_test-c(House 1 Plot Plus +100,House 2 Plot Plus +100,House 3 Plot
Plus -100,House 4 Plot Plus
True, I suppose. The specification from Ma could have been more
explicit. Should she want to only have duplicates in sequence, then
perhaps the use of construction within the rle function would be useful.
act2[ act2$Rep[-1] != act2$Rep[-length(act2$Rep)], ]
DateDtime Hour Min
Hi,
After hearing good things about the R community over on the west coast, I
have started a New York based R meetup group. We intend to have our first
get together on Thursday April 2nd, at 7pm. For more details pop on over to
http://meetup.com/nyhackr
Regards,
Josh
--
Joshua Reich
646 256
Hi,
Im using the lm() function where the formula is quite big (300 arguments) and
the data is a frame of 3000 values.
This is running in a loop where in each step the formula is reduced by one
argument, and the lm command is called again (to check which arguments are
useful) .
This takes
Hi,
Im using the lm() function where the formula is quite big (300 arguments)
and the data is a frame of 3000 values.
This is running in a loop where in each step the formula is reduced by one
argument, and the lm command is called again (to check which arguments are
useful) .
This takes 1-2
Jason Rupert wrote:
I would like to use grep and gsub to manipulate a vector to make the names used consistent, i.e. reduce a level or two.
This is dangerous: grep and gsub use regular expressions, so a lot of
characters (see ?regexp for the list) have special meanings.
For the kind of
You should try installing fGarch package. Then you can simulate various garch
models using garchsim command. Hope it help.
jurica
SUMANTA BASAK-4 wrote:
Hi All,
I,m trying to do a GARCH simulation in R 2.3.0 release
in Windows XP. I've seen garchsim function but that is
for garch
yes, indeed, you can certainly speed things up, by just changing the
design matrix X and feeding it back to lm.fit().
In addition, if you just need the least squares estimates, then you gain
a bit more by using constructs of the form:
XtX - crossprod(X)
Xty - crossprod(X, y)
betas -
Look at:
?update
For example:
lm.obj - lm (y ~ x1 + ... + x300)
lm.obj1 - update(lm.obj, . ~ . - x1)
lm.obj2 - update(lm.obj1, . ~ . - x2)
Ravi.
Ravi Varadhan, Ph.D.
Assistant Professor,
Division of Geriatric Medicine
On Mar 12, 2009, at 3:02 PM, David Winsemius wrote:
True, I suppose. The specification from Ma could have been more
explicit. Should she want to only have duplicates in sequence, then
perhaps the use of construction within the rle function would be
useful.
act2[ act2$Rep[-1] !=
I think you will find that many readers of this list would rather try
to dissuade you from this misguided strategy. You are unlikely to get
to a sensible solution in using step-down procedures with this sort of
situation (large number of predictors with modest size of data).
--
David
On Thu, 12 Mar 2009, Duncan Murdoch wrote:
Jason Rupert wrote:
I would like to use grep and gsub to manipulate a vector to make the names
used consistent, i.e. reduce a level or two.
This is dangerous: grep and gsub use regular expressions, so a lot of
characters (see ?regexp for the
Hi,
Is there a package I can use to write to multiple sheets on xls file,
other than using the paid version of xlsReadWrite package
(xlsReadWritePro)?
Thank you,
Ferry
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
Unfortunately, the suggested way with JauntePE does not work completely.
Maybe I have found a better alternative now:
There is a little free application by Sysinternals (now owned by MS) which
can create junction points, i.e. a system link to a directory under Windows.
It is available here:
Jens Oehlschlägel wrote:
Sean,
would like to receive expert opinion to avoid potential trouble
[..]
i think the following is the most secure way if one really
really has to do assignment in a function call
f({a=3})
and if one keeps this convention, - can be dropped
I would like to find the maximum values on the density function of a
random variable. For example, I have a random variable
rv - rbinom(1,1,0.1) + rnorm(1)
Its density function is given by density(rv) and can be displayed by
plot(density(rv)). How to calculate its maximum values?
A
I think that Gabor sent something earlier today ( or atleast something
related ) so check the archives for his solution.
On Thu, Mar 12, 2009 at 4:29 PM, Ferry wrote:
Hi,
Is there a package I can use to write to multiple sheets on xls file,
other than using the paid version of
Sueli Rodrigues wrote:
Hi,
I'm studying the Michael J. Crawley book (Statistics An Introduction using
R)and I'm trying to reproduce one example (page 168) for plot a barplot
with error bars, but when I put the function I receive the following
message. Please, could someone told me what's wrong?
I have two data.frames.
One is the subset of the other. I need to determine what rows were removed
from the superset data.frame to form the subset data.frame.
I was hoping to turn up functionality similar to setdiff for data.frames, but
did not have any luck. It seems I need to find a
Jorge was kind enough to point out that I made a mistake in referencing
gabor's solution.
See below. Thanks Jorge.
Begin forwarded message
Subject: Re: [R] Writing xls - multiple sheets
Date: 3/12/2009 6:03:13 PM
From: Jorge Ivan Velez
To: markle...@verizon.net
Hi Mark,
Try this:
library(sqldf)
BODsub - BOD[1:3,]
sqldf(select * from BOD except select * from BODsub)
For more see: http://sqldf.googlecode.com
On Thu, Mar 12, 2009 at 5:58 PM, Jason Rupert jasonkrup...@yahoo.com wrote:
I have two data.frames.
One is the subset of the other. I need to determine
That was to read worksheet names but its true that the same packages
rcom and rdcom can control Excel and therefore write spreadsheets.
On Thu, Mar 12, 2009 at 5:52 PM, markle...@verizon.net wrote:
I think that Gabor sent something earlier today ( or atleast something
related ) so check the
I get the following error when I run qplot()
qplot(grade, read,data = hhm.long.m, geom = c(point, smooth))
Error in smooth.construct.cr.smooth.spec(object, data, knots) :
x has insufficient unique values to support 10 knots: reduce k.
I am not sure how to tackle this problem. When I take a
Hello All,
Does anyone know of a package that performs constraint-based clusters?
Ideally the package could perform Time-Ordered Clustering, a technique
applied in a recent journal article by Runger, Nelson, Harnish (using MS
Excel). Quote, in our specific implementation of constrained
1 - 100 of 113 matches
Mail list logo