[R] Planned comparison ANOVA

2011-09-17 Thread robcinm
I am trying to do a priori ANOVA analysis for a class assignment. The professor uses SPSS and does not know R. I want to do a simple planned comparison but have been unable to find a function or specific help. There is a grouping variable with five levels and a subsequent response variable. I also

[R] "graph"

2011-09-17 Thread Sandy Mitchell
Hello, I have downloaded "graph" package from http://www.bioconductor.org/packages/release/bioc/html/graph.html Is there anyone who can teach me how to install "graph" into R? Thanks [[alternative HTML version deleted]] _

[R] Could R run Cox proportional hazard model with multiple failure time data

2011-09-17 Thread XUT
If my survival data includes more than two outcomes, for example death, disability and censored. Does R have corresponding packages or programes to run Cox proportional hazard model with multiple failure time data. Millions of thanks!! -- View this message in context: http://r.789695.n4.nabble.co

Re: [R] Tobit Fixed Effects

2011-09-17 Thread Arne Henningsen
On 18 September 2011 04:58, Felipe Nunes wrote: > Do you guys know how could I increase the time limit for the iterations in > censReg? I'm receiving the following message whatever method I use: > Newton-Raphson maximisation, 150 iterations > Return code 4: Iteration limit exceeded. > Log-likeliho

Re: [R] Tobit Fixed Effects

2011-09-17 Thread Felipe Nunes
Do you guys know how could I increase the time limit for the iterations in censReg? I'm receiving the following message whatever method I use: Newton-Raphson maximisation, 150 iterations Return code 4: Iteration limit exceeded. Log-likelihood: -67680.41 on 42 Df thanks, *Felipe Nunes* CAPES/Fulb

Re: [R] Plot history

2011-09-17 Thread R. Michael Weylandt
Could you repaste your data using dput() so we can see exactly what we are working with here. As things stand currently, I can't tell why it doesn't suffice to just combine the dates and times into POSIX objects and subtract them all from the first row. For more information on that approach, try ?

Re: [R] Plot history

2011-09-17 Thread Ayyappa Chaturvedula
Dear Group, Can you help me solve the problem with calculating Time Post dose using date time fomat? Here is information for subject 1 and I appreciate your help for code to generate time post dose from the first dose in each individual. Thank you in advance. Regards, Ayyappa ID DA

Re: [R] Replacing matching values by related values

2011-09-17 Thread R. Michael Weylandt
Try playing with match(). Something like d[match(v,d[,1]),2] Should work (untested bc I'm writing from my phone though) Michael Weylandt On Sep 17, 2011, at 4:33 PM, "Janssen, K.J.M." wrote: > > I am trying to replace values of a vector (consisting of 15 values) by a > value that is relate

Re: [R] apply question

2011-09-17 Thread Jorge I Velez
Hi clsnyder, You forgot the MARGIN argument in you function call. Take a look at ?apply apply(iris, 2, length) apply(iris, 2, mean) HTH, Jorge * * On Sat, Sep 17, 2011 at 6:42 PM, clsnyder <> wrote: > Why does one method work, and the other not? > > attach(iris) > head(iris) > > mean(Sepal.L

[R] apply question

2011-09-17 Thread clsnyder
Why does one method work, and the other not? attach(iris) head(iris) mean(Sepal.Length) [1] 5.84 apply(iris, Sepal.Length, mean) "Error in if (d2 == 0L) { : missing value where TRUE/FALSE needed" I seem to get this error frequently with the apply functions but other "manual" methods work no

[R] Replacing matching values by related values

2011-09-17 Thread Janssen, K.J.M.
I am trying to replace values of a vector (consisting of 15 values) by a value that is related to a matching value in a dataset (consisting of 17 rows). Here's an example The vector: v <- c(f,a,e,d,m,o,e,f,i,n,e,i,b,a,o) The dataset's columns consist of the following values d[,1] <- c(a,b,c,d,e

Re: [R] £50 for help in my masters dissertation kalman filter forecasting

2011-09-17 Thread nclfinance
Dear R users, Just to clarify. I am not offering to pay someone to do my Dissertation. These 4-5 commands on Kalman Filter would be only a tiny part of my 10,000 words dissertation. A part that even after trying for a few days, I am still stuck on. I am offering £50, just to say thanks. Regards

Re: [R] completing missing samples

2011-09-17 Thread David Winsemius
On Sep 17, 2011, at 5:04 PM, Joshua Wiley wrote: On Sat, Sep 17, 2011 at 1:43 PM, Joshua Wiley wrote: Hi Eran, You have already gotten some suggestions from Michael, but I think that Rich is correct to question the rational. I hope we should question the rationale, and doubt that Rich has

Re: [R] completing missing samples

2011-09-17 Thread Joshua Wiley
On Sat, Sep 17, 2011 at 1:43 PM, Joshua Wiley wrote: > Hi Eran, > > You have already gotten some suggestions from Michael, but I think > that Rich is correct to question the rational.  Any mechanism you > choose to replace the missing values will impose its structure on the > data. Veritate ab abs

Re: [R] completing missing samples

2011-09-17 Thread Joshua Wiley
Hi Eran, You have already gotten some suggestions from Michael, but I think that Rich is correct to question the rational. Any mechanism you choose to replace the missing values will impose its structure on the data. Veritate ab absurdo: ## data x <- sin(seq(1, 17, .1)) + seq(-.5, .5, length.out

Re: [R] onet.permutation()

2011-09-17 Thread song_gpqg
Oh I see the problem. Many thanks! -- View this message in context: http://r.789695.n4.nabble.com/onet-permutation-tp3820372p3820721.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org mailing list https://stat.ethz.c

[R] R kiteChart (plotrix) arguments

2011-09-17 Thread antqueen
Hey, I hope someone is able to help. I've created a graph using kiteChart (in the plotrix package) but need to know how to perform some basic text alterations to the graph. I would like to re-orientate the text on the y-axis from vertical to horizontal. I understand that the argument normally

Re: [R] completing missing samples

2011-09-17 Thread Rich Shepard
On Sat, 17 Sep 2011, R. Michael Weylandt wrote: How are your time samples missing? Michael/Eran, Pardon me for jumping in since I did not read the original message, but I have a more fundamental question: are the missing time samples meaningful? For example, a time series of precipitatio

Re: [R] completing missing samples

2011-09-17 Thread R. Michael Weylandt
How are your time samples missing? If they are recorded as NA, the na.locf() function will fill them with the previous value (zero-order hold) and with the reversability arguments can give linear interpolation: library(xts) x = c(1:5,NA,6:10) x = xts(x,Sys.Date()+0:10) na.locf(x) (na.locf(x) + n

Re: [R] implicit data frame reference

2011-09-17 Thread R. Michael Weylandt
This one seems to have slipped through the cracks, but I'll take a stab at it: I don't think there's an (easy) general solution but there are two direct ways to get at the idea. If you are generally only seeking to do things like: x[order(x)] would sort(x) work for you? More generally, this can

Re: [R] Referring to an object by a variable containing its name: 6 failures

2011-09-17 Thread Duncan Murdoch
On 17/09/2011 2:15 PM, Joshua Wiley wrote: Hi Andrew, On Sat, Sep 17, 2011 at 12:04 AM, andrewH wrote: > Dear Folks -- > The anonymous poster ("rmailbox") is perfectly correct. I had forgotten you > could use names in this way. When referring to rows or columns by name > rather than by nu

Re: [R] Extracting a a chunk of text from a pdf file

2011-09-17 Thread Jeff Newmiller
Doesn't seen like an R task, but see pdf2text? (From pdftools, UNIX command line tools) --- Jeff Newmiller The . . Go Live... DCN: Basics: ##.#. ##.#. Live Go... Live: OO#.. Dead: OO#.. Playing Research Engineer (Solar

Re: [R] onet.permutation()

2011-09-17 Thread B77S
** I am not a statistician, but... The default arguments are: onet.permutation(x, nsim=2000, plotit=TRUE) you are making nsim=4 (this parameter is not the size of each sample, it is the number of samples taken) You are building the distribution for a statistical test through random sampling, s

Re: [R] Referring to an object by a variable containing its name: 6 failures

2011-09-17 Thread Joshua Wiley
Hi Andrew, On Sat, Sep 17, 2011 at 12:04 AM, andrewH wrote: > Dear Folks -- > The anonymous poster ("rmailbox") is perfectly correct.  I had forgotten you > could use names in this way.  When referring to rows or columns by name > rather than by number, I usually use either attach() or the $ oper

[R] Extracting a a chunk of text from a pdf file

2011-09-17 Thread Victor
In an R script I need to extract some figures from many web pages in pdf format. As an example see http://www.terna.it/LinkClick.aspx?fileticket=TTQuOPUf%2fs0%3d&tabid=435&mid=3072 from which I would like to extract the "Totale: 1,025,823"). Is there any solution? Ciao Vittorio __

[R] onet.permutation()

2011-09-17 Thread song_gpqg
I saw the manual of this function but not sure what to do. I have a array contain 7 numbers and want to choose 4 to do permutation test. But using this function with parameters as onet.permutation(scores,4), it returns 0. Instead, with no parameter, onet.permutation(), it returns something but eve

[R] Warning in 'probtrans'-function ('mstate'-package)

2011-09-17 Thread Johannes Johow
Dear all, in order to estimate transition-specific probabilities in a multi-state model i applied the 'probtrans()' function from the 'mstate'-package. Now, i am at loss with the following message (see attached example): Warning message: In probtrans(msf.0, predt = 0) : Negative diagonal elemen

[R] Fw: clustering with finite mixture model

2011-09-17 Thread karuna m
Dear R-help, I am doing clustering via finite mixture model. I am trying with mixtools. Please suggest some packages in R to find clusters via finite mixture model with continuous variables. And also I wish to verify model fit and the distributional properties of the mixture distributions by

[R] Fw: likelihood ratio test

2011-09-17 Thread karuna m
Ms.Karunambigai M PhD Scholar Dept. of Biostatistics NIMHANS Bangalore India - Forwarded Message From: karuna m To: R-help@r-project.org Sent: Thu, 26 May, 2011 10:55:02 AM Subject: likelihood ratio test Dear R-help, Can anybody tell me which R package has Lo-Mendell Rubin LR test

Re: [R] cutree() and rect.hclust(): different labelling of classes

2011-09-17 Thread Rui Tian
Hi, Agustin, I also encountered similar problems when cutting trees. However, I think one himself/ herself has to figure out the components of each cluster. The labels that are generated by the current packages simply only distinguish one cluster from another, but do not have more meanings. Rui O

Re: [R] Weibull point process

2011-09-17 Thread Torbjørn Ergon
On Thu, 15 Sep 2011 22:12:54 +1200, Rolf Turner wrote: On 15/09/11 19:24, Torbjørn Ergon wrote: On Thu, 15 Sep 2011 09:47:35 +1200, Rolf Turner wrote: On 15/09/11 07:21, Torbjørn Ergon wrote: Dear list, I'm looking for a function to generate (simulate) a random Weibull point process. Can a

[R] Problem using SNOW with data frame as a function argument

2011-09-17 Thread Paul.Rustomji
Hello I would like to use SNOW to parallelise some computations to be made on columns of a data frame, using different parameter values for each SNOW "worker". I gather(?) clusterMap() is the appropriate SNOW function to do something like this. I suspect the problem lies in the fact that I am

Re: [R] Name the dots! ("...")

2011-09-17 Thread Duncan Murdoch
On 11-09-17 5:27 AM, andrewH wrote: Dear Folk-- Suppose I have some objects A, B& C, and a function getDots<- function(...) {args<- list(...) etc.} If I do a call to getDots(A, B, C) then the variable args will be assigned to a list which contains the objects to which A, B& C refer, but whic

Re: [R] question about glm vs. loglin()

2011-09-17 Thread Ted Harding
On 17-Sep-11 01:20:53, Rolf Turner wrote: > On 17/09/11 01:19, Michael Friendly wrote: > > >> Trying to interpret associations in complex loglinear models >> from tables of parameter estimates is like trying to extract >> sunlight from a cucumber. You have to squeeze very hard, and >> then are

[R] Name the dots! ("...")

2011-09-17 Thread andrewH
Dear Folk-- Suppose I have some objects A, B & C, and a function getDots <- function(...) {args <- list(...) etc.} If I do a call to getDots(A, B, C) then the variable args will be assigned to a list which contains the objects to which A, B & C refer, but which will not (except by happenstance)

Re: [R] Referring to an object by a variable containing its name: 6 failures

2011-09-17 Thread andrewH
Dear Folks -- The anonymous poster ("rmailbox") is perfectly correct. I had forgotten you could use names in this way. When referring to rows or columns by name rather than by number, I usually use either attach() or the $ operator, neither of which works here. If anyone understands why data.df[