Gabor Grothendieck gmail.com> writes:
>
> I am not sure what the point of all this is.
The point of all this is:
1. Some users of R here on the list apparently would like to have more powerful
CAS functionalities than Yacas can provide.
2. Many of the solution hints to Ryacas on this list w
Mao Jianfeng gmail.com> writes:
> I am new to R. And, I want to perform a multiple nested anova on a large
> datasets (with 9448 observations).
>
> I tried the function ("lmer()" in package "lme4"). But, I failed. Can anyone
> help me?
>
> my dataset("SeedL.txt") was attached. Data are not sort
Hi Nguyen,
I also experienced this problem, running the development version of R. The
way I fixed it was to uninstall cwhmisc using remove.packages("cwhmisc") and
reinstall it. Note, however, that I manually removed the cwhmisc folder
before reinstalling.
Regards, Mark.
> sessionInfo()
R versio
Marcin Kozak gmail.com> writes:
> How can I deal with truncated labels in the mosaicplot()? Look at the example:
> mosaicplot(~ gear + carb, data = mtcars, color = TRUE)
>
> Look at "carb" (the number of carburetors): the label "8" is
> truncated. How this might be dealt with?
As Michael noted
megh wrote:
>
> suppose I have three vectors like :
>
> l1 = 1:4
> l2 = 4:9
> l3 = 16:67
>
> now I want to construct a loop like :
>
> for (i in 1:3)
>{
> count1[i] = length(li) # i.e. it will take l1, l2, l3 according to
> value of i
>}
>
Try this. There's probably a more ele
Marcin Kozak gmail.com> writes:
>
> How can I deal with truncated labels in the mosaicplot()? Look at the example:
> mosaicplot(~ gear + carb, data = mtcars, color = TRUE)
>
> Look at "carb" (the number of carburetors): the label "8" is
> truncated. How this might be dealt with?
>
Ummm... not
Hi,
suppose I have three vectors like :
l1 = 1:4
l2 = 4:9
l3 = 16:67
now I want to construct a loop like :
for (i in 1:3)
{
count1[i] = length(li) # i.e. it will take l1, l2, l3 according to
value of i
}
Can anyone please tell me how to do that?
Regards,
--
View this message in
You want to do:
temp2 <- matrix(rnorm(10),nc=1,nrow=10)
rownames(temp2) <- c("a","b","c","d","e","f","g","h","i","j")
print(temp2)
o <- order(temp2[,1])
temp2 <- temp2[o,,drop=FALSE]
print(temp2)
Note that it makes no difference if you drop the dimension of a single
column vector or not when pass
try this:
temp2 <- temp2[order(temp2[,1]),] # PROBLEM IS HERE
dim(temp2) <- c(5,1)
print(temp2)
On Wed, Feb 11, 2009 at 2:23 PM, wrote:
> i'm sorry. i had an error in my previous code because i left out a letter in
> the rownames.
> while fixing that, i also found a solution. so i'm sorry fo
i'm sorry. i had an error in my previous code because i left out a
letter in the rownames.
while fixing that, i also found a solution. so i'm sorry for the
confusion.
below is my fix.
temp2 <- matrix(rnorm(10),nc=1,nrow=10)
rownames(temp2) <- c("a","b","c","d","e","f","g","h","i","j")
print(te
this is a bad question but I can't figure it out and i've tried. if i
sort the 2 column
matrix , temp1, by the first column, then things work as expected. But,
if I sort the 1 column matrix, temp2, then it gets turned coerced to a
vector. I realize that I
need to use drop=FALSE but i've put it i
Tinn-R is for MS Windows only.
You may take a look at: http://www.sciviews.org/_rgui/projects/Editors.html
Regards,
Yihui
--
Yihui Xie
Phone: +86-(0)10-82509086 Fax: +86-(0)10-82509086
Mobile: +86-15810805877
Homepage: http://www.yihui.name
School of Statistics, Room 1037, Mingde Main Building,
Welcome to cos.name :-) And thanks for your "advertising", Ronggui!
Regards,
Yihui
--
Yihui Xie
Phone: +86-(0)10-82509086 Fax: +86-(0)10-82509086
Mobile: +86-15810805877
Homepage: http://www.yihui.name
School of Statistics, Room 1037, Mingde Main Building,
Renmin University of China, Beijing, 100
Dear Sir or Madame
At our university we use the program "R" and "Tinn-R", but they don't
know how/where to download this program for Mac-users...could you
please help me? Is Tinn-R for Mac called by an other name?
I'm looking forward to hearing soon from you part.
Thank a lot
Ternia Naduva
How can I deal with truncated labels in the mosaicplot()? Look at the example:
mosaicplot(~ gear + carb, data = mtcars, color = TRUE)
Look at "carb" (the number of carburetors): the label "8" is
truncated. How this might be dealt with?
Thanks in advance,
Marcin
__
Or if the computation is time-consuming, you may try the 'animation'
package, which has provided an HTML/JavaScript approach to show your
graphs one by one, either automatically or manually. See
http://animation.yihui.name/animation:start for details.
Regards,
Yihui
--
Yihui Xie
Phone: +86-(0)10-
Gundala Viswanath gmail.com> writes:
>
> Dear all,
>
> Is there a way to generate K numbers of integer (K = 10^6).
> The maximum value of the integer is 200,000 and minimum is 1.
>
> And the occurrences of this integer follows
> a lognormal distribution.
>
> - Gundala Viswanath
> Jakarta - I
You may read the article in the latest volume of R News:
Hadley Wickham, Michael Lawrence, Duncan Temple Lang, and Deborah F.
Swayne. An introduction to rggobi. R News, 8(2):3-7, October 2008.
Regards,
Yihui
--
Yihui Xie
Phone: +86-(0)10-82509086 Fax: +86-(0)10-82509086
Mobile: +86-15810805877
H
Respectfully, I must disagree. (And it's not my cats, but those of
Turner's colleague.) I particularly disagree with using a Fisher's
exact test as did Turner, as it would double the sample size
improperly (even though the FET is known to be conservative.)
Your strategy appears very much i
This splits them into a 3 column matrix:
> x <- c("ASK", "DGH", "ASG", "AUJ", "URT")
> library(gsubfn)
> z <- strapply(x, ".", simplify = rbind)
> z
[,1] [,2] [,3]
[1,] "A" "S" "K"
[2,] "D" "G" "H"
[3,] "A" "S" "G"
[4,] "A" "U" "J"
[5,] "U" "R" "T"
z[,1], z[,2] and z[,3] are your
-BEGIN PGP SIGNED MESSAGE-
Hash: SHA1
I've written a function in R which takes a symmetrical matrix as input and
processes all triplicate combinations of values from the matrix. The function
looks something like:
my_fun <- function(m) {
if( nrow(mat) != ncol(mat) ) {
stop("'m' must
There are a number of others considerations as well.
Were all the cats given the same dose? If we can establish
that the cats that vomited had a higher dose than the
ones that did not then this would be further evidence. Or if
the cats that vomitted had a lower dose perhaps the drug
is protectiv
Dear all,
Is there a way to generate K numbers of integer (K = 10^6).
The maximum value of the integer is 200,000 and minimum is 1.
And the occurrences of this integer follows
a lognormal distribution.
- Gundala Viswanath
Jakarta - Indonesia
__
R-help
I am not sure what the point of all this is.
This is an R list, not a CAS list. The recommendation to use yacas is based on
the fact that there is an interface between R and yacas. There is no
interface between
R and Maxima so Maxima is not in the running. Anyone who has used
Maxima knows that
David,
If you really want to do a test on this data, I would suggest
a Fisher's Exact test, but you want to use hypergeometric
probabilities. You would probably want to try the CMH
test, if the function allows a single table and actually uses
hypergeometric probabilities.
My suggestion, would be
Dear Larry,
That's odd: The only lmer of which I'm aware is a function in the lme4
package, and it creates objects of class "mer," I believe. What this has to
do with the Rcmdr package is beyond me, since the Rcmdr doesn't depend upon
lme4, but, on the other hand, shouldn't conflict with it. I don
On 11/02/2009, at 1:06 PM, Bert Gunter wrote:
The only question at issue (i.e. capable of being addressed) is: is
giving
the drug to non-vomiting cats associated with vomiting? (I would
strongly
suspect that cats that were vomiting beforehand would have been
excluded
from the study, as the
Thanks.
I seem to have gotten a bit further, but still not completely successful.
BTW, I should have said that I am running R 2.8.1 GUI 1.27 Tiger build
32-bit (5301).
X11 starts fine from the icon.
library(tcltk) loads fine.
Now I get this when I load Rcmdr:
--
> li
There is one in google group (http://groups.google.com/group/R-Forum),
but the traffic is very low. Personally, I don't see the why we need a
Chinese mail list.
Besides, there are two online forum of R, both of which is quite active.
http://rbbs.biosino.org/Rbbs/forums/list.page
http://cos.name/bb
Wow, table() works wonderfully fast! Thank you for pointing it out to me.
I still need to associate those counts with specific
parameter/station/site combinations, and I'm as stumped by that as I am by
the object returned by by(). It looks like I can do the following:
tmp <- with(df1, table(p
Hi Bert: I think what you said about a prior guess for the NULL is
also similar to what Chuck said about people looking with a blank stare.
Thanks for the clarification.
On Tue, Feb 10, 2009 at 7:06 PM, Bert Gunter wrote:
The only question at issue (i.e. capable of being addressed) is:
The only question at issue (i.e. capable of being addressed) is: is giving
the drug to non-vomiting cats associated with vomiting? (I would strongly
suspect that cats that were vomiting beforehand would have been excluded
from the study, as the researcher would have felt that one couldn't then
tell
A colleague from China is visiting our lab is working with us on some
R-based remote sensing projects, and he was curious if there is a r-help
listserv in Chinese and, if not, if it would be possible to set one up?
I didn't see one listed on the Mailing List page on the cran-r website.
Thanks
Still seems that McNemar's test is the appropriate test for the
matched design, but my first answer sent the input to the function
incorrectly, and it's not clear that a normal theory test would be
accurate in all instances. The matrix should have 61 cats with no
vomiting under either situa
Hi,
Is there a way to define the height of an y-axis not only by defining
the tick-marks
but by defining a height in inches or cm?
So that the absolute axis height is the same over different plots...
__
R-help@r-project.org mailing list
https://sta
To quickly see if any duplicates exist you could use table() and
checking for entries that are more than 1. Use na.omit()
on the entire data.frame before passing it to table. E.g.,
tmp <- with(na.omit(df1), table(parameter, station, site))
sum(tmp>1) # number of parameter/station/site keys
Dear Larry,
The tcltk package is part of the standard R distribution; it doesn't live on
CRAN, and should already be installed. Have you tried loading tcltk directly
-- library(tcltk)?
It's possible that you don't have X11 Tcl/Tk installed on your Mac. Take a
look at the Rcmdr installation notes
kayj wrote:
> Hi ,
>
>
> Does any one know how to split a character vector , I have a vector X that
> looks like this and each row has 3 characters
>
> X
> ASK
> DGH
> ASG
> AUJ
> FRT
>
> I would like to split the vector into 3 vectors that look like this
>
> X1X2 X3
> A S K
>
On 11/02/2009, at 11:52 AM, Douglas Bates wrote:
Things are different under Windows.
A fortune candidate?
cheers,
Rolf Turner
##
Attention:\ This e-mail message is privileged and confid..
Hi: Bert: can you do that because the null is that they are equal
before and after,
not that the proportion is zero ? Thank for any clarification to my
lack of understanding.
On Tue, Feb 10, 2009 at 5:43 PM, Bert Gunter wrote:
Ah, experimental units,again ... a subject little taught by
On Sun, Feb 8, 2009 at 1:36 PM, Tom Backer Johnsen wrote:
> Dear me. Is the installation of R under Ubuntu really that complex? I have
> a dual boot machine (Linux / Windows, where I use the latter the most) and
> have plans to try R under Linux, but have not done so yet. Is it possible
> to si
I just installed Rcmd on my OSX-10.5 32-bit powerpc box, with R 2.8.1, and it
seems to work (see below). Beyond that, I can't help, since I've not used
this package. Good luck!
> library(Rcmdr)
Loading required package: tcltk
Loading Tcl/Tk interface ... done
Loading required package: car
Rcm
on 02/10/2009 03:33 PM Rolf Turner wrote:
>
> I am appealing to the general collective wisdom of this
> list in respect of a statistics (rather than R) question. This question
> comes to me from a friend who is a veterinary oncologist. In a study that
> she is writing up there were 73 cats who w
Folks,
I'm checking the structure of a dataframe for duplicate parameters at a
site station (i.e depth should be measured once, not twice), using
aggregate to count each parameter within a site station. The fake data
below has only 26000 rows, and takes roughly 14 seconds. My real data has
7
On Tue, Feb 10, 2009 at 7:51 AM, Gabor Grothendieck
wrote:
> ...Also while Maxima is more sophisticated in terms of algorithms,
Glad to hear it... (I first worked on Maxima in 1971...)
> yacas is actually more sophisticated from the viewpoint of its language which
> borrows ideas from both impe
Ah, experimental units,again ... a subject little taught by statisticians
that is often the crux of the matter. As here.
The cat is the experimental unit. There are 73 of them. 12 of them
experienced vomiting after treatment. What's a confidence interval for the
true proportion based on our sampl
73 cats were treated. None barfing before and 12 after.
This gives the table:
| After | Yes | No | Total |
|+-++---|
| Before Yes | 0 | 0 | 0 |
| Before No | 12 | 61 |73 |
|+-++---|
| Total | 12 | 61 |73 |
and a McNem
Hello.
Somebody knows some R package to work with queues simulation?
Thanks in advance,
Arnau.
Arnau Mir Torres
Edifici A. Turmeda
Campus UIB
Ctra. Valldemossa, km. 7,5
07122 Palma de Mca.
tel: (+34) 971172987
fax: (+34) 971173003
e
Take it in a couple of steps. 'ifelse' will take the evaluation of a
logical vector (first parameter) and return it second parameter if
TRUE or the third parameter if FALSE:
> X<-c(2,2,1,1,0,0)
> X > 0
[1] TRUE TRUE TRUE TRUE FALSE FALSE
> ifelse(X > 0, 1,0)
[1] 1 1 1 1 0 0
>
Now the second
> -Original Message-
> From: r-help-boun...@r-project.org
> [mailto:r-help-boun...@r-project.org] On Behalf Of kayj
> Sent: Tuesday, February 10, 2009 1:44 PM
> To: r-help@r-project.org
> Subject: [R] ifelse()
>
>
> I have a problem with ifelse(), I do not understand how it works.
>
>
kayj asks:
> I have a problem with ifelse(), I do not understand how it works.
>
> > X<-c(2,2,1,1,0,0)
> > str(X)
> num [1:6] 2 2 1 1 0 0
> > Y<-ifelse(X>0,1,0)
> > Y
> [1] 1 1 1 1 0 0
Since X is a vector, the operation X>0 is also a vector. The function
ifelse() is correctly providing output
Did you read the helpfile?
If your condition is true, the first option is returned. If it is false, the
second option is returned.
For the first four elements of X, all of which are greater than zero,
1 is returned. For the last two, which are not greater than zero,
0 is returned.
Sarah
On Tue,
Gabor Grothendieck gmail.com> writes:
>
> Yacas was completely rewritten in java (Ryacas interfaces to the
> C version) since the Ryacas project started so I would not exactly
> characterize yacas as dead. The work that is going on in yacas
> may not have high visibility but that does not mean
Dear Dimitri,
> -Original Message-
> From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org]
On
> Behalf Of Dimitri Liakhovitski
> Sent: February-10-09 4:21 PM
> To: R-Help List
> Subject: [R] Putting values and axis X labels on the charts based on
> allEffects
>
> Dear e
On Tue, Feb 10, 2009 at 01:44:17PM -0800, kayj wrote:
>
> I have a problem with ifelse(), I do not understand how it works.
>
> > X<-c(2,2,1,1,0,0)
> > str(X)
> num [1:6] 2 2 1 1 0 0
> > Y<-ifelse(X>0,1,0)
> > Y
> [1] 1 1 1 1 0 0
> >
>
> Can some one explain what is going on, I do not understa
Hello,
ifelse checks condition whether X (or, ina fact every value in X) is greater
than 0. If it is 1 is assigned and if it isn't then 0.
2 and 1 are greater than 0, therefore the first four values of Y are 1. 0
isn't greater than 0, therefore the last two values of Y are 0.
2009/2/10 kayj
>
It is checking the series of values of X and returning a series of 1
and 0's;
1 if it meets your condition (X>0) and
0 if it doesn't.
What did you expect?
A more complex invocation might look like this:
> ifelse( X > 0, 1:3, -1:-2)
[1] 1 2 3 1 -1 -2
Note the recycling of the elements
Dear Colleagues,
When I try to install Rcmdr the following happens:
--
> library(Rcmdr)
Error in structure(.External("dotTclObjv", objv, PACKAGE = "tcltk"), class =
"tclObj") :
[tcl] invalid command name "font".
Error : .onAttach failed in 'attachNamespace'
Error
I have a problem with ifelse(), I do not understand how it works.
> X<-c(2,2,1,1,0,0)
> str(X)
num [1:6] 2 2 1 1 0 0
> Y<-ifelse(X>0,1,0)
> Y
[1] 1 1 1 1 0 0
>
Can some one explain what is going on, I do not understand what ifelse is
doing in this case. Can someone explain the output Y.
Thank
I have a problem with ifelse(), I do not understand how it works.
> X<-c(2,2,1,1,0,0)
> str(X)
num [1:6] 2 2 1 1 0 0
> Y<-ifelse(X>0,1,0)
> Y
[1] 1 1 1 1 0 0
>
Can some one explain what is going on, I do not understand what ifelse is
doing in this case. Can someone explain the output Y.
Thank
In the biomedical arena, at least as I learned from Rosner's
introductory text, the usual approach to analyzing paired 2 x 2 tables
is McNemar's test.
?mcnemar.test
> mcnemar.test(matrix(c(73,0,61,12),2,2))
McNemar's Chi-squared test with continuity correction
data: matrix(c(73, 0
Hi everyone,
We're very pleased to offer a two-day looking-at-data course.
July 30-31
Washington DC
Day one: static graphics with ggplot2
Day two: interactive graphics with rggobi and GGobi.
You can attend one day (for $295) or both days (for $550).
Student discounts are available.
All proceed
I am appealing to the general collective wisdom of this
list in respect of a statistics (rather than R) question. This question
comes to me from a friend who is a veterinary oncologist. In a study
that
she is writing up there were 73 cats who were treated with a drug called
piroxicam. None
Dear everybody!
Need help with graphics. I am runnig a simple lm and then using
allEffects from 'effects' package:
require(effects)
model<-lm(Y~A+B, data=mydataframe)
I am trying to build (for each predictor - A and then B) a plot of means on Y.
I was successful doing it like this - in one swoop:
Hi all,
I have data from an experiment with 3 independent variables, 2 are
within and 1 is between. In addition to the dependent variable, I have
a covariate that is a single measure per subject. Below I provide an
example generated data set and my approach to implementing the ANCOVA.
However the
Marino, Mark wrote:
Dear R users,
Is there any way to control the size of the box around the mean when
creating a Forest plot using the forestplot function?
H,
> forestplot
Error: object 'forestplot' not found
> PLEASE do read the posting guide
http://www.R-project.org/posting-guide
OK,
obviously that was a silly question...
SOLUTION: just specify the path in the write function...!
Sorry !
Applejus wrote:
>
> Hi,
>
> I have this very simple two-line script :
>
> test.r:
>
> print ("abc");
> write("def","test.txt");
>
> When run with rscript.exe --vanilla test
Hi,
I have this very simple two-line script :
test.r:
print ("abc");
write("def","test.txt");
When run with rscript.exe --vanilla test.r (or just rscript.exe test.r), I
get the "abc" output in the shell but no file is created in the
directory
I noticed that the only time the test.
Henrik,
Ah, I did briefly search the list archive, but hadn't come across your
earlier post, sorry.
I can confirm the conditions you outlined back then (ie. it's L-ALT+TAB, and
only that, that causes the problem, that on return to rterm I can use ctrl-C
to interrupt a running R call but then no f
hi: maybe it's minimizing the negative of the likeihood rather than
maximizing the likelihood ?
other than that, i don't see that being possible because of what you
said. see what
happens if you estimate an arima(1,1,3) ?
On Tue, Feb 10, 2009 at 12:51 PM, Stephen Collins wrote:
All -
I a
Another option is:
read.table(textConnection(gsub("(\\w)", "\\1;", X)), sep = ";")
On Tue, Feb 10, 2009 at 4:50 PM, kayj wrote:
>
> Hi ,
>
>
> Does any one know how to split a character vector , I have a vector X that
> looks like this and each row has 3 characters
>
> X
> ASK
> DGH
> ASG
> AU
kayj wrote:
>
> Hi ,
>
>
> Does any one know how to split a character vector , I have a vector X that
> looks like this and each row has 3 characters
>
> X
> ASK
> DGH
> ASG
> AUJ
> FRT
>
> I would like to split the vector into 3 vectors that look like this
>
> X1X2 X3
> A S
Nguyen Dinh Nguyen garvan.org.au> writes:
>
> Dear Christian,
> Every single time check update package, “cwhmisc” always requests updating.
> I’m aware that the package was latest updated in CRAN on 20Nov2008.
> Is there anything wrong with my R library or somethingelse?
> I use R 2.8.1 on Windo
See ?strsplit
Hope this is helpful,
Dan
Daniel J. Nordlund
Washington State Department of Social and Health Services
Planning, Performance, and Accountability
Research and Data Analysis Division
Olympia, WA 98504-5204
> -Original Message-
> From: r-help-boun...@r-project.org
> [ma
Dear kayj,
Try this:
X<-c("ASK", "DGH", "ASG", "AUJ", "FRT")
res<-data.frame(do.call(rbind,strsplit(X,"")))
colnames(res)<-paste("X",1:ncol(res),sep="")
res
HTH,
Jorge
On Tue, Feb 10, 2009 at 1:50 PM, kayj wrote:
>
> Hi ,
>
>
> Does any one know how to split a character vector , I have a ve
Dear R users,
Is there any way to control the size of the box around the mean when
creating a Forest plot using the forestplot function?
regards
Mark
Mark T. Marino, MD
VP, Early Clinical Development
Mannkind Corp.
61 S. Paramus Road
Paramus, NJ 07652
201-983-5238 Office
203-512-4008
Hi ,
Does any one know how to split a character vector , I have a vector X that
looks like this and each row has 3 characters
X
ASK
DGH
ASG
AUJ
FRT
I would like to split the vector into 3 vectors that look like this
X1 X2 X3
A S K
D G H
A S G
A
Hi,
this sounds awfully similar to what I reported in r-help thread '[R]
Windows Vista, Rterm & LeftAlt + Tab issue' on Aug 21, 2008, cf.
https://stat.ethz.ch/pipermail/r-help/2008-August/171548.html
At least we are not alone ;) ...and I though it was an issue with
Windows Vista only, because
Yes you can run the Rserve application on a remote computer in e.g network
and then you can connect
it from your Desktop computer. I've done this already with a remote Linux
computer (Rserve is running
very stable on a Linux computer).
For example i clustered an R,G,B image on a linux computer wit
Harsh wrote:
At the cost of sounding far-fetched and almost incredulous, I would
like to know if any R user is remotely considering the use of R on
Mobile devices, and Android in particular.
In addition to the other objections in the other replies, I add this
one: the assignment operator and
I've been trying to use Rterm.exe, mainly because of Rgui.exe's lack of
support for incremental history searching. But on my system (Win XP, SP3),
it (rterm from R-2.8.1) frequently just stops responding to keyboard input.
This happens often (enough to make it unusable), and generally after a
switc
Have you checked that that time exists in the time zone you are using? From
?ISOdatetime:
Note
...
Remember that in most timezones some times do not occur and some occur twice because of transitions to/from summer time. What happens in those cases is OS-specific.
You could try working out wh
All -
I am evaluating an arima(2,1,3) and arima(3,1,3) and notice the
log-likelihood of the restricted model is higher than the log-likelihood
of the unrestricted.
Since these are nested models, I thought the unrestricted model would have
a log-likelihood at least as large as that of the res
I, for one, don't fully understand the question, and while there are many on
this list that are much smarter than I am, they may also have trouble
understanding the question.
At first glance it looks like you are asking if R runs on Windows XP (the
answer is yes, but this seems to simple to be
Did you try the forums on Knime? That might be a better bet.
-Original Message-
From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org]
On Behalf Of Antje
Sent: Tuesday, February 10, 2009 12:24 PM
To: r-h...@stat.math.ethz.ch
Subject: Re: [R] R running as a server on wind
Antje wrote:
Nobody can help with this question (or tell me where I can find help) ?
Antje schrieb:
Hi there,
I'm not sure whether this is the right mailing list to put this
question... if not, please tell me where I can find help.
I have a spare dual-core machine which run Windows XP Pro
On Tue, 10 Feb 2009, Keith Jewell wrote:
sapply(trgroup, function(x) which(x>=starts & x <=ends))
works for your example data, but there's probably a better way
For large scale problems this works well:
http://finzi.psych.upenn.edu/R/Rhelp02a/archive/60815.html
HTH,
Chuck
HT
In this case, following your instructions and a bit of noodling with
the example "solves" the problem with:
persp(teta1, teta2, matrixMuSigma) -> res
points(trans3d( teta1, sigmaMax, apply(matrixMuSigma, 1, max),
pmat=res), col="red")
--
David Winsemius
On Feb 10, 2009, at 11:50 AM, Uwe Li
sapply(trgroup, function(x) which(x>=starts & x <=ends))
works for your example data, but there's probably a better way
HTH
KJ
"T Joshi" wrote in message
news:11417a880902100902w53664a3dq11aee64a963d8...@mail.gmail.com...
> Hi,
> I have scenario in which I wish to check whether numeric values
Nobody can help with this question (or tell me where I can find help) ?
Antje schrieb:
Hi there,
I'm not sure whether this is the right mailing list to put this
question... if not, please tell me where I can find help.
I have a spare dual-core machine which run Windows XP Professional and
i
on 02/10/2009 11:02 AM T Joshi wrote:
> Hi,
> I have scenario in which I wish to check whether numeric values in one array
> falls within the range of numbers defined over corresponding values in two
> other vectors:
>
> starts = c(12,45,67,110)
> ends=c(24, 58,102,150)
>
> trgroup=c(18,87)
>
>
or simultaneously estimate the coefficients and variance structure via
nlme::gls and its 'weights' argument...
On Tue, Feb 10, 2009 at 7:57 AM, John Fox wrote:
> Dear Kishore,
>
> Yes, White's heteroscedasticity-consistent standard errors are just that --
> standard errors for the OLS coefficient
?predict
On Tue, Feb 10, 2009 at 11:06 AM, lauramorg...@bluewin.ch
wrote:
> Hello,
> to plot the result of a singular non linear regression (using nls) I usually
> use the function plotfit, for example:
>
> r.PTG.V<-nls(PTG.P~ fz1(Portata, a,b), data=dati, start=list(a=10, b=10),
> nls.control(
lauramorgana bluewin.ch bluewin.ch> writes:
>
> Hello,
> I'm using the following for loop to find regression curves using a list of
functions (formList), a list of
.. long non-reproducible code removed
> And I got the following error message:
> Error in chol2inv(object$m$Rmat()) : l'elemento
Hi,
I have scenario in which I wish to check whether numeric values in one array
falls within the range of numbers defined over corresponding values in two
other vectors:
starts = c(12,45,67,110)
ends=c(24, 58,102,150)
trgroup=c(18,87)
The result should be "1,3" , indices of vector starts/ends.
See ?persp, particularly the example that starts with
# (2) Add to existing persp plot - using trans3d() :
Uwe Ligges
Domenico Vistocco wrote:
Dear All,
is there a way to superimpose points and/or lines on a surface plot?
Below I try to explain my problem.
Suppose I have the following surfac
My guess is that your two machines have different setting for locale.
What does this produce on each of them:
sessionInfo()$locale
(Note: It would have been part of the information that you were asked
to provide per the posting guide.)
---
For questions about unexpected behavior or a possib
Thanks a lot. I will need to change the time to GMT, then, I guess (this is
indeed GMT time).
Pedro
Daylight savings time, I guess. The clock went directly from 1:59 to
3:00 that day. Other systems behave strangely (er, "handle the anomaly")
in different ways:
> ISOdatetime(1995,03,26,1,10,0)
Hi All,
I am watching a strange behaviour of ISOdatetime. In my work
computer, I get NA when I try to do
> ISOdatetime(1995,03,26,2,0,0)
[1] NA
But on other dates and/or times (hour) works OK
> ISOdatetime(1995,03,25,2,0,0)
[1] "1995-03-25 02:00:00 GMT"
In my home computer, I do not have this
Thanks. Now I know how to handle it. It was killing me
Pedro
Prof Brian Ripley wrote:
>
> On Tue, 10 Feb 2009, pbarros wrote:
>
>>
>> Hi all,
>>
>> I am using ISOdatetime, and I just found out that when I do
>> ISOdatetime(1995,03,26,2,10,0) (or any other value under minutes or
>> seconds)
1 - 100 of 156 matches
Mail list logo
