Hi
Mikkel Grum wrote:
> Dear useRs,
>
> Why are the rotated blue and yellow boxes in the example below clipped
> outside of 6 x 6 inch window in the middle of the page?? Where does the 6 x 6
> inch window come from? I would like to make use of the entire page.
6x6 corresponds to the default
Hi, I am new to R. I don't have strong background of statistics. I am
a student of Geotechnical Engineering. I tried to run a nonlinear
regression for a three-variable function, that is
N = f(CSR, ev) # N is a function of CSR and ev, and N = CSR/(A
+B*CSR), wherer (A,B) are function of ev.
Thank you Charilaos, I think I misunderstood the 'selfStart'. After
checking 'selfStart', I think I don't have the ability to create a
selfStart object. I don't have good sense to guess the initial
values. So I tried 'poor guess' by nls. But I got erroe message as:
Error in c * (ev^2) : non-
Dear R-project
I'm sure this is a very simple problem, but I've not been able to find a good
example to guide me (and I'm very new to R).
I have been looking at time constants of inactivation (rate of turn-off) of
sodium channel currents in neurons at 5 different potentials, looking to see if
Hello, R users.
I need help. When I run one of my own functions, I got follwoing error
message.
Error: cannot allocate vector of size 350493 Kb
So I check the memory size.
memory.limit()
[1] 1610612736
Is it enough for vector of size 350493 Kb? Am I missing something?
Thank you for your
Hi Mike,
You function is discontinuous at -0.8, so you can
expect everything :-}! But this is not the only
problem. The algorithm for optimize never gets there.
In general there exists no universal method to find
the global maximum of a function (unless it satisfies
certain conditions). You can al
This would probably also be interesting to some:
On 2007-October-01 , at 00:48 , hadley wickham wrote:
>> That's great!
>> In fact I think I found exactly what I was looking for. I can just
>> do:
>> p = ggplot() + coord_equal()
>> p$aspect.ratio = 1
>> to set up the plot, and t
This was meant to be sent on the list:
On 2007-September-30 , at 23:12 , jiho wrote:
> On 2007-September-30 , at 21:01 , hadley wickham wrote:
[...]
>>> As expected there is nothing in the data part of the p object
p$data
>>> NULL
>>>
>>> But there is no data specification either in th
Hi all,
Consider the following function:
my.func = function(x){
y=ifelse(x>-.5,0,ifelse(x< -.8,abs(x)/2,abs(x)))
print(c(x,y)) #print what was tested and what the result is
return(y)
}
curve(my.func,from=-1,1)
When I attempt to find the maximum of this function,
Hi again
Here is the sample:
> yuk
Jan Feb Mar Apr May Jun Jul Aug Sep Oct Nov Dec
2005 41.26 40.02 38.24 35.37 39.35 38.90 43.51 40.32 38.14 41.04 41.78 40.48
2006 40.55 42.15 42.30 39.93 38.12 35.79 34.71 34.29 36.27 37.33 37.97 40.16
2007 40.74 39.59 36.74 37.87 38.8
This is close to what you want. I created the list by hand, but you
can create it in your processing loop. Once you have the list
created, you can create your own print routine.
> x <- runif(100)
> z <- runif(100)
> y <- runif(100)
>
>
> # I am doing this by hand, but you could easily automate i
Try this:
x <- ts(11:30, start = 2000, freq = 12)
data.frame(x = c(x), time = c(time(x)))
The time will be in years plus fraction of a year.
On 10/1/07, Edna Bell <[EMAIL PROTECTED]> wrote:
> Dear R gurus
>
> I would like to take a monthly time series and convert it to a data
> frame without los
Hi Edna,
Can you send a small subset of the data as an example and the
function call you used to read the data in originally? It might be
helpful in understanding why you're losing the "time element".
Jeff.
On Oct 1, 2007, at 12:27 AM, Edna Bell wrote:
> Dear R gurus
>
> I would like to ta
Dear R gurus
I would like to take a monthly time series and convert it to a data
frame without losing the tsp items, pleae
I've tried as.data.frame and data.frame but I get the series without
the time element.
Any suggestions, please?
tia
Edna Bell
Jim,
You are indeed trying to help, again my thanks.
What I want to do is make a single structure - a table is an apt description
that will summarize all the regressions, something like:
Estimate Std. Error t value Pr(>|t|)
(Intercept) lm(formula = y ~ x) 4
Thank you Moshe,
I understand you point, but I would hope that I could use summary to save my
self some work. I need to do what I described in my original Email to the list
server on tens of regressions.
John
John
John Sorkin M.D., Ph.D.
Chief, Biostatistics and Informatics
University of Maryl
Jim,
Again thank you for your quick reply. Your suggestion does not give me exactly
what I want:
> whatIwant<-list(,summary(fitdelete)$call,summary(fitdelete)$coefficients)
> whatIwant
[[1]]
lm(formula = y ~ x)
[[2]]
Estimate Std. Error t value Pr(>|t|)
(Intercept) 4.927791 2.6
try using a 'list':
whatIwant<-list(call=summary(myreg)$call, coef=summary(myreg)$coefficients)
On 9/30/07, John Sorkin <[EMAIL PROTECTED]> wrote:
> Widows XP
> R 2.3.1
>
> I have been trying to make a data structure that will contain both the
> coefficients from a linear regression along with c
There are many. Check out the 'rgl' package to see an interactive
version of 3D plots.
On 9/30/07, John Sorkin <[EMAIL PROTECTED]> wrote:
> Windows XP
> R 2.3.1
>
> I have a funciton
> fit1<-lm(y~x+z)
> Is there a function that will produce a 3-dimensional plot of y,x,z?
>
> I looked at the help
Widows XP
R 2.3.1
I have been trying to make a data structure that will contain both the
coefficients from a linear regression along with column and row titles AND the
call, i.e.
myreg<-lm(y~x+y+z)
whatIwant<-cbind(c(summary(myreg)$call,"",""),summary(myreg)$coefficients)
Neither the statement
Hi Irene,
In this case the computer is right - your gradient is
really singular!
If you scale all you parameters (g,a and b) by any
nonzero constant c nothing changes, meaning that there
is a "degree of freedom" and this causes the gradient
to be singular.
You can check whether g = 0 (and then y =
Windows XP
R 2.3.1
I have a funciton
fit1<-lm(y~x+z)
Is there a function that will produce a 3-dimensional plot of y,x,z?
I looked at the help files, but did not find a clean answer to my question.
Thanks,
John
John Sorkin M.D., Ph.D.
Chief, Biostatistics and Informatics
University of Maryland
Luis Guillermo DÃaz Monroy
Profesor Asociado Departamento de EstadÃstica
Universidad Nacional de Colombia
--- Begin Message ---
Dear members of R community,
I want to know about the R procedures (or routines) to carry out the
selection of covariates in an Aalen or
cox mo
The model is not identifiable since if (a, b, g) is a solution then
so is every multiple of it.
On 9/30/07, Irene Mantzouni <[EMAIL PROTECTED]> wrote:
> Dear all,
>
> I would like to fit a non-linear model of the form:
> y=g*x/(a+b*x)
> with nls().
> However this model is somehow overparameterized
Dear all,
I would like to fit a non-linear model of the form:
y=g*x/(a+b*x)
with nls().
However this model is somehow overparameterized and I get the error message
about
singular gradient matrix at initial parameter estimates.
What I am interested in is to make inference about parameters b and g
On 9/30/07, hadley wickham <[EMAIL PROTECTED]> wrote:
> On 9/30/07, jiho <[EMAIL PROTECTED]> wrote:
> > On 2007-September-30 , at 22:40 , hadley wickham wrote:
> > >> hadley wickham wrote:
> > >>> [...]
> > >> PS if one specifies "errorbars" without specifying min and max one
> > >> gets
> > >> th
On 9/30/07, jiho <[EMAIL PROTECTED]> wrote:
> On 2007-September-30 , at 22:40 , hadley wickham wrote:
> >> hadley wickham wrote:
> >>> [...]
> >> PS if one specifies "errorbars" without specifying min and max one
> >> gets
> >> the error
> >>
> >> Error in rbind(max, max, max, min, min, min) :
> >
> Thanks for trying this, Hadley, because the comparison
> is instructive in terms of the difference between the
> communication goals of analysis and presentation graphs.
Yes, and I think it's a difference that not enough people are familiar with.
> Actually, one should regard income as the inde
On 2007-September-30 , at 22:40 , hadley wickham wrote:
>> hadley wickham wrote:
>>> [...]
>> PS if one specifies "errorbars" without specifying min and max one
>> gets
>> the error
>>
>> Error in rbind(max, max, max, min, min, min) :
>> cannot coerce type closure to list vector
>>
>>
> hadley wickham wrote:
> > I'd do this a little differently, using the reshape
> > (http://had.co.nz/reshape) and ggplot2 (http://had.co.nz/ggplot2)
> > packages:
> >
> > library(reshape)
> > library(ggplot2)
> >
> > # Get data in format required for reshape
> > df <- rename(df, c("OD" = "value"))
hadley wickham wrote:
> I'd do this a little differently, using the reshape
> (http://had.co.nz/reshape) and ggplot2 (http://had.co.nz/ggplot2)
> packages:
>
> library(reshape)
> library(ggplot2)
>
> # Get data in format required for reshape
> df <- rename(df, c("OD" = "value"))
>
> # Summarise and
On Sun, 2007-09-30 at 20:56 +0200, Tom Cohen wrote:
> Dear list,
> I have following table
>
> ee<-table(ID,Day)
> ee
> Day
> ID2 3 4 5 6 7 9 10 14 16
> 35 5 0 0 3 1 0 0 5 0 0
> 36 0 0 0 0 0 0 0 1 0 0
> 43 13 15 15 0 0 13 13 15 13 15
> 46 0
On 9/30/07, Ben Bolker <[EMAIL PROTECTED]> wrote:
> Marcelo Laia gmail.com> writes:
>
> >
> > Hi,
> >
> > I have a data set like this:
> >
>
> [snip]
>
> > I need to plot a graph OD over the time for each one mutant with error bars.
> >
> >
>
>
> ## I put your data in a temporary file, this reads
> > There are a few ways you could describe the graph you want. Here's
> > the one that I'd probably choose:
> >
> > ggplot(mapping = aes(x = log, y = lat)) +
> > geom_path(data = coast) +
> > geom_point(data = coords) +
> > coord_equal()
> >
> > We don't define a default dataset in the ggplot cal
Dear list,
I have following table
ee<-table(ID,Day)
ee
Day
ID2 3 4 5 6 7 9 10 14 16
35 5 0 0 3 1 0 0 5 0 0
36 0 0 0 0 0 0 0 1 0 0
43 13 15 15 0 0 13 13 15 13 15
46 0 1 0 0 0 0 0 0 0 0
58 0 0 0 0 0 0 4 4 0 0
and want to
Hi,
I tried to load a .RData object on unix system using R, it gives error:
Error: restore file may be empty -- no data loaded
In addition: Warning message:
file 'junk3.RData' has magic number ''
Use of save versions prior to 2 is deprecated
This happens only for using MY user account for th
Marcelo Laia gmail.com> writes:
>
> Hi,
>
> I have a data set like this:
>
[snip]
> I need to plot a graph OD over the time for each one mutant with error bars.
>
>
## I put your data in a temporary file, this reads it
x = read.table("tempdata.txt",header=TRUE)
## compute means and sta
Omar Baqueiro wrote:
> Hello,
>
> I have tested a distribution for normality using the Shapiro-Welch
> statistic. The result of this is the following:
>
>
> Shapiro-Wilk normality test
>
> data: mydata
> W = 0.9989, p-value = 0.8791
>
>
> I know that the p-value > 0.05 (for my purpos
On 9/30/07, Prof Brian Ripley <[EMAIL PROTECTED]> wrote:
> >>> Consider the following piece of code:
> >>>
> >>> pdf(file="figure.pdf", family="Palatino")
> >>> plot(0,0,type='n', xlim=c(-20,20), ylim=c(0,2),xlab="",ylab="",axes=F)
> >>> text(-1.4,1.168,expression(italic("The font looks different w
On Sun, 30 Sep 2007, bbolker wrote:
> Michael Friendly wrote:
>>
>> In a long session, producing multiple graphs, I sometimes repeatedly
>> change par() settings, particularly with multi-row/col displays.
>> If I'm using a script, I'll do
>>
>> op <- par(newsettings)
>> ... plots ...
>> par(op)
>>
On Sun, 30 Sep 2007, Paul Smith wrote:
> On 9/30/07, Marc Schwartz <[EMAIL PROTECTED]> wrote:
>>> Consider the following piece of code:
>>>
>>> pdf(file="figure.pdf", family="Palatino")
>>> plot(0,0,type='n', xlim=c(-20,20), ylim=c(0,2),xlab="",ylab="",axes=F)
>>> text(-1.4,1.168,expression(italic
On 2007-September-30 , at 18:35 , hadley wickham wrote:
>> The ggplot book specifies that "[ggplot] makes it easy to combine
>> data from multiple sources". Since I use ggplot2 as much as I can
>> (thanks it's really really great!) I thought I would try producing
>> such a plot with ggplot2.
>>
>>
On 9/30/07, Marc Schwartz <[EMAIL PROTECTED]> wrote:
> > Consider the following piece of code:
> >
> > pdf(file="figure.pdf", family="Palatino")
> > plot(0,0,type='n', xlim=c(-20,20), ylim=c(0,2),xlab="",ylab="",axes=F)
> > text(-1.4,1.168,expression(italic("The font looks different when this
> > i
Dear Letticia,
Are you using R-help for your homework?
See your previous postings:
1. 15th September:
With a single R command complete the following:
create a vector calles seqvec that repeats the sequence 1, 3,6,
10,15,21.( I was trying to use c() but this does not work)
create a 5-row, 6-col
Hi JiHO,
> The ggplot book specifies that "[ggplot] makes it easy to combine
> data from multiple sources". Since I use ggplot2 as much as I can
> (thanks it's really really great!) I thought I would try producing
> such a plot with ggplot2.
>
> NB: If this is possible/easy with an other plotting
On Sun, 2007-09-30 at 16:49 +0100, Paul Smith wrote:
> Dear All,
>
> Consider the following piece of code:
>
> pdf(file="figure.pdf", family="Palatino")
> plot(0,0,type='n', xlim=c(-20,20), ylim=c(0,2),xlab="",ylab="",axes=F)
> text(-1.4,1.168,expression(italic("The font looks different when this
Dear All,
Consider the following piece of code:
pdf(file="figure.pdf", family="Palatino")
plot(0,0,type='n', xlim=c(-20,20), ylim=c(0,2),xlab="",ylab="",axes=F)
text(-1.4,1.168,expression(italic("The font looks different when this
is seen with Acrobat Reader!")),xpd=T)
dev.off()
When viewing the
Michael Friendly wrote:
>
> In a long session, producing multiple graphs, I sometimes repeatedly
> change par() settings, particularly with multi-row/col displays.
> If I'm using a script, I'll do
>
> op <- par(newsettings)
> ... plots ...
> par(op)
>
> but sometimes I do things on the fly an
Dear useRs,
Why are the rotated blue and yellow boxes in the example below clipped outside
of 6 x 6 inch window in the middle of the page?? Where does the 6 x 6 inch
window come from? I would like to make use of the entire page.
> library(grid)
> pdf(file = "FarmMaps.pdf", paper = "a4")
> pushV
Hi,
In my package RGtk2, there's a directory called 'src/RGtk2' that contains,
like all the other directories in the package, a '.svn' directory. It seems
that R CMD build is somehow missing that one '.svn' in 'src/RGtk2', even
though it excludes all the other instances of '.svn'.
I've tried putt
50 matches
Mail list logo
