Is there a way in GAMM to increase the amount of smoothing (something like
min.sp in GAM)?
This is a large-data problem with binomial response. I can get something
like the desired effect by
using very few knots, but this seems like a kludge.
thanks!
Scott Olsson
[[alternative HTML
Dear users,
I'm trying to produce a 3d bar plot but the x and y dimensions have categorical
data -- so I only want 3 points on each axis. So I try:
require(scatterplot3d)
mymat-data.frame(
x=c(1,1,1,2,2,2,3,3,3),
y=c(1,2,3,1,2,3,1,2,3),
z=c(1,2,3,4,5,6,7,8,9))
Hi, The following code, from Angelo Canty article on line
Resampling Methods in R: the boot Package, 2002, works fine for Angelo Canty
using R 2.6.0 on Windows XP.
It also works for me using R 1.2.1 and S-PLUS 2000 on Windows XP after
installing the S-PLUS
Hi, The following code, from Angelo Canty article on line Resampling
Methods in R: the boot Package, 2002, works fine for Angelo Canty using R
2.6.0 on Windows XP.
It also works for me using R 1.2.1 and S-PLUS 2000 on Windows XP after
installing the S-PLUS bootstrap library,
Dear R-helpers,
I would like to have my students read into R an online Stata dataset:
'http://www.stat.ucla.edu/projects/datasets/risk_project.dta'
I was able to read it into R after downloading it and converting it
with StatTransfer (http://www.stattransfer.com/).
Here is what happens when
Hi,
The following code, from Angelo Canty article on line Resampling Methods in
R: the boot Package, 2002, works fine for Angelo Canty using R 2.6.0 on
Windows XP.
It also works for me using R 1.2.1 and S-PLUS 2000 on Windows XP after
installing the S-PLUS bootstrap library,
Dear Alyaa,
you've now sent this message in many different forms to the
R-help mailing list.
The fact that you did not get a response may very well be
related to the way you asked your question,
but also that it probably is not an easy question to answer.
But sending your e-mail repeatedly to
Hi
I am using ggplot2 at the moment and I must say it is definitely better
then ggplot - good work.
My problem is that I am using facet_grid() in the following way:
p - ggplot(ssq, aes(x=year, y=-log(ssq)))
p + geom_point() + facet_grid(me*gi~cs*rz)
and it works nicely, except that I
Hallo,
Ihre E-Mailadresse [EMAIL PROTECTED] wurde in den Newsletterverteiler
hinzugefügt.
URL: http://www.fueralles.de/cgi-bin/nsl_pro/newsletter.cgi?id=steinfeuerfred
Mit freundlichen Grüssen Ihr Newsletter-Team
__
R-help@r-project.org mailing list
Hello
Is there an analogous function to ”rowsum” that count the numbers
according to a given vector (preferably larger than a given value)
instead of summing them?
E.g. rowsum(x, group)
X is the dataframe, A B C
5 0 0
Colleagues,
I am preparing scripts that will be used by others on both Windows and
Linux/OSX platforms. The scripts call an existing Fortran
application. The user may have any of a a variety of Fortran
compilers - my goal is to determine whether or not the test command
returns no input
Is this what you are looking for in doing computations on the columns
of a dataframe?
x -
data.frame(group=sample(1:2,10,TRUE),A=sample(1:4,10,TRUE),B=sample(1:5,10,TRUE),
+ C=sample(1:3,10,TRUE))
x
group A B C
1 1 2 4 1
2 2 4 1 2
3 1 3 2 1
4 2 3 2 2
5 2 4 2 2
On 26/01/2008 9:17 AM, Dennis Fisher wrote:
Colleagues,
I am preparing scripts that will be used by others on both Windows and
Linux/OSX platforms. The scripts call an existing Fortran
application. The user may have any of a a variety of Fortran
compilers - my goal is to determine
Dear R-help,
In nlme package, anova () can be used to test the difference between
two coefficients as shown on page 225
of Mixed Effects Models in S and S-Plus:
anova(fm2BW.lme, L = c(Time:Diet2 = 1, Time:Diet3 = -1))
Now my question is instead of test the difference between two
coefficients,
Greetings, R-ians:
I would like to know which version or R was used to create a given package.
I think I remember seeing that topic discussed recently but cannot find it
among my notes. Can anyone tell me how to determine which version of R
created a package?
Thanks.
Charles Annis,
library(help=lattice)
[...]
Built: R 2.6.0; i486-pc-linux-gnu; 2008-01-23 13:52:49; unix
[...]
G.
On Sat, Jan 26, 2008 at 11:44:53AM -0500, Charles Annis, P.E. wrote:
Greetings, R-ians:
I would like to know which version or R was used to create a given package.
I think I
These also works should you need to capture it:
packageDescription(lattice)$Built
On Jan 26, 2008 12:04 PM, Gabor Csardi [EMAIL PROTECTED] wrote:
library(help=lattice)
[...]
Built: R 2.6.0; i486-pc-linux-gnu; 2008-01-23 13:52:49; unix
[...]
G.
On Sat, Jan 26, 2008 at 11:44:53AM
On 26/01/2008 11:44 AM, Charles Annis, P.E. wrote:
Greetings, R-ians:
I would like to know which version or R was used to create a given package.
I think I remember seeing that topic discussed recently but cannot find it
among my notes. Can anyone tell me how to determine which version
hits=-2.6 tests=BAYES_00
X-USF-Spam-Flag: NO
Dear Geoff,
scatterplot3d was never intended to draw a plot for categorical data,
but you can tweak it as follows:
- Note that the help page says (stolen from ?par) The values of x and
y give the (approximate) number of tickmarks on the x and y
Many thanks to all. These methods all provide the needed info:
library(help= lattice)
packageDescription(lattice , fields=Built)
packageDescription(lattice)$Built
The first method produces nore complete information, should that also be
interesting.
Charles Annis, P.E.
[EMAIL PROTECTED]
Hi,
Not going into the reason why would you like to do such a thing anyway, I
would suggest:
1. Creating an R file that would contain all the functionality coded in
functions.
2. The only executed code would be a function that calls the menu() and
executes the appropriate functions.
Sources
Try like this:
switch(menu(c(Normal, Poisson)), rnorm(5), rpois(5, 3))
On 25/01/2008, WCD [EMAIL PROTECTED] wrote:
Hello, I would like to make my R script more interactive. Well, I want the
script to work like this:
I run the script from R Console and it offers me some options (press 1 for
Try this also:
noquote(tapply(x$Abundance, list(x$Plot, x$Species), paste))
On 25/01/2008, [EMAIL PROTECTED] [EMAIL PROTECTED] wrote:
Hello.
Does anyone know of an existing function that takes a list in the form of:
Plot1 Species1Abundance1
Plot1 Species2Abundance2
Plot2
On 26/01/2008 1:01 PM, Cleber Nogueira Borges wrote:
hello all,
when I start up the R and I execute o follow code:
ls()
character(0)
x=123
assign(test_x, x, envir = .GlobalEnv )
ls()
[1] test_x x
setwd('C:\\R\\etc')
save.image('TEST.RData')
q('no')
I have two
dear all,
is there a tutorial specifically teaching how to use R under unix? it
seems most of them are about using r under window.
thank you so much!
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the
Hi, I have a scalar valued function with several variables. One of the
variables is restricted to be non-negative. For example,
f(x,y)=sqrt(x)*exp(y), then x should be non-negative. I need the
gradient and hessian for some vector (0,y), i.e., I need the gradient and
hessian at the boudary of
When you double click it does it start up the same version of R
as when you manually start R?
If you manually start R, load it, save it under another name and
then exit R and double click the new name does that work?
On Jan 26, 2008 1:01 PM, Cleber Nogueira Borges [EMAIL PROTECTED] wrote:
hello
Try this:
tmpfile - tempfile()
download.file('http://www.stat.ucla.edu/projects/datasets/risk_project.dta',
tmpfile)
risk2 - read.dta(tmpfile)
unlink(tmpfile)
On 26/01/2008, Michael Kubovy [EMAIL PROTECTED] wrote:
Dear R-helpers,
I would like to have my students read into R an online Stata
Yes, the switch function works fine.
To store generated values in some variable I just use:
switch(menu(c(Normal, Uniform)), a-rnorm(5), a-runif(5))
But why does not work this? (I know there is waste of variables, but in
principle):
b-numeric(10)
switch(menu(c(Normal, Uniform)), a-1, a-2)
Hello,
On 1/26/08, Wensui Liu [EMAIL PROTECTED] wrote:
dear all,
is there a tutorial specifically teaching how to use R under unix? it
seems most of them are about using r under window.
thank you so much!
What exactly are you interested in: installing or using R? If it is
the latter, R is
Hi netters,
Suppose I have two data frames X and Y. X has three colnames A, B and C. Y has
three colnames A,B and D.
I want to combine them into one matrix, joining the rows having the same A and
B values (X$A==Y$A and X$B = Y$B). So the resulting dataframe has four
variables/columns: A,B,C
Henrique Dallazuanna wrote:
Try this:
tmpfile - tempfile()
download.file('http://www.stat.ucla.edu/projects/datasets/risk_project.dta',
tmpfile)
risk2 - read.dta(tmpfile)
unlink(tmpfile)
The stata.get function in the Hmisc package will automatically fetch
from http://... and provides
dear liviu,
first of all, i appreciate your insight.
i've been used r on windows for years and might use it on unix. i
understand r is similar across OSs. However, per my understanding,
still there is nontrivial difference between running on unix and
running on windows.
have a nice weekend.
On
Try this:
merge(x, y, all=T)
On 26/01/2008, zhihuali [EMAIL PROTECTED] wrote:
Hi netters,
Suppose I have two data frames X and Y. X has three colnames A, B and C. Y
has three colnames A,B and D.
I want to combine them into one matrix, joining the rows having the same A
and B values
Hi, Huali,
there is another way to do so
library(sqldf)
wanted - sqldf(select * from X inner join Y on X.A = Y.A and X.B = Y.A)
2008/1/26 zhihuali [EMAIL PROTECTED]:
Hi netters,
Suppose I have two data frames X and Y. X has three colnames A, B and C. Y
has three colnames A,B and D.
I want
Dear R-forumites,
I want to apply a Wilcoxon test on subsets of the data frame mydata,
splitted using the myindice variable. When I send :
wilcoxtest - by(mydata, mydata$myindice, function(x)
{wilcox_test(x$value~x$fact)})
I get :
Error in eval(expr, envir, enclos) : object x not
Henrique Dallazuanna [EMAIL PROTECTED] wrote in
news:[EMAIL PROTECTED]:
Try this:
merge(x, y, all=T)
That will give NA's in the non-matching rows, which may be what the
OP wanted. If, on the otherhand, only the rows with matches were
desired, the usage might be:
A.B - merge(X, Y, by =
On 26/01/2008, David Winsemius [EMAIL PROTECTED] wrote:
Henrique Dallazuanna [EMAIL PROTECTED] wrote in
news:[EMAIL PROTECTED]:
Try this:
merge(x, y, all=T)
That will give NA's in the non-matching rows, which may be what the
OP wanted. If, on the otherhand, only the rows with matches
hello Gabor,
the result in combobox is:
R version 2.6.1 (2007-11-26) - C:\R
I not have Perl installed! I remember that I deleted it!
Just double click on Rversions.hta from Windows Explorer.
Alternately place it anywhere in your path and then enter
its name in the Windows console without
OK. The system only knows about one R version.
What about the answer to my original second question.
If you start up R and then load the .Rdata file (which you
previously indicated works) and then
save your workspace using a different name from the original
Rdata file and then double click on
RSQLite 0.6-7 has been uploaded to CRAN and should hit a mirror near
you in the next few days.
This version changes the behavior of the dbGetQuery method to make it
more consistent with dbSendQuery. Specifically:
1. dbGetQuery now closes a complete result set as dbSendQuery does.
2. If there
Lately R has been behaving strange on my Linux (Ubuntu 7.10) machine,
with occasional segfaults. Today something else and reproducible
happened:
If I type the code below (meant for calibrating data), I get the error
message that the C stack usage is too close to the limit.
calcurve -
hi,
I tried this approach, but don't work too!
(Duncan Murdoch make this suggestion)
Cleber
OK. The system only knows about one R version.
What about the answer to my original second question.
If you start up R and then load the .Rdata file (which you
previously indicated works) and then
Have a look at rbind.fill() in the reshape library.
This will return a data.frame not a matrix but you can
always convert it to a matrix.
--- zhihuali [EMAIL PROTECTED] wrote:
Hi netters,
Suppose I have two data frames X and Y. X has three
colnames A, B and C. Y has three colnames A,B and
Gabor Grothendieck wrote:
Please clarify. What happens?
ok!
I start up the R, load my file by win menu: Load Workspace,
ls() show-me variables!
Then, I save my workspace with other name (test2.rdata, for example) and
close my E session.
When I start up R by double-click in my new RData file
Please clarify. What happens?
On Jan 26, 2008 5:02 PM, Cleber Nogueira Borges [EMAIL PROTECTED] wrote:
hi,
I tried this approach, but don't work too!
(Duncan Murdoch make this suggestion)
Cleber
OK. The system only knows about one R version.
What about the answer to my original
Dupouey Jean-Luc [EMAIL PROTECTED] wrote in
news:[EMAIL PROTECTED]:
On my system wilcox.test() exists, whereas wilcox_test() did not.
Help.search(wilcox_test) does not produce any clues. If you are using
a test from another package, you should specify which package the test
comes from. Web
Your function has 4 parameters and you are only calling with two. The
first statement:
caldist - function(cage=Cage, error=Error, sdev=Sdev, times=Times, By=By)
{
theta - seq(min(calcurve[,1]), max(calcurve[,1]), by=By);
use By which is not defined.
On Jan 26, 2008 4:29 PM, Maarten Blaauw
You could try
Rgui --vanilla
so that it does not run the scripts you have set up.
On Jan 26, 2008 5:31 PM, Cleber Nogueira Borges [EMAIL PROTECTED] wrote:
Gabor Grothendieck wrote:
Please clarify. What happens?
ok!
I start up the R, load my file by win menu: Load Workspace,
ls() show-me
hi Gabor,
I try it:
in *MS-DOS*:
C:\R\etcRgui --vanilla
In *R*:
a=1
b=2
c=3
save.image(c:\\R\\etc\\test.rdata)
q('no')
then I double-clicked' in 'test.rdata' (in c:\R\etc folder)
The R start without my variables: a,b,c
If I make again: R -- vanilla, and then load file, the variables are there!
Does anyone know any more than is in the following press release
about REvolution Computing and their commercialization of R?
http://www.intel.com/capital/news/releases/080122.htm
Intel Capital, the global investment arm of Intel Corporation, today
announced that it has invested in the Series
Try entering this at the Windows console:
assoc .Rdata
assuming the answer is:
.Rdata=RWorkspace
try this:
ftype RWorkspace
That will show you what command is executed when you double
click on an Rdata file. Perhaps you need to change that.
On Jan 26, 2008 6:40 PM, Cleber Nogueira Borges
Hi,
My R just froze. I can't get it to do anything. It gives Error: band value
message to everything I type. Does anyone know if R has a safe mode that I
could check for errors and perform diagnostics? I am using R 2.6.1 on
Windows XP.
ls()
Error: bad value
search()
Error: bad value
?ls
Error:
Sorry, indeed I forgot to put some of the factors in the code. Here it
is again, now updated:
calcurve - cbind(1:2e4, 1:2e4, rep(100, length=2e4));
caldist - function(cage, error, sdev=2, times=5, By=1)
{
calcurve - calcurve[which((calcurve[,2]+calcurve[,3]) =
cage-(times*error)),];
huali,
if i were you, i will create a view on the MySql server to aggregate
the data first and then use R to pull the data through this created
view. This is not only applicable to R but also a general guideline in
similar situation.
Per my understanding and experience, R is able to do data
How does the it compare if you read it into R and then do your
aggregate with sqldf:
library(sqldf)
# example using builtin data set CO2
CO2agg - sqldf(select Plant, Type, Treatment, avg(conc) from CO2
group by Plant, Type, Treatment)
# or using your data:
Xagg - sqldf(select Group, Age, Type,
When and how does this happen? In Rgui or Rterm? How do you launch
R? Have you tried the obvious and restarted R? If you share the
startup messages, sessionInfo() and everything else you think could be
helpful, someone might be able to help you. There is a vanilla
options to start R, e.g.
hi Gabor,
Now, seems an light, because...
C:\R\etcassoc .Rdata
.Rdata=RWorkspace
C:\R\etcftype RWorkspace
Tipo de arquivo 'RWorkspace' não encontrado ou nenhum comando open
associado a ele.
*File type 'RWorkspace' not found or none command open associated for it*.
seems something with error
On 26/01/2008 5:31 PM, Cleber Nogueira Borges wrote:
Gabor Grothendieck wrote:
Please clarify. What happens?
ok!
I start up the R, load my file by win menu: Load Workspace,
ls() show-me variables!
Then, I save my workspace with other name (test2.rdata, for example) and
close my E
http://www.revolution-computing.com/revolution_company.html
seems to lead to something albeit vague.
Lots of apparently big names but no indication of what
added value the company might bring to an R user.
--- John Maindonald [EMAIL PROTECTED]
wrote:
Does anyone know any more than is in the
On 26/01/2008, WCD [EMAIL PROTECTED] wrote:
Yes, the switch function works fine.
To store generated values in some variable I just use:
switch(menu(c(Normal, Uniform)), a-rnorm(5), a-runif(5))
But why does not work this? (I know there is waste of variables, but in
principle):
How large is your dataframe? How much memory do you have on your
system? Are you paging? Here is a test I ran with a data frame with
1,000,000 entries and it seems to be fast:
n - 100
x - data.frame(A=sample(LETTERS,n,TRUE), B=sample(letters[1:4],n,TRUE),
+ C=sample(LETTERS[1:4], n,
Sudhindra,
You should read the Posting Guide before asking for help:
http://www.R-project.org/posting-guide.html
It has greats hints on how to solve problems like this yourself, and you
should try them before posting.
It might suggest that this is a good thing to run in R if you need
Sudhindra Gadagkar [EMAIL PROTECTED] wrote in
news:[EMAIL PROTECTED]:
Hi, I can't believe I am unable to do this. I just installed R
(never used it before, but want to learn it). At present, though.
all I want to do is export a data file (specifically crabs in the
MASS library) to a text
first of all, which is your code to export the data in R? did you use
write.table?
secondly, it is totally a waste of time to export data from R for use
in spss. for whatever can be done by spss, r can handle better.
On Jan 26, 2008 9:38 PM, Sudhindra Gadagkar [EMAIL PROTECTED] wrote:
Hi, I
Hi,
I’m using the gls function to fit a model with autoregressive
residuals. Does anyone know how to do weighted least squares with
autoregressive residuals in gls?
Thanks,
Sincerely,
Xiaopeng
__
R-help@r-project.org mailing list
Hi, all.
I need help on improving the efficiency of my R
simulations. Below is a function that simulates a
change point model. It first generates a sequence of
three dimensional ARMA(1,1) observations, then
calculates the one step ahead prediction errors, some
statistic is calcualted and
Dear List,
I am not sure how should i optimize a log-likelihood numerically:
Here is a Text book example from Statistical Inference by George Casella, 2nd
Edition Casella and Berger, Roger L. Berger (2002, pp. 355, ex. 7.4 # 7.2.b):
data = x = c(20.0, 23.9, 20.9, 23.8, 25.0, 24.0, 21.7, 23.8,
You asked for a hint.
library(MASS)
x - c(20.0, 23.9, 20.9, 23.8, 25.0, 24.0,
21.7, 23.8, 22.8, 23.1, 23.1, 23.5,
23.0, 23.0)
fitdistr(x, gamma)
shape rate
316.56387213.780766
(119.585534) ( 5.209952)
To do it with more general and
69 matches
Mail list logo