;
> >> or
> >>
> >> a<-array(dim=c(1,3,10))
> >>
> >> error message:reached getOption("max.print") -- omitted 6667 row(s)
> and 6
> >> matrix slice(s) ]
>
> PS> That is not an error messa
NA NA
[7,] NA NA NA
[8,] NA NA NA
[9,] NA NA NA
[10,] NA NA NA
[11,] NA NA NA
>
See
?getOption
Paul
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting
Dear Ken.
You could also try and use RSPerl (http://www.omegahat.org/RSPerl/). It
allows one to use R commands in Perl and vice-versa.
regards,
Paul
Pierce, Ken schreef:
> Is there a way to run a simple perl script from R?
>
>
> Kenneth B. Pierce Jr.
>
> Research Ecolo
On 9/8/07, Gabor Csardi <[EMAIL PROTECTED]> wrote:
> Paul, what is the question? If the question is why you get this
> warning message, the reason is that the 'lib' argument is missing
> and install.packages is using '/usr/lib/R/library'.
>
> If you want
2007-06-27)
>
R is here installed on Fedora 7.
Thanks in advance,
Paul
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
0 0 -0.5
-0.50 0.1 -0.4
-0.50 0.3 -0.2
0.5 0 -0.50
-0.50 0.5 0
0.5 0 -0.40.1
0.5 0 -0.20.3
0.5 0 0 0.5
are all solutions for your problem.
Paul
_
vn/trunk/R/sqldf.R";)
> > sqldf(c("delete from BOD where demand = 8.3", "select * from BOD"))
> Time__1 demand
> 1 2 10.3
> 2 3 19.0
> 3 4 16.0
> 4 5 15.6
> 5 7 19.8
I see, Gabor, but I would expect as more natu
n there, but it
> should be simple... Ummm...
Thanks, Barry and Gabor. Please, look at the following:
> library(sqldf)
Loading required package: RSQLite
Loading required package: DBI
Loading required package: gsubfn
Loading required package: proto
> packageDescription("sqldf", fields="Depends")
[1] "R (>= 2.5.1), RSQLite (>= 0.5-5), gsubfn"
>
packageDescription does not mention the packages DBI and proto.
Paul
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
Dear All,
Is sqldf equipped with delete queries? I have tried delete queries but
with no success.
Thanks in advance,
Paul
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http
package:graphics" "package:grDevices" "package:utils"
[10] "package:datasets" "package:methods" "Autoloads"
[13] "package:base"
The packages
RSQLite
DBI
gsubfn
proto
were not detached.
Thanks in advance,
Paul
__
Interval Number_reads
> 0-100 685
> 100-200 1082
> 200-300 744
> 300-400 4213
What abo
5847847
>
> I would like to know how to compute the length of each SEQUENCE.
Maybe the following code?
> data
var1 var2
1 This is a string 12
2 This is another string 34
> nchar(data[,1])
[1] 16 22
>
Paul
_
large redesign problem that is not anywhere
near the top of anyone's todo list (to my knowledge).
As Michael mentioned, the gridSVG package allows you to draw using the
grid package and access some of the fancier SVG features (including
embedding scripts). It has its own problems of course,
and the ASCII file into the spatial classes
presented by sp (SpatialGrid for the CO2 data (I assume that it is a
grid) and SpatialPoints for the ASCII data). This could be done using
the rgdal-package.
good luck!
Paul
Yogesh Tiwari schreef:
> Hello 'R' Users,
>
> I have
() with "yes". The manual says that q() is
equivalent to save(list = ls(all=TRUE), file = ".RData"). There seems to be no
way to set ascii or compression of save through q function, unless the q
function is replaced explicitly with save(list = ls(all=TRUE), file = ".R
On 9/2/07, Paul Smith <[EMAIL PROTECTED]> wrote:
> > I have a function like this:
> >
> > fun <- function (x, y) {
> > a <- log(10)*y
> > b <- log(15)*x
> > extr <- a-b
> > extr
> >
tr
}
x <- c(1,2,3)
y <- c(4,5,6)
combs <- expand.grid(x,y)
for (i in 1:nrow(combs))
cat(fun(combs[i,1],combs[i,2]),"\n")
Paul
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read t
without mixing compilers on this machine, i.e., I don't want
to use gcc and Studio 11 together. Is there some way to get Studio 11 to
compile the standalone library by altering the makefile?
Thanks for any advice/explanation people can offer. I'm a newcomer to compiling
code on UNI
experience?
Paul.
- Original Message
From: Barry Rowlingson <[EMAIL PROTECTED]>
To: Eric Turkheimer <[EMAIL PROTECTED]>
Cc: r-help@stat.math.ethz.ch
Sent: Wednesday, August 22, 2007 9:43:57 AM
Subject: Re: [R] Synchronzing workspaces
Eric Turkheimer wrote:
> How do people go abou
Dear Chris,
I use Python (http://www.python.org) in combination with Rpy
(http://rpy.sourceforge.net/). Rpy enables you to use R commands inside
Python, not the other way around. Works quite well, also for different R
versions (I currently run R 2.5.1 under Linux).
cheers,
Paul
Chris Parkin
ode.
>
Dear Pratap,
You could try something like this:
x = c(1,0,1,1)
y = c(1,2,3,NA)
y_rev = rev(y)
ifelse(x == 1, y_rev, y)
hope this helps,
Paul
--
Drs. Paul Hiemstra
Department of Physical Geography
Faculty of Geosciences
University of Utrecht
Heidelberglaan
On 8/30/07, Vladimir Eremeev <[EMAIL PROTECTED]> wrote:
>
> Use rscript
>
> Rscript myscript.R
> or
> Rscript -e 'cat("Hello!\n")'
>
> will show Hello! on the console.
>
> R CMD BATCH writes its output to the file myscript.Rout
Thanks, Vl
; 'Hello!'.
Thanks, Barry. Indeed, the file myscript.Rout exists and contains the
output of cat. I was expecting a behavior similar to the bash scripts.
And by the way, cannot a R script write only on the console and just
what one tells it to write, likewise bash scripts?
Paul
__
Dear All,
I am trying to write my first R script. The code is simply
cat("Hello!\n")
However, when I run
$ R CMD BATCH myscript.R
I do not see "Hello!" on the console. I am using Fedora 7 (Linux) and R-2.5.1.
Any ideas?
Than
The instance of R running will be immediately killed and then you can
start R again.
Paul
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
anyone know a method to do this or do I have to write something. :(
If I do any starter code :D jj. If I've read something wrong or
misunderstood my apologies.
cheers,
Paul
--
Research Technician
Mass Spectrometry
o The
/
o Scripps
\
o Resear
Hi
Antje wrote:
> Hi Paul,
>
> > You are getting the error because you are setting the figure region to
> > be larger than the current device (typically 6 or 7 inches wide/high).
> > You SHOULD be getting the error when you try par(fin), BUT there is a
> > c
there IS a check on the values, the
values are larger than the current device and you get the error. Now,
because there is an error in resetting par(fig), that parameter is not
reset, so when you type par()$fin (or, equivalently, par("fin")) after
the heatmap() call, you get the last s
n the symbol font).
I have added a feature to plotmath in the development version (to be R
2.6.0) to make it easier to get at any symbol you want. In the next
version of R, you will be able to do something like this ...
grid.newpage()
grid.text(expression("whatever "*symbol("\042
> www.users.skynet.be/fa244930/fotos/outlined.jpg
You mean the diameter measured in number of pixels?
Paul
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.or
n
> calling f...)
> I've looked under formals, alist etc. but so far without luck.
I hope that the following helps:
> f <- function(x) {sum(x)}
> f(c(2,3))
[1] 5
> f(c(2,3,5))
[1] 10
>
Paul
__
R-help@stat.math.ethz.ch mailing
e code:
>
> \newpage
> \includegraphics[scale=1.0]{\Sexpr{print(location)}}
> \newpage
>
> (in tex this looks like):
> \newpage
> \includegraphics[scale=1.0]{/home/arjun/sample.pdf}
> \newpage
Try putting your image in a figure* environment (should go full width of
t
m within R. It can be
used for computer algebra, exact arithmetic, ASCII pretty printing and
R to TeX output.»
Paul
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-proje
Hi Wesley,
Try the text() function. An example:
a = rep(10,10)
b = seq(1,10)
plot(a,b)
text(a,b, labels = b, pos = 4, offset = 0.7)
?text
hth,
Paul
Wesley Roberts schreef:
> Hi All,
>
> I am busy using R to do some regression modelling and have been using
> plot(x,y,"&
ist("sp.points",
geocode)) # Plot the predictions (grid) and the points (geocode)
hope this helps,
Paul
Lawrence D. Brenninkmeyer schreef:
> Hi-
>
> I am trying to find a way to diffuse GIS data on a European map. I have a
> dataset consisting of particular locat
ple Behrens-Fisher-Test`
cmpeffect lower.cl upper.cl p.value.1s p.value.2s
1 1-2 0.643 0.4610459 0.8256208 0.08595894 0.14750647
2 1-3 0.444 0.2576352 0.6312537 0.99636221 0.75376639
3 2-3 0.328 0.1602449 0.4964218 1. 0.04476692
What elementary error am I ma
I don't think you can define a matrix this large in R, even if you have
the memory. Then, of course, inverting it there may be other programs
that have limitations.
Paul
Jiao Yang wrote:
> Can R invert a 16x16 matrix with all positive numbers? Than
dvice.
Paul
--
Dr. Paul Cleary,
Specialist Registrar in Public Health,
Chester Health Protection Unit.
[[alternative HTML version deleted]]
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do rea
determine where to draw each
piece of the puzzle.
(ii) In another way, this also sounds a bit like the gridSVG package
(http://www.stat.auckland.ac.nz/~paul/index.html), but instead of
converting grid objects to SVG, you convert them to pic. This does
avoid going through the R graphics engi
Hi
Emilio Gagliardi wrote:
> haha Paul,
>
>
> It's important not only to post code, but also to make sure that other
> people can run it (i.e., include real data or have the code generate
> data or use one of R's predefined data sets).
>
>
&
? Text in the corner!\nWhere do we want it?
Here!",
x=unit(2, "mm"), y=unit(1, "npc") - unit(2, "mm"),
just=c("left", "top"))
Paul
> Thanks
>
> Benilton Carvalho wrote:
>> maybe this is what you want?
&g
Hi
Emilio Gagliardi wrote:
> Hi Paul,
>
> I'm sorry for not posting code, I wasn't sure if it would be helpful without
> the data...should I post the code and a sample of the data? I will remember
> to do that next time!
It's important not only to post code, bu
.rect.345" and then
draw in just the sides you want, which would require getting to the
right viewport, for which you'll need to study the viewport tree (see
current.vpTree())
> size of the font for the axis labels independently of the tick marks. I'd
That's the one we
setMethod("read.metlin", "xcmsRaw", function(xml) {
#Parsing the METLIN XML File
reading<-readLines(xml)
#do rest of script
})
Any help as to why I'm getting the inherited method error would be great.
Cheers,
Paul
___
Dear All,
Is there some package to do optimal control?
Thanks in advance,
Paul
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and
ow all of the common work arounds and such. At the moment the
only method I can think of is to split the data using condition
statements and manually running the regression on each set. Thanks or
your help
-Paul
_
CONFIDENTIALITY NOTICE\ IMPORTANT: This communicatio
On 8/2/07, Paul Smith <[EMAIL PROTECTED]> wrote:
> > I'm using the function constrOptim together with the "SANN" method and
> > my objective function (f) has two parameters. One of the parameters
> > needs be into (2^(-10), 2^4) range and the other into (2
> range. How can I do it using constrOptim??
I think that you, André, do not need constrOptim; optim is enough. See
?optim
Paul
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide
Dear R-helpers,
I have several versions of R installed on my computer, and I cannot do
without any of these.
However RCDOM seems to authorize only one version installed. Do you know any
means to overcome this problem ?
Thank you very much for your response.
Paul Poncet
swap off. There are
debugging situations where this turns out to be useful.
HTH,
Paul
zhihua li wrote:
> Hi netters,
>
> I'm using the 64-bit R-2.5.0 on a x86-64 cpu, with an RAM of 2 GB. The
> operating system is SUSE 10.
> The system information is: -uname -a
> Linux someo
x1>=0;
> x1<=1;
> x2>=0;
> x2<=1;
> x3>=0;
> x3<=1;
>
> Any suggests ?
What do you mean by 'radq', Massimiliano?
Paul
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-hel
Hi Deepayan,
"Deepayan Sarkar" <[EMAIL PROTECTED]> schrieb am 06.07.2007
02:05:02:
> On 7/5/07, Paul Matthias Diderichsen
> <[EMAIL PROTECTED]> wrote:
>> library(lattice)
>>> xyplot(speed~dist|speed, data=cars, layout=c(3,3))
> If this is your use ca
Dear Prof. Ripley,
Prof Brian Ripley <[EMAIL PROTECTED]> schrieb am 05.07.2007 21:46:20:
> > Dear S-users.
> This is the help forum for R users
Indeed. (How embarrasing not to be able to spell a one-letter word
correctly...)
>> How do I change pages on an X11 graphics
>> device under linux
Dear S-users.
This should be an easy one: How do I change pages on an X11 graphics
device under linux?
I thought that the page-up/page-down keys were supposed to do the trick,
but the frame (window) seems to be kind of immune to any kind of keyboard
input. The only reaction I ever see is that t
On 7/4/07, Paul Smith <[EMAIL PROTECTED]> wrote:
> On 7/4/07, RAVI VARADHAN <[EMAIL PROTECTED]> wrote:
> > My point is that it might be better to try multiple (feasible) starting
> > values for constrOptim to ensure that you have a good local minimum, since
>
t constrOptim find boundary solutions
correctly? The set of solution that I got is below.
Paul
2.67682495728743e-080.676401684216637 5.18627076390355e-09
0.00206463986063195 0.871859686128364.32039325909089e-11
0.9996234
3.
> However you can get similar result with using points
>
Yes I new that, but I wanted to try and go without an if() for
deciding between the first and consecutive columns.
Thnx for helping out!
Paul Lemmens
__
R-help@stat.math.ethz.ch mailing l
orrespond to an impossible linear
system, but it works in many cases.
Paul
> - Original Message -
> From: Paul Smith <[EMAIL PROTECTED]>
> Date: Tuesday, July 3, 2007 7:32 pm
> Subject: Re: [R] Fine tunning rgenoud
> To: R-help
>
>
> > On 7/4/07, Ravi Vara
<- par(new=TRUE) either before or immediately
after the first plot statement (the type='n' one) in the above code
fragment, the resulting graph still only shows one column of data.
Have I misinterpreted the instructions or the functionality
Hi
Dieter Menne wrote:
> Dear UseRs called Hadley, or Paul,
>
> I am trying to print an edited ggplot2/grid graphics to a metafile. With the
> commented line below it works, but when I edit the plot by uncommenting the
> line, it fails, because it's illegal to have 2 graphi
ints are satisfied, but the solution
given by constrOptim is not interior; the gradient is not equal to
zero.
Paul
> -Original Message-
> From: [EMAIL PROTECTED]
> [mailto:[EMAIL PROTECTED] On Behalf Of Paul Smith
> Sent: Tuesday, July 03, 2007 5:10 PM
> To: R-help
> Subjec
t; the constraints become equality); and finally, if the solution is an
> interior point, check to see whether the gradient there is close to zero.
> Note that if the solution is one of the vertices of the polyhedron, then the
> gradient may not be zero.
That is a very good idea, Ravi! Thank
m able to get a solution for my problem
with constrOptim, but I am not enough confident that the solution is
right. That is why I am trying to get a solution with rgenoud, but
unsuccessfully until now.
Paul
> -Original Message-
> From: [EMAIL PROTECTED]
> [mailto:[EMAIL PROTECTED]
Dear All,
I am trying to solve the following maximization problem, but I cannot
have rgenoud giving me a reliable solution.
Any ideas?
Thanks in advance,
Paul
library(rgenoud)
v <- 0.90
O1 <- 10
O2 <- 20
O0 <- v*O1+(1-v)*O2
myfunc <- function(x
value, but rgenoud seems not being converging.
Any ideas?
Thanks in advance,
Paul
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and
that does
what I want, and Perl has the wantarray function
detecting the context in which a function is called.
Perhaps match.call() can be made to do what I want,
but, if so, I don't see it in reading the
documentation.
Sincerely,
Paul Laub
__
R-h
So is this just to make the
graph look better, or should I always use -loc for the y-axis of a
similar plot for a completely different data set.
Kind regards,
Paul Lemmens
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r
Hello.
Does anyone out there know if there is a function (or functions) that
will enable me to calculate prediction and confidence intervals from a
non-linear least squares model? (predict.nls does not do this.)
Thank you,
Paul
__
R-help
knows an easier way, then
I'd still be glad to hear from you.
Thanks,
Paul
On 21/06/07, Paul Suckling <[EMAIL PROTECTED]> wrote:
> Hello.
>
> Does anyone out there know if there is a function (or functions) that
> will enable me to calculate prediction and confidence interv
that deals with using grid to modify such
> graphics.
A couple of suggestions:
- look at Hadley Wickham's online book draft for ggplot2
http://had.co.nz/ggplot2/
- look at the online chapter on grid from my book
http://www.stat.auckland.ac.nz/~paul/RGraphics/chapter5.pdf
Paul
> Vikas
o another format as to where
each block could be worked on separately. ie so in one block would be
rows1:10 the next block would be rows11:40 and so on.
Thanks,
Paul
--
Research Technician
Mass Spectrometry
o The
/
o Scripps
\
o Research
/
o Institute
__
ckage might be useful to you. See "Viewing Binary Files
with the hexView Package" in R News 7/1.
Paul
> Todd Remund
>
> __
> R-help@stat.math.ethz.ch mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEAS
(probably after a bit of transformation) to extract the values to very
good accuracy. If you can provide me with an example file, I'd be happy
to play around and see if I could get this to work.
Paul
(Ted Harding) wrote:
> On 15-Jun-07 16:29:53, Ted Harding wrote:
>> [...]
&g
ual level random
error. However, it still appears that there is a substantial place
effect, for Nonwhites only.
Do I understand that right?
Well, thanks in advance, as usual.
--
Paul E. Johnson
Professor, Political Science
1541 Lilac Lane, Room 504
University of Kansas
___
is
more than one editor window open)?
I'm using R 2.4.1 on Windows XP.
Paul
[[alternative HTML version deleted]]
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide
ic to have a good R-implementation of this...
Unfortunately my own skills are several log units below what's required. I
have written to Chow and Ruskey before but unfortunately not heard anything.
Anyone up for the job??
Best wishes
Paul
Nina Hubner wrote:
>
> Hello,
>
>
dim(xs)
[1] 400 4
> xs[1:5,]
x y Z1 Z2
1 27172.4 19062.4 0128
2 27000.9 19077.8 0 0
3 27016.8 19077.5 0 0
4 27029.5 19077.3 0 0
5 27045.4 19077.0 0 0
Cheers,
Paul
--
Research Technician
Mass Spectrometr
On 5/21/07, Alberto Monteiro <[EMAIL PROTECTED]> wrote:
> Paul Lynch wrote:
> >
> > I don't think it makes sense to compare models with
> > and without an intercept term. (Also, I don't know what the point of
> > using a model without an intercept term
On Tue, 22 May 2007, hadley wickham wrote:
> On 5/22/07, Rolf Turner <[EMAIL PROTECTED]> wrote:
>
> > On 22/5/07 6:48 AM, "hadley wickham" <[EMAIL PROTECTED]> wrote:
> >
> > > Two possible solutions:
> > >
> > > * DISPLAY=0.0.0.0:0 R - and then X11() should work without having to use
> > > xterm
ows, it is easy to
construct a "fit" (e.g. y = 10,000,000*x) so that SSR > SST if one is
not deriving the fit from the regular linear regression process.
--Paul
On 5/19/07, 李俊杰 <[EMAIL PROTECTED]> wrote:
> I know that "-1" indicates to remove the inte
,],add=T,axes=F)
axis(side=2,at=c(-150,-100,-50,0,10,20,30))
But I still does not work for the axis??
I would appriciate any kind of hint
Greetings
Paul Magdon
--
___
BSc. Paul Magdon
-Research Assistant-
Institute of Fo
On 5/10/07, Frank E Harrell Jr <[EMAIL PROTECTED]> wrote:
> Paul Johnson wrote:
> > This is a follow up to the message I posted 3 days ago about how to
> > estimate mixed ordinal logit models. I hope you don't mind that I am
> > just pasting in the code and
; > a good idea.
> >
> > 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")
> >
> > Paul Smith wrote:
> >
> >
> >> Dear A
Thanks a lot, Jasjeet. That is it.
Paul
On 5/10/07, Jasjeet Singh Sekhon <[EMAIL PROTECTED]> wrote:
>
> Hi Paul,
>
> I see. You want to increase the population size (pop.size)
> option---of lesser importance are the max.generations,
> wait.generations and P9 options
was
$par
[1] 0.7090278 0.7051806
which is not very precise comparing with the values of the
(analytical) solution. Is it possible to increase the degree of
closeness of the rgenoud solutions with the analytical ones?
Paul
On 5/10/07, Jasjeet Singh Sekhon <[EMAIL PROTECTED]> wrote:
>
?
Thanks in advance,
Paul
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible
d "continuation ratio" logit models before, so I
### looked up a few articles on estimation of ordinal models by
### re-coding the output as a sequence of binary comparisons (stop
### ratios, continuation ratios, etc). The article that is most clear
### on how this can be done to estimate a
Hi
Paul Gilbert wrote:
> Tony
>
> Thanks for the summary.
>
> My ad hoc system is pretty good for catching flagged errors, and
> numerical errors when I have a check. Could you (or someone else)
> comment on how easy it would be with one of these more formal frameworks
otted output that
matters. Is it possible to set up automatic tests of plotting? I can
already test that plots run. I want to know if they "look very
different". And no, I don't have a clue where to start on this one.
Paul Gilbert
[EMAIL PROTECTED] wrote:
> Greetings -
>
>
nce=0.01)
How can one increase the precision of the solution
$par
[1] 0.7072442 0.7069694
?
I have tried solution.tolerance but without a significant improvement.
Any ideas?
Thanks in advance,
Paul
__
R-help@stat.math.ethz.ch mailing list
https://
Thanks, Ravi, for your clear explanation!
Paul
On 5/8/07, RAVI VARADHAN <[EMAIL PROTECTED]> wrote:
> Paul,
>
> The problem lies neither with R nor with numercial methods. The onus is
> always on the user to understand what the numerical schemes can do and what
> they
Ravi Varadhan wrote:
> Paul,
>
> Your solution based on SVD does not work
Ooops. I really am getting rusty. The idea is based on eigenvalues
which, of course, are not always the same as singular values.
Paul
even for the matrix in your example
> (the reason it worked for e=3, wa
may take a wrong solution as a correct
one.
Paul
On 5/8/07, Ravi Varadhan <[EMAIL PROTECTED]> wrote:
> Your function, (x1-x2)^2, has zero gradient at all the starting values such
> that x1 = x2, which means that the gradient-based search methods will
> terminate there because th
on multi-category qualitative dependent variables and
also because the distribution does not seem to fall into the
exponential family. However, some textbooks do include the
multicategory models in the GLM treatment.
--
Paul E. Johnson
Professor, Political Science
1541 Li
On 5/7/07, Paul Smith <[EMAIL PROTECTED]> wrote:
> > I think the problem is the starting point. I do not remember the details
> > of the BFGS method, but I am almost sure the (.5, .5) starting point is
> > suspect, since the abs function is not differentiable at 0.
> the starting point even slightly you will have no problem.
>
> "Paul Smith"
> <[EMAIL PROTECTED]
> > To
> Sent by:
Dear All
I am trying to perform the below optimization problem, but getting
(0.5,0.5) as optimal solution, which is wrong; the correct solution
should be (1,0) or (0,1).
Am I doing something wrong? I am using R 2.5.0 on Fedora Core 6 (Linux).
Thanks in advance,
Paul
3+0.0141531i 0.5287-0.0141531i
[2,] 0.5287-0.0141531i 0.4713+0.0141531i
>
Paul Gilbert
Doran, Harold wrote:
> Suppose I have a square matrix P
>
> P <- matrix(c(.3,.7, .7, .3), ncol=2)
>
> I know that
>
>
>> P * P
>
> Returns the element by element
ws().
... and in grid.xspline() and grid.curve().
Paul
> Thank you.
>
> Paulo Barata
>
> ---
> Paul Murrell wrote:
>> Hi
>>
>>
>> Paulo Barata wrote:
>>> Dr. Snow and Prof. Ripley,
curveto).
In summary: there is some support for smooth curves, but we could still
benefit from a specific arc() function with the standard
centre-radius-angle parameterisation and we could also benefit from
exposing the native strengths of different graphics devices (rather than
the current
1 - 100 of 982 matches
Mail list logo