Thanks!
I found an other, less elegant, way as well. (This always happens, just
after asking someone else)
a-data.frame(x=c(1:5),y=c(2,2,3,3,4))
require(lattice)
png(filename=let's change the font.png, width=480, height=300,
pointsize=12)
xyplot(y~x,data=a,
ylab=list(y-axis,font=6),
It seems that this addition works, but has created just one object
called 'peak' with all the data from those 100 files. I'd like each
file to have a corresponding object containing the data.
Thanks for your help!
On Tue, Apr 28, 2009 at 19:43, Zeljko Vrba zv...@ifi.uio.no wrote:
On Tue, Apr
On Tue, 28 Apr 2009, K. Elo wrote:
mathallan wrote:
How can I from the summary function, decide which glm (fit1, fit2 or fit3)
fits to data best? I don't know what to look after, so I would please
explain the important output.
Start with the AIC value (Akaike Information Criterion). The
By doing peak - list() you initialize a list and then you are adding the
data objects to the list in the for loop.To access each data object just go:
peaks[[i]]
Where i is in [1,100]. This will return the data object you want.
Adrian
On Tue, Apr 28, 2009 at 11:23 PM, Taylor Hermes
H == HighSchool2005 gemeaux...@yahoo.fr
on Tue, 28 Apr 2009 06:29:07 -0700 (PDT) writes:
H Hello dear R users,
H My objective is to change the size of this graphic : plot(lm(a~b), 4)
H (Cook's distance)
H I have found the help on the internet saying to change the size of
Tena Sakai tsakai at gallo.ucsf.edu writes:
I learned 3 new tricks. (Not bad for a newbie?)
$ R --silent --no-save barebone.R
$ R --quiet --no-save barebone.R
$ R --slave barebone.R
And don't forget
R --vanilla
which I like most because of the taste.
Dieter
carol white wrote:
Hi,
I wanted to plot the histogram of a vector and then, plot the density function
of subsets of the vector on the histogram. So I use truehist in MASS package
and lines(density) as follows:
length(b) = 1000
truehist(b)
lines(density(b[1:100]))
however the density plot of
Hillary Cooper a écrit :
Hi-
I'm trying to make my PCA (princomp) colored. In my csv excel sheet, I have
the first column numbered according to the groupings I want to assign to the
PCA. I've played around with trying to set this first column as the color
vector, but haven't had any luck.
Any
Wacek Kusnierczyk wrote:
Dieter Menne wrote:
R --vanilla
which I like most because of the taste.
once we talk about preferences, here's the version which i like most
because of (a) least typing, (b) the non-pompous little r:
... and (c) superior performance (in this
Dieter Menne wrote:
Tena Sakai tsakai at gallo.ucsf.edu writes:
I learned 3 new tricks. (Not bad for a newbie?)
$ R --silent --no-save barebone.R
$ R --quiet --no-save barebone.R
$ R --slave barebone.R
And don't forget
R --vanilla
which I like most because of the
Bert Gunter gunter.berton at gene.com writes:
Martin's reply provides an appropriate response, so nothing to add. But my
questions dig deeper: Why do so many (presumably nonstatisticians, but ?)
belong to this R^2 religion? Is it because:
1) This is what they are taught in their Stat 101
Based on your solution, could we not fix this issue for the symbol
function?
Le 27 avr. 09 à 23:47, Duncan Murdoch a écrit :
On 27/04/2009 5:23 PM, Christophe Dutang wrote:
so it could be a limitation of graphics on windows?
It's not graphics on Windows, but it appears to be a limitation
Hi,
That's pretty impressive performance, but wait.
Where does this little r come from? And how
does it differ from its big brother?
Regards,
Tena Sakai
tsa...@gallo.ucsf.edu
-Original Message-
From: r-help-boun...@r-project.org on behalf of Wacek Kusnierczyk
Sent: Wed 4/29/2009
Thanks! I have not used system before. May it works if command ist the
call for the program, but how can I integrate the transfer of the file names
(strings) from R to the program?
Ranjan Maitra wrote:
Does ?system help?
Ranjan
On Tue, 28 Apr 2009 13:00:28 -0700 (PDT) thoeb
Hi, I am trying to do an analysis of variance for Randomized-complete blocks
design. I have 4 treatments and 3 replication, without intraction.then I am
going to use LSD test for comparison between medium values of treatments.
I already did that in SAS, but I am not familiar with R very well.
Tena Sakai wrote:
Hi,
That's pretty impressive performance, but wait.
Where does this little r come from? And how
does it differ from its big brother?
the little r comes from littler [1]. it doesn't claim to be larger than
it is.
[1] http://dirk.eddelbuettel.com/code/littler.html
Hi all,
If I have a huge matrix/ dataframe and I want to create a new matrix/
dataframe with every second (or third, or fourth etc.) row of the original
matrix, how can I do it? Any help much appreciated.
Thanks,
Umesh
[[alternative HTML version deleted]]
Try this,
a = matrix(rnorm(10*10),ncol=10)
a[, seq(1,ncol(a),by=2)]
baptiste
On 29 Apr 2009, at 09:28, Umesh Srinivasan wrote:
Hi all,
If I have a huge matrix/ dataframe and I want to create a new matrix/
dataframe with every second (or third, or fourth etc.) row of the
original
Richardson, Patrick wrote:
Uwe,
I deleted gregmisc from by library folder and ran update.packages() again and
I still get the same error:
package 'fBasics' successfully unpacked and MD5 sums checked
Error in unpackPkg(foundpkgs[okp, 2L], foundpkgs[okp, 1L], lib) :
malformed bundle
probably you want to use seq(), e.g.,
mat - matrix(1:500, 100, 5)
mat[seq(1, nrow(mat), 2), ]
mat[seq(1, nrow(mat), 3), ]
mat[seq(1, nrow(mat), 4), ]
I hope it helps.
Best,
Dimitris
Umesh Srinivasan wrote:
Hi all,
If I have a huge matrix/ dataframe and I want to create a new matrix/
Thanks a lot Baptiste Dimitris,
It's working!
Cheers,
Umesh
On Wed, Apr 29, 2009 at 2:06 PM, Dimitris Rizopoulos
d.rizopou...@erasmusmc.nl wrote:
probably you want to use seq(), e.g.,
mat - matrix(1:500, 100, 5)
mat[seq(1, nrow(mat), 2), ]
mat[seq(1, nrow(mat), 3), ]
mat[seq(1,
Thank you very much again. Now I just need to find out how to update R on
Ubuntu with only aptitude available. haha!
Martin Maechler-4 wrote:
H == HighSchool2005 gemeaux...@yahoo.fr
on Tue, 28 Apr 2009 06:29:07 -0700 (PDT) writes:
H Hello dear R users,
H My objective is
Hi.
I have some problem so, I want you to help me.
I have to classify some data, and show the result as plot diagram to my
boss.
I decide to use plsr function in pls package to analysis that.
My data are some similar oliveoli data in that package. so, I make my data
Bendix,
you can also take a source package an run
R CMD Rd2dvi --no-clean packageName
on it and you will get a temporary directory with the TeX sources in it.
Best wishes,
Uwe
BXC (Bendix Carstensen) wrote:
In version 2.8.1, running Rcmd check on the package foo would leave the file
Hi R-users,
I'm not sure what's wrong and how do I check the error?
I have 3 sets of data
aug_dt1
[1] 22.8 49.6 77.9 67.6 91.2 48.9 124.6 54.6 85.8 62.8 63.3 174.7
32.3 22.7 99.9 123.2 4.8 149.0 66.0
[20] 73.9 37.4 22.4 69.5 87.3 66.9 87.2 81.5 109.1 47.4 22.8
Hello Derek,
Unfortunately, I can't help you with R2HTML.
But if you want to generate html output of R commands, you could also try ascii
package http://eusebe.github.com/ascii/ (available on CRAN).
Here an example file (example.Rnw, with asciidoc
http://www.methods.co.nz/asciidocmarkup) :
thoeb wrote:
Hello,
does anybody know about how to send a command or a text line from R to
another program? I have written a script in which several calculations are
made and outputfiles (csv) are generated. Afterwards I open another program
(Fortran) via shell.exec. This program asks for the
Hi,
Consider also that you may easily modify system names from inside your R
script as:
system(mv oldname newname)
Best Regards,
Javier
...
thoeb wrote:
Hello,
does anybody know about how to send a command or a text line from R to
another program? I have written a script in which several
Hi all,
I tried to a nice legend with small boxes filled with the colors used
for the plots. But it does nor work, boxes are always filled with black.
An example is here
plot(1:4,1:4)
lines(1:4,4:1, col=blue)
legend(top,leg=c(a,b),col=c(black,blue), fill=TRUE)
How could I specify the colors?
Christophe Dutang a écrit :
Hi all,
I tried to a nice legend with small boxes filled with the colors used
for the plots. But it does nor work, boxes are always filled with black.
An example is here
plot(1:4,1:4)
lines(1:4,4:1, col=blue)
#try something like this
Jim Lemon wrote:
thoeb wrote:
Hello,
does anybody know about how to send a command or a text line from R to
another program? I have written a script in which several
calculations are
made and outputfiles (csv) are generated. Afterwards I open another
program
(Fortran) via shell.exec. This
Dieter Menne wrote:
Bert Gunter gunter.berton at gene.com writes:
Martin's reply provides an appropriate response, so nothing to add. But my
questions dig deeper: Why do so many (presumably nonstatisticians, but ?)
belong to this R^2 religion? Is it because:
1) This is what they are taught
Christophe Dutang wrote:
Hi all,
I tried to a nice legend with small boxes filled with the colors used
for the plots. But it does nor work, boxes are always filled with black.
An example is here
plot(1:4,1:4)
lines(1:4,4:1, col=blue)
legend(top,leg=c(a,b),col=c(black,blue), fill=TRUE)
How
Hi Christophe,
just try :
legend(top,leg=c(a,b),fill=c(black,blue))
instead of your :
legend(top,leg=c(a,b),col=c(black,blue), fill=TRUE)
Regards, Olivier
--
Olivier ETERRADOSSI
Maître-Assistant
CMGD / Equipe Propriétés Psycho-Sensorielles des Matériaux
Ecole des Mines d'Alès
Hélioparc, 2
Dieter Menne dieter.menne at menne-biomed.de writes:
q(paper) = 10* n(pvalues) + 5*n(R^2) + 3.5*n(Error Bars)
Values above 300 qualify for immediate acceptance, and Journals
like Lancet, New English and British Journal of XXX provide
professional advice.
I noted the and was misleading.
Hi,
I have a comma seperated data file which has blanks in it. I am trinying to
import it to R using
data1-read.csv(oa_2006.csv, header = TRUE, sep = ,, quote=\,
dec=.)
I want the missing values to be NA. instead R reads them as U. Any idea
why this happens ?
Thanks in advance.
Regards,
On Wed, Apr 29, 2009 at 01:58:06PM +0530, Umesh Srinivasan wrote:
Hi all,
If I have a huge matrix/ dataframe and I want to create a new matrix/
dataframe with every second (or third, or fourth etc.) row of the original
matrix, how can I do it? Any help much appreciated.
Others ahve already
ID: T8dfe471e360ac80264c6c
Dear,
I'm using R 2.8.1 with the package reshape version 0.8.2 .
It seems that the combine_factor function has problems with combining levels to
an empty level:
#EXAMPLE WORKING
test-sample(c(1:10),1000,replace=T)
testf-factor(test, c(1:10), LETTERS[1:10])
c - structure(c(2L, 2L, 1L, 3L, 4L, 2L, 3L, 2L, 3L, 2L, 5L), .Label = c(foo,
+ bar, a really really long variable label mostly here to show the need of
word-wrapping text in labels,
+ a not so important value, baz), class = factor)
plot(c)
Is there a way to get the long variable labels to
I thought I have tried to pass a vector of colors and that it did
not work..
Thanks for all!
Christophe
Le 29 avr. 09 à 12:16, Jim Lemon a écrit :
Christophe Dutang wrote:
Hi all,
I tried to a nice legend with small boxes filled with the colors
used for the plots. But it does nor
Christophe Dutang wrote:
Based on your solution, could we not fix this issue for the symbol
function?
Sure, send a patch.
Duncan Murdoch
Le 27 avr. 09 à 23:47, Duncan Murdoch a écrit :
On 27/04/2009 5:23 PM, Christophe Dutang wrote:
so it could be a limitation of graphics on
G'day Uwe,
On Wed, 29 Apr 2009 11:03:43 +0200
Uwe Ligges lig...@statistik.tu-dortmund.de wrote:
you can also take a source package an run
R CMD Rd2dvi --no-clean packageName
on it and you will get a temporary directory with the TeX sources in
it.
Which is fine for manual processing and
Terry Therneau wrote:
It is likely a problem with survival, since
2.9 merged in a large number of changes that had occured in my source
tree that had not propogated to the R tree
my test suite doesn't have a test for this particular case (2 factors)
and -
Hi Hans,
strwrap is your friend. \n inserts a cr in a line of text.
wordwrap-function(x,len) paste(strwrap(x,width=len),collapse=\n)
par(mar=c(11,3,2,1))
tmp-plot(c,axes=F)
axis(2)
axis(1,at=tmp,labels=sapply(levels(c),wordwrap,len=15),padj=1)
box()
For unique abbreviations see ?abbreviate
Hans Ekbrand wrote:
c - structure(c(2L, 2L, 1L, 3L, 4L, 2L, 3L, 2L, 3L, 2L, 5L), .Label = c(foo,
+ bar, a really really long variable label mostly here to show the need of word-wrapping text in labels,
+ a not so important value, baz), class = factor)
plot(c)
Is there a way to get the long
Nattu wrote:
Hi,
I have a comma seperated data file which has blanks in it. I am trinying to
import it to R using
data1-read.csv(oa_2006.csv, header = TRUE, sep = ,, quote=\,
dec=.)
I want the missing values to be NA. instead R reads them as U. Any idea
why this happens ?
Hi Nataraju,
Dieter Menne dieter.me...@menne-biomed.de wrote
I noted the and was misleading. Read: Good journals like Lancet,
New English and many British Journal of XXX really help you to do
better.
I am one of the statistical editors for PLoS Medicine, and I try to help
people do better; often, the
the little r comes from littler [1]. it doesn't claim to be larger
than
it is.
[1] http://dirk.eddelbuettel.com/code/littler.html
Unless you're working on a Mac or Windows (if I'm not mistaken about
case handling there), where it will just overwrite the standard R
interpreter. :-(
Hi Tena,
I recommend rapache for building websites with R. See
http://biostat.mc.vanderbilt.edu/rapache/
-Ista
-- Forwarded message --
From: Tena Sakai tsa...@gallo.ucsf.edu
To: ted.hard...@manchester.ac.uk, r-help@r-project.org
Date: Tue, 28 Apr 2009 15:04:53 -0700
Subject:
Hi,
I tried to change my settings of the mailing-list, so that I get my
own message again from the mailing-list.
But it doesn't work - I don't get a copy of my own email.
Does anyone know how you can change that?
thanks for any help!
__
Dear all,
Is there an R package which calculates the modified Struve function
(StruveL)?
I was not able to find anything apt with a RSiteSearch so it does not
appear to have been asked in the list.
Many thanks and best wishes,
Ranjan
__
It's gmail, not the mailing list. Gmail stores your mailing list messages under
Sent mail and I at least haven't found any way to convince it to
filter them or
store them elsewhere. If you reply to an earlier message, your reply should
appear in the proper place, but there doesn't seem to be any
When you read in a .csv file with some empty cells using read.csv,
there is no need to specify anything - empty cells will be recorded as
NAs automatically:
try:
data1-read.csv(oa_2006.csv)
Dimitri
On Wed, Apr 29, 2009 at 7:54 AM, Jim Lemon j...@bitwrit.com.au wrote:
Nattu wrote:
Hi,
I
Just another opinion about R^2 coming from the field of US Psychology
research and business:
The first and foremost technique taught in Psychology Departments in
subfields where experimental designs are rarely possible (i.e., social
psychology, personality psychology, developmental psychology,
Hi there,
I'm trying to solve a boxplot problem (should be simple, but I cannot find the
solution...):
data - list(long_name1_xxx = rnorm(100),
long_name2_yy = rnorm(200))
boxplot(data , las = 2, mar = c(20,4,4,4))
Why does the margin not change
Thanks to Jim and Eik!
I really appreciate your help, and I think can use your suggestions
and perhaps write a wrapper for plot that integrates them.
--
Hans Ekbrand (http://sociologi.cjb.net) h...@sociologi.cjb.net
Q. What is that strange attachment in this mail?
A. My digital signature, see
Hello:
First, RExcel completely rocks! Hat tip to the development team.
However, I don't understand something very basic about it.
Suppose I put the following code in two cells and choose the Run Code option
plot(runif(20),runif(20))
text(.5,.5,Hello)
I get the expected plot.
HOWEVER, if I
boxplot() isn't passing extra parameters in the way you'd expect. If
you set par() first, it works fine.
mar.orig - par()$mar # save the original values
par(mar = c(20,4,4,4)) # set your new values
boxplot(data, las = 2)
par(mar = mar.orig) # put the original values back
Sarah
On Wed, Apr 29,
Hi friends,
Please help me with this bug.
*Bug in my code:*
In this variable sub_grp_whr_cls_data[sbgrp_no,1] I store the where
clause.every sub group has a where condition linked with it.
Database1
Where clause was not found for a particular subgroup,
sub_grp_whr_cls_data[sbgrp_no,1]
I knew the solution would be kind of simple...
Thanks a lot, Sarah!
Sarah Goslee schrieb:
boxplot() isn't passing extra parameters in the way you'd expect. If
you set par() first, it works fine.
mar.orig - par()$mar # save the original values
par(mar = c(20,4,4,4)) # set your new values
Even if you have the mailing list options set so that you are
sent copies of your own posts, gmail hides them for you.
http://groups.google.com/group/Gmail-Help-Message-Delivery-en/browse_thread/thread/e7716ab04941c383?pli=1
It has nothing to do with the mailing list manager.
The original
I got a distribution function and a empirical distribution function. How do I
make to Kolmogorov-Smirnov test in R.
Lets call the empirical distribution function Fn on [0,1]
and the distribution function F on [0,1]
ks.test( )
thanks for the help
--
View this
Berwin A Turlach wrote:
G'day Uwe,
On Wed, 29 Apr 2009 11:03:43 +0200
Uwe Ligges lig...@statistik.tu-dortmund.de wrote:
you can also take a source package an run
R CMD Rd2dvi --no-clean packageName
on it and you will get a temporary directory with the TeX sources in
it.
Which is fine
see
?is.null
e.g.
if( is.null(sub_grp_whr_cls_data[sbgrp_no, 1]) )
{
your code
}
Moumita Das wrote:
Hi friends,
Please help me with this bug.
*Bug in my code:*
In this variable sub_grp_whr_cls_data[sbgrp_no,1] I store the where
clause.every sub group has a where condition linked with
Hi
you put your problem in absolutely messy state
r-help-boun...@r-project.org napsal dne 29.04.2009 16:16:55:
Hi friends,
Please help me with this bug.
*Bug in my code:*
In this variable sub_grp_whr_cls_data[sbgrp_no,1] I store the where
clause.every sub group has a where condition
help.search(kolmogorov)
?ks.test
andydol...@gmail.com
2009/4/29 mathallan mathanm...@gmail.com
I got a distribution function and a empirical distribution function. How do
I
make to Kolmogorov-Smirnov test in R.
Lets call the empirical distribution function Fn on [0,1]
This is the third homework question you have asked the list to do for you. How
many more should we expect?
The posting guide is pretty clear in that: Basic statistics and classroom
homework: R-help is not intended for these.
-Original Message-
From: r-help-boun...@r-project.org
Hi
r-help-boun...@r-project.org napsal dne 29.04.2009 17:05:01:
see
?is.null
e.g.
if( is.null(sub_grp_whr_cls_data[sbgrp_no, 1]) )
{
your code
}
It probably will not work as
sub_grp_whr_cls_data[sbgrp_no,1]==NULL
implicates that there is character value NULL. I am not sure if you
My program consists of a number of functions and a main script.
It loops through many time series analyzing one at a time.
I have declared a number of arrays and scalar variables in the main script
namespace (using C++ terminology)
because those data structures ar shared by many functions.
test - list(NULL)
Bert Gunter
Nonclinical Biostatistics
467-7374
-Original Message-
From: r-help-boun...@r-project.org
[mailto:r-help-boun...@r-project.org] On Behalf Of Petr PIKAL
Sent: Wednesday, April 29, 2009 8:21 AM
To: ml-r-h...@epigenomics.com
Cc:
Petr PIKAL wrote:
Hi
r-help-boun...@r-project.org napsal dne 29.04.2009 17:05:01:
see
?is.null
e.g.
if( is.null(sub_grp_whr_cls_data[sbgrp_no, 1]) )
{
your code
}
It probably will not work as
sub_grp_whr_cls_data[sbgrp_no,1]==NULL
implicates that there is character value
Tena Sakai wrote:
Hi,
Many thanks to Wacek Kusnierczyk and Ted Harding.
I learned 3 new tricks. (Not bad for a newbie?)
$ R --silent --no-save barebone.R
$ R --quiet --no-save barebone.R
$ R --slave barebone.R
With slight differences, they all do what I wanted.
Moving
Hi all
I seem to have fallen at the first hurdle with my analysis, I have a set
of binary disease outbreak data linked to a large number of landscape
metrics variables and environmental variables which I would like to as
predictor variables in a Least Angle Logistic Regression using the
Dear R People:
I upgraded to Open Suse 11.1 and am having trouble with the rgl
package, which I need for other applications.
Here is the install output:
install.packages(rgl,depen=TRUE)
--- Please select a CRAN mirror for use in this session ---
Loading Tcl/Tk interface ... done
trying URL
Hi R community,
I am trying to obtain a sample from a matrix but sample(my.matrix) doesn't
do what I need. I have a matrix of 1287 interactions between the species in
columns and the species in rows and I want to obtain a smaller matrix with
say, 800 interactions, that may or may not have the
On Tue, 28 Apr 2009, Neotropical bat risk assessments wrote:
Hi again all,
I can read all of the data is in CSV format with a header row and with
1,200,240 rows.
I can do all the plots and ggplot2 has no problems.
Seems to be a problem with the MASS package???
It seems that the
If it is a fixed-effects ANOVA you can just use biglm(). Otherwise
you could use biglm() to fit the necessary sequences of models to give the
RSS that you need for the F-tests.
-thomas
On Tue, 28 Apr 2009, Hardi wrote:
Hi,
I'm using aov() to analyze the data and get the rank of
If we knew two pth quantiles for a normal distribution,
is it possible that we can find mean and sigma for the normal distribution
using R?
Let x ~ norm(mean, sigma).
Suppose that qnorm(0.9,mean,sigma) and qnorm(0.1,mean,sigma) are known.
Can we find mean and sigma using R?
Thanks,
-james
If you just want to shuffle the row around, this should do it. By
limiting the number of samples, you can get a smaller matrix:
pla- c(10, 9, 6, 5, 3) #abundance of pla species
pol- c(14, 10, 9, 4, 2) #abundance of pol species
m-pla%*%t(pol) #matrix of interactions
m[sample(nrow(m)),]
Hi,
I have an experiment with 2 independant factors which I have been trying to
analyse in R. The problem is that there are several data points recorded on
the same animal. However, no combination of treatments is repeated on the
same animal. All possible combinations of treatments are done in
Hello everybody,
using the lmer function, I have fitted the following logistic mixed
regression model on an experimental data set containing one fixed factor
(Cond) and three random variables (Sito, Area, Trans):
model-lmer(Caul~Cond+(1|Sito)+(1|Area)+(1|Trans), data=dataset,
Thank you very much.
I'll try.
Jamie
Kevin Middleton-3 wrote:
Jamie,
Try it as:
data2 - read.table(pipe('pbpaste'), header = T)
I have a read.clipboard() function defined in ~/.Rprofile:
read.clipboard - function(){read.table(pipe('pbpaste'), header = T)}
Also see the
I'm looking for library which let mi convert dates
for example like this:
00-06-05 00:00
00-08-06 00:00
00-08-16 00:00
00-05-23 00:00
00-01-14 00:00
00-10-28 00:00
and as a result I want to get a 3 levels
--
View this message in context:
Thank you very much!
See, I'm a new user both for R and mac.. :-)
Jamie.
Prof Brian Ripley wrote:
On Tue, 28 Apr 2009, Jamie.lannister wrote:
Hi I'm a mac user. I have problems loading data from mac excel in R.
I'm using these script:
library(utils)
data2 -
Hi R community,
I am trying to obtain a sample from a matrix but sample(my.matrix) doesn't
do what I need. I have a matrix of 1287 interactions between the species in
columns and the species in rows and I want to obtain a smaller matrix with
say, 800 interactions, that may or may not have the
Hi, my name is Marcel R. Lopes. My problem is,
I made a code to calculate the estimates of a Cox model with random effects.
Used to optimize the R command for this. The estimates were calculated
correctly, but the Hessian matrix does not have good values. The same thing
was done in SAS and gave
Hi,
I need your help!!
I imported a big coloumn vector from a txt file but it results as mode
:list I want to change it in numeric otherwise I can't do my analysis.
This is what i get:
mode (data) - double
Error in eval(expr, envir, enclos) :
(list) object cannot be coerced to type 'double'
I am fitting logistic regression models, by defining my own link
function, and would like to get AICc values. Using the glm command
gives a value for AIC, but I haven't been able to get R to convert
that to AICc. Is there a code that has already been written for
this? Right now I am
Dear R-Help,
I have a data frame and a vector
df=data.frame(Letter=c(Z,Q,R,A,E,F), Number=c(11,32,4,1,9,3))
v=c(C,Q,R,A,E)
From df, I'd like to construct a subset of the field Number, with
deletions dictated by the vector, v, of letters. I've succeeded in doing
this for a single deletion (for
Hi all,
I'm reposting this with a more appropriate subject.
Do I need to define limits as the error message seems to suggest? If so, how?
The error message, my code, the output and the first few lines of my data are
all below.
Thank you!
Error in Getlim(at, allow.null = TRUE, need.all =
look at zoo and chron
On Wed, Apr 29, 2009 at 8:54 AM, Grześ gregori...@gmail.com wrote:
I'm looking for library which let mi convert dates
for example like this:
00-06-05 00:00
00-08-06 00:00
00-08-16 00:00
00-05-23 00:00
00-01-14 00:00
00-10-28 00:00
and as a result I want to get a
Hi useRs,
I had posted about Adobe FLEX talking to R for rich visualisation.
Reply from Jeffery Horner contained links to the
revolution-computing.com webpage which had information pertaining to
the Bay Users R group Meetup on Web Dashboards with R.
I have a very specific project that I need to
HI:
Is it possible to use the RODBC package within MS Access. I have been using
from R but was just wondering if it could be used along with R(D)COM. Something
like this:
Dim dbs As DAO.Database
Dim rst As DAO.Recordset
Dim myApp As StatConnector
Set myApp = New StatConnector
Hi James
its just simple calculus, since with
q90-qnorm(.9,me,sd)
q10-qnorm(.1,me,sd)
mean-(q90+q10)/2# the normal distribution is symmetric around the mean
sd-(q90-q10)/ (qnorm(.9)-qnorm(.1)) # between 10th and 90th are
qnorm(.9)-qnorm(.1)=2.563103sds
hth.
g...@ucalgary.ca schrieb:
If
Dear R-users,
I would like to know if is it possible to set a function argument as an
evaluated expression. I have tried several syntaxes, including the
following example, but could not get it anything to run. The plot
function is used here but I would like to later apply the same approach
Are you using optim() function in R?
Although I am not sure what your problem really is (as you haven't provided
sufficient information), it may be that the hessian computation in optim()
is not accurate enough for you.
Try the function hessian() in the numDeriv package. It is very accurate.
Dear R experts,
I am trying to estimate an ARMA 2,2 model with garch errors.
I used the following code on R 2.9.
#library
library(fGarch)
#data
data1-ts(read.table(C:/Users/falcon/Desktop/Time
Series/exports/goods1.csv), start=c(1992,1), frequency=12)
head(data1)
#garch
garchFit(formula.mean=
Dear List,
I have a wrapper function that draws a graph that I'd like to use in a
vector-like manner. The for-loop version I currently use is below.
library(ggplot2)
data(economics)
h - 600
w - 800
#--
draw_metric_by_date - function(
I've been using sna to work with some networks, and am trying to visualize
them easily. My networks are hierarchical (food webs). All of the layout
engines I've tried with gplot don't seem to plot hierarchical networks, as
one would using dot from graphviz. While I could do all of this by
Hi,
Thank you for littler webpage. It is quite interesting.
Tena Sakai
tsa...@gallo.ucsf.edu
-Original Message-
From: Wacek Kusnierczyk [mailto:waclaw.marcin.kusnierc...@idi.ntnu.no]
Sent: Wed 4/29/2009 1:22 AM
To: Tena Sakai
Cc: r-h...@stat.math.ethz.ch
Subject: Re: [R] Newbie R
1 - 100 of 130 matches
Mail list logo