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
Deepayan Sarkar 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
__
R-help@r
?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 sim
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.
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:
http://www.nabb
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
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 a
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 mess
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"
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 on
On 10/15/08, Dieter Menne <[EMAIL PROTECTED]> wrote:
> Gabor Grothendieck 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 answer the
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, fixed
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]
>
Gabor Grothendieck 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 http://www.R-project.org/post
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 PROTECT
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
batchfile
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 http://www.mayin.org/ajays
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(mymo
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 t
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() o
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
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 a
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())
>
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 structur
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.
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"
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
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/
Greg Snow 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 lattice par
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 certai
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 th
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 t
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 ma
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
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
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
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
> r
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
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 ge
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
differ
Meesters, Erik 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
library(lattice)
show.sett
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 p1>0,
p2>0,
p3>0, p1>p3.
2) The code is giving the following warning.
Warning: optimization did not converge (code 1)
?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:
>
>scale_y_con
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 sho
Benoit Boulinguiez 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
__
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,
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 see
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 colu
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
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 c
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
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 becaus
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
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
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 p
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
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
l
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=
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 i
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 wo
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 overly elaborate,
can anyone suggest a better / shorter / neater way?
The below doe
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 <- sa
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 (try5.
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 Hub
Dear friends
There is a list of arrays comprising different no of rows and columns even
sometimes NULL, such as [[2]] given below. How can we ignore [[2]] or others
like this in the complete list. Any help in this regard is needed. Thanks
[[1]]
[,1] [,2]
[1,]31
[2,]3
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
veld
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 appli
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
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 t
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 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. I
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=
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 k
Alternatively if a number of such operations will be
done and memory is not an issue, then you could do:
singlefoo <- do.call('rbind', foo)
max(singlefoo[, 1])
Patrick Burns
[EMAIL PROTECTED]
+44 (0)20 8525 0696
http://www.burns-stat.com
(home of S Poetry and "A Guide for the Unwilling S User")
Dear Katrin,
Store the "old" data in a dataframe instead of vectors and supply the
name of that dataframe to your model. eg
Model <- glm(species ~ temp + prec + elev, data = your.data.frame,
family = binomial(link = logit))
Notice that I've added some spaces which makes your code more easy to
r
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 in
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 in a data frame, with on
Dear R-helpers,
I would appreicate if someone can help me on the transfer parameter in ARIMAX
and also see what I am doing is correct.
I am using ARIMAX with 2 Exogeneous Variables and 10 years data are as follows:
DepVar Period, depVar, IndepVar1 Period, indepVar1, IndepVar2 Period, indepVar2
101 - 178 of 178 matches
Mail list logo