[R] Running R on a hosting o server

(8) - daemon to execute scheduled commands (Vixie Cron) So you schedule your task. -- Marco Arthur @ (M)arco Creatives __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEA

Re: [R] replace NA's with row means for specific columns

tc..) and the row id itself (so consider pass the left hand of assignment entirely). I would use: fmean <- function(row, col_selection) { # homework for you here } Best Regards, -- Marco Arthur @ (M)arco Creatives __ R-help@r-project.org mailing

Re: [R] problem in R

gt; > to programming in R. I see no future in this enterprise. > > > > *Gotta* be a fortune!!! > > cheers, > > Rolf > Jim++ Rolf++ must be. -- Marco Arthur @ (M)arco Creatives __ R-help@r-project.org mailing list -- To

[R] Dichotomous variables

y for this long post marco -- View this message in context: http://r.789695.n4.nabble.com/Dichotomous-variables-tp3721550p3721550.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org mailing list https://stat.ethz.ch/ma

Re: [R] Dichotomous variables

oes not solve the problem: same message. So: this is the blind alley. thanks for any feedback marco -- View this message in context: http://r.789695.n4.nabble.com/Dichotomous-variables-tp3721550p3724639.html Sent from the R help mailing list archive at Nabble.com. __

[R] converting matrix in array

Hi everyone, have a small problem trying to converting a dataset in matrix form to an array. Specifically: data include 3D measurement -x,y,z of 59 points in 36 objects. They are stored as a matrix (x) of 2124 rows and 3 columns. What I want to do is to extract each subject's dataset using an arra

[R] How to merge a matrix and a dataframe with different types of columns

7-15" [14,] "2011-07-16" [15,] "2011-07-17" [16,] "2011-07-18" [17,] "2011-07-19" [18,] "2011-07-20" [19,] "2011-07-21" [20,] "2011-07-22" [21,] NA [22,] NA [23,] NA

[R] How to merge a matrix and a dataframe with different types of columns

7-15" [14,] "2011-07-16" [15,] "2011-07-17" [16,] "2011-07-18" [17,] "2011-07-19" [18,] "2011-07-20" [19,] "2011-07-21" [20,] "2011-07-22" [21,] NA [22,] NA [23,] NA

[R] smallest space analysis

Hi everyone just a brief question: I'm trying to perform a smallest space analysis on nonmetric variables with R. It seems there are at lear a couple of packages to perform non metric MDS but still not found any info concerning SSA. Any clue? thanks in adv. Marco -- View this message in co

[R] smooth scatterplot and geo map

n the map, I would like to have a "smoothed" scatterplot of the same superimposed on the map. Tried with the smoothScatter function, but really have no idea how to combine the map with the resulting graph. tx in adv marco -- View this message in context: http://r.789695.n4.nabble.com/smo

Re: [R] smooth scatterplot and geo map

thanks a lot!! -- View this message in context: http://r.789695.n4.nabble.com/smooth-scatterplot-and-geo-map-tp3702374p3702424.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org mailing list https://stat.ethz.ch/mail

[R] SSA in R

Hi everyone I'm here with a brief question: Does anyone have some idea on how to perform Smallest Space Analysis in R? I'm aware of several packages/methods for non metric multidimensional scaling, but no clue about this specific technique. Thanks for any feedback Marco -- View th

[R] [R-pkgs] lmvar

arry out cross-validations. The vignette 'Introduction to the package' provides an introduction, the vignette 'Math details' provides mathematical background. Version 1.2.1 is now available. Please send bugs, requests and comments to the package maintainer. Regards,

Re: [R] bnlearn and cpquery

distributions from the data or 2) use posterior estimates for the parameters. Cheers, Marco On 13 July 2017 at 06:29, Ross Chapman wrote: > Hi all > > > > I have built a Bayesian network using discrete data using the bnlearn > package. > > > > When I try to run the cp

Re: [R] I have been unable to run R from within QGIS on my Mac

din your system, or it is not correctly configured to be > used from QGIS > Click here to know more about how to install and configure R to be > used with QGIS" > > [[alternative HTML version deleted]] > Maybe you should ask in a QGIS forum, your QGIS question. But, j

[R] bnlearn cpquery

rvals in the evidence. Cheers, Marco -- Marco Scutari, Ph.D. Lecturer in Statistics, Department of Statistics University of Oxford, United Kingdom [[alternative HTML version deleted]] __ R-help@r-project.org mailing list -- To UNSUBSCRIB

Re: [R] cpquery problem

that these classes should both typically > return a value for ABW that is very much higher than the threshold value. That may be, but much depends on the specific sample the model was fitted from. How does the fitted network look like? Cheers, Marco -- Marco Scutari, Ph.D. L

Re: [R] cpquery problem

7;y', TR = c(9, max(data$TR)), BU = c(15819, max(data$BU)), RF = c(2989, max(data$RF)), n=10^6, method = "lw") 3) look at the parameters in your fitted network and diagnose why this is happening. Cheers, Marco -- Marco Scutari, Ph.D. Lecturer in Statistics, Department

[R] type of objects in last session

t;"character""character" op rain t t2 ts_ "character""character""character""character""character" And I noticed everything is type "cha

Re: [R] cpquery problem

an average of ~= 13 without considering FFB. If I assume FFB is positive, then I can easily see E(y) ~= 15 and E(y) - 1.96 * 0.96 s.d. ~= 13.5. So ABW < 11 has zero or almost zero probability mass. Cheers, Marco -- Marco Scutari, Ph.D. Lecturer in Statistics, Department of Statistic

Re: [R] using a regular expression

s for deletion of the role pattern space, so if you are trying to substitute use: s/pattern//g # effectively removing pattern from the text Best Luck, Marco -- Marco Arthur @ (M)arco Creatives __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and

[R] Mixed-effects model for pre-post randomization design

ite where it is explicitly reported to be "confidential"), so I have to review a bit of likelihood theory along with some implementation details. In the meantime, I wonder if anyone here could point out any related documentation to me. Thank you. Marco. __

Re: [R] help in anova

I suspect that your "treatment" levels are perfectly nested inside the "soilgroup" levels. If this were the case you can either use one of the factor in your analysis or the other, depending on what you want to analyze. Il giorno Wed, 11 Feb 2015 14:47:18 + Mir Salam ha scritto: > Dear all,

Re: [R] Mixed-effects model for pre-post randomization design

Thank you very much, Ben. Your answer has been very useful. __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and

[R] Siegel nonparametric regression / mblm package

standard regression and Theil-Sen regression are practically identical. But the Siegel regression seems to have a bias that I cannot understand. May I ask for a possible explanation? The bias may be related to the number of ties in the set of data? Here's the code and the image. Best regards. M

Re: [R] Siegel nonparametric regression / mblm package

that we haven’t overlooked > something more basic, is that such > high breakdown estimates sacrifice some efficiency, that is to say, they > are more variable than other methods > when the data is well behaved, and of course, the Galton data is famously > “almost Gaussian”. > > On Feb

[R] coplot: extract regression output

cex=1, number=3, columns=3, col="black", pch=16, overlap = 0.1, panel=function(x,y,...) { points(x,y,pch=16) panel.smooth(x,y,span=.8,iter=5,...) abline(lm(y ~ x), col="blue")} ) Thank you. Marco __ R-help@r-project.org mailing li

Re: [R] Strange error while passing string as an argument to the function in bnlearn package

e(text="(M=='s')")), evidence=list(lag1.M1='s'), method = "lw") passing str2 as a list. Cheers, Marco -- Marco Scutari, Ph.D. Lecturer in Statistics, Department of S

Re: [R] Strange error while passing string as an argument to the function in bnlearn package

ot;method")= chr "lw" To compute P(A = a | whatever you conditioned on), just sum the corresponding weights over the total weight mass. Since no language trickery is involved, this works reliably. Cheers, Marco -- Marco Scutari, Ph.D. Lecturer in Statistics, Department of Sta

[R] Metafor and forest(); not showing 'ilab' and text

Hello folks, I have a couple of issues with the metafor package, specifically with the forest graphs. I am currently conducting a Meta-Analysis in economics throughout the metafor package. My meta-analysis has the specific of having different cases from single studies, and this proven to be chall

Re: [R] Metafor and forest(); not showing 'ilab' and text

ce and your time, and many compliments for the package, is guiding me through the use of R for the first time - as you might have guessed. Marco On 24 August 2015 at 16:50, Viechtbauer Wolfgang (STAT) wrote: > I cannot reproduce the issue with 'ilab' not being shown when using '

Re: [R] Metafor and forest(); not showing 'ilab' and text

ssage- >> From: Michael Dewey [mailto:li...@dewey.myzen.co.uk] >> Sent: Tuesday, August 25, 2015 13:23 >> To: Marco Colagrossi; Viechtbauer Wolfgang (STAT) >> Cc: r-help@r-project.org >> Subject: Re: [R] Metafor and forest(); not showing 'ilab' and text &

Re: [R] Metafor and forest(); not showing 'ilab' and text

e there the formulas for it somewhere? So far I only found the formulas for the estimators included in the metafor package. 3. metareg1 <- rma.mv(pc, var, random = ~ 1 | author, mods = ~ pub + SIMiv, data=codebook) Again, if I'm correct this should be a multilevel meta regression (correct me

Re: [R] Metafor and forest(); not showing 'ilab' and text

The solution that you proposed works perfectly, thank you very much. I'll wait for Wolfgang answer as I'm having few doubts about the models. Thanks On 31 August 2015 at 18:34, Michael Dewey wrote: > Comments in line > > On 31/08/2015 16:08, Marco Colagrossi wrote: >>

Re: [R] Metafor and forest(); not showing 'ilab' and text

and correct that with struct="HAR" ? Sorry for the long mail, I wanted to explain as better as possible, and thank you again for your help, It is incredibly appreciated. Marco On 31 August 2015 at 19:28, Viechtbauer Wolfgang (STAT) wrote: > Have you read help(rma.mv)? It describes i

[R] find all unit vectors which are orthogonal to a given vector

Is there a function in R to get the set of all unit vectors which are orthogonal to a given vector? __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http:

Re: [R] find all unit vectors which are orthogonal to a given vector

Ok, thanks! On 2015-10-29 09:03 PM, Duncan Murdoch wrote: On 29/10/2015 6:38 PM, Marco Inacio wrote: Is there a function in R to get the set of all unit vectors which are orthogonal to a given vector? No. Duncan Murdoch __ R-help@r-project.org

Re: [R] find all unit vectors which are orthogonal to a given vector

I think that's just what I needed. Thanks! On 2015-10-29 09:55 PM, peter dalgaard wrote: On 30 Oct 2015, at 00:03 , Duncan Murdoch wrote: On 29/10/2015 6:38 PM, Marco Inacio wrote: Is there a function in R to get the set of all unit vectors which are orthogonal to a given vector

Re: [R] Query in R.

als first. Cheers, Marco -- Marco Scutari, Ph.D. Lecturer in Statistics, Department of Statistics University of Oxford, United Kingdom __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEAS

Re: [R] bnlearn and TAN network

produce a graph + barplots plot from Rgraphviz. I have never managed to make it work, though. Cheers, Marco -- Marco Scutari, Ph.D. Lecturer in Statistics, Department of Statistics University of Oxford, United Kingdom __ R-help@r-project.org mai

Re: [R] bnlearn and TAN network

aining? TAN does not do any sort of feature selection, so all the nodes should be there. As for "showing the network after training", what kind of plot are you looking for? Cheers, Marco -- Marco Scutari, Ph.D. Lecturer in Statistics, Department of S

[R] Expand 80 columns terminal output

Hi R'ers, I noted that R breaks lines for any output greater than 80 chars, Is there a way to expand this limit ? Thks -- Marco Arthur @ (M)arco Creatives __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.et

Re: [R] web scraping tables generated in multiple server pages

> 4700 Keele StreetWeb:http://www.datavis.ca > Toronto, ONT M3J 1P3 CANADA > > # what is missing to you ?gsub # aliasing df <- sp.tab[, wanted] # convert to double as.double( # convert to double gsub(',', '.',

Re: [R] bnlearn impute problem

ve to upgrade to use it. Cheers, Marco -- Marco Scutari, Ph.D. Lecturer in Statistics, Department of Statistics University of Oxford, United Kingdom __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman

[R] Graph

| -- - | 7.01% | 78 | - Marco __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the

[R] chisq

Hi, the chisq value for .05 level of signif. and 3 DF is 7.81 (one-sided). Is there a way to calculated it in R (giving level of signif. and DF)? Marco __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the

[R] expression

t(paste(expression(paste("R"^2)),round(marco2[1,i],digits=3)," N° of proteins:",marco3[i]),side=4,cex=.6) the output is: paste("R"^2) Thank you in advance Marco / __ R-help@r-project.org mailing list https://stat.ethz.ch/m

[R] hier.part

joint all the vectors into a unique data frame (xcan in the usage of hier.part) ? I´ve been trying to run it, but i keep getting different errors. Thanks in advance, Marco [[alternative HTML version deleted]] __ R-help@r-project.org mailing list

[R] error_hier.part

one please give me some clues on what that is? Thanks in advance Marco [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.o

[R] why this function does not run correctly?

Hi, my name is Marco Bressan i'm working to improve ADati package. I study psicology ad Padua University (Italy). I have this problem: why bartlett.test function running good and my anova.welch function no? Ciao, il mio nome è Marco Bressan e sto lavorando per migliorare il pacchetto

[R] Following progress in a lapply() function

(x) return(list(p.value=tTest$p.value,t.stat=tTest$stat)) } Using an lapply statement, as in: myResults <- lapply(aList, myFnc) The goal would be to report on the progress of the lapply() function during processing. Any suggestion would be greatly appreciated. Thanks Marco -- Marco

[R] calculate weeks for 2005 and aggregate them

<- cbind(date, hcrime) I have following questions: a) how can I find out the weeks of every date b) aggregate the dates falling in the same week c) sum up the different crimes per week for better understanding here is an example for a) and b): dateweek 27-08-2005

[R] apply command

control if the result of the apply is correct. Do anyone find a solution to this issue? thanks Marco [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting

[R] CRAN Packages By Date

Hi, why the list of "Available CRAN Packages By Date of Publication" is stopped at 2022-12-25 ? Regards Marco Scazzocchio [[alternative HTML version deleted]] __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more,

[R] clusplot() and fviz_cluster() functions

t=50)#library(cluster)windows()clusplot(z, fit$cluster, labels=3, lines=0, cex=0.6, col.txt="black", col.p="black")#library(factoextra)windows()fviz_cluster(fit, data=z)#* Thanks if someone has a suggestion! ----

[R] Error with stepwise models selection using NADA library

Hi, I'm working with environmental data which are left-censored and I found the R package NADA which seems to do the job. After fitting a complete model, using the cenreg function,I'd like to do a stepwise procedure to minimize the AIC using the stepAIC function. Unfortunately, after fitting the m

[R] Help with unzip error -1

I'm running Rx64 3.01 on Windows 7 and until the other day I could easily unzip zip files with unzip (package utils 3.01). Since yesterday however I get an error. Let's say my zipfile file path is called zfile (I'm omitting the output directory for this example). When I try to unzip it I get the fo

[R] IRanges::unlist in package

tried using only Imports: as suggested by Chambers but it breaks. Using Depends does not help. Am I having clashing namespace? Here is my Depends: (or Imports:) line: Depends: Rsamtools, GenomicFeatures, parallel, rtracklayer, edgeR Am I simply missing something? Thanks -- Marco Blanchette,

Re: [R] If cycle takes to much time...

mmetric matrix, where raws and colums are artists and value I get 1 (or true) if they had an exhibition in the same and in the same place. My unelegant code is working, but for 96000 observation is requiring months and months. Gerrit is very elegant, but i run out of memory... the problem is t

[R] best ordination method for binary variables

are > 0, wascores) && any(comm < 0)) { : missing value where TRUE/FALSE needed* *In addition: Warning message:* *In Ops.factor(left, right) : < not meaningful for factor"* Do you think Principal coordinate analysis on a binary distance matrix is a decent strategy? Thanks

Re: [R] bnlearn: how to compute boot strength with mmhc and a blacklist

different network every time it is run on > the same data set? This is not surprising, because mmhc() does not have a "start" argument, so it's starting from the same network over and over. There is no way to provide a random seed to mmhc(), so the only way

Re: [R] Cross validation for Naive Bayes and Bayes Networks

to reassemble observed and predicted class labels and compute your metrics. > I also tried the *e1071* package, but I could not find a way to do > cross-validation. You might be able to trick the tune() function to do it, but I am not sure. Marco -- Marco Scutari, Ph.D. Research Assoc

[R] Goodness-of-fit tests for Complex Survey Logistic Regression

~dnum+snum, weights=~pw, data=apiclus2) a<-svyglm(sch.wide~1, design=dclus2,family=quasibinomial()) #Null model c<-svyglm(sch.wide~ell*meals, design=dclus2,family=quasibinomial()) anova(a,c) Many thanks for your help Marco [[alternative HTML version deleted]]

[R] colour coded dotchart

ealt with a similar task and found a solution? Thank you very much in advance for any help. Marco [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting

[R] [R-pkgs] lqmm package 1.02

available on https://r-forge.r-project.org/R/?group_id=1396 prior to submission to CRAN. -- Marco Geraci PhD CStat Lecturer in Biostatistics Institute of Child Health, University College London London WC1N 1EH, UK Web: http://marcogeraci.wordpress.com [[alternative HTML version deleted

Re: [R] Specifying a response variable in a Bayesian network

ne of the BN classifiers in bnlearn, naive.bayes()/tree.bayes(), which handle the concept of a response variable more naturally than general-purpose BNs. Hope it helps, Marco -- Marco Scutari, Ph.D. Research Associate, Genetics Institute (UGI) University C

[R] problem with the installation of r commander on a mac

Tcl/Tk and XQuark as recommended. I’ve also tried uninstalling R, reinstalling the latest version (2.15.1 signed) and all packages, and following the same official steps again. It didn’t work neither. Could you please help me? R Commander always worked fine on my iMac, until I updated to Mountain Lion.

[R] plotting principal components on geographic map

that the dataset includes relatively few points. Therefore, I would also like to interpolate the values for the uncovered geographical areas. Thanks in advance for any feedback, best marco [[alternative HTML version deleted]] __ R-help@r-projec

[R] recall variable in a loop

what I have to write in place of "y" in the last loop in order to refer to xx1, xx2 or xx3?? I want to write something like: for (i in 1:3) { xxi<-1} Thank you for any advice. Marco __ R-help@r-project.org mailing list https

Re: [R] NA/NaN values in bnlearn package R

owever, you can estimate the CPTs from incomplete data using table() and prop.table() and assemble them in a fitted BN with custom.fit(). On the other hand, maybe it would be better to write an EM wrapper around bn.fit() to make the best of the dependence structure of the data? Cheers,

[R] survplot invert number at risk labels

o solve the problem or to understand it. Thank you Marco. [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posti

Re: [R] survplot invert number at risk labels

my version), and this is why I put time.inc=1 in the previous call (This other issue is still present in 3.0.2) I'm going to upgrade my R and check that this little error is actually gone. Thank you Marco 2014-07-03 0:17 GMT+02:00 Frank Harrell : > Please try hard to cr

Re: [R] Median expected survival

Hi, Lars. I don't understand well your question. why don't you simply type pred_leuk Call: survfit(formula = leuk.cox, newdata = leuk_new) records n.max n.start events median 0.95LCL 0.95UCL 1 4242 42 30 1 1 23 2 4242 42 30 7

[R] fitting a gamma frailty model (coxph)

or in if (history[2, 3] < (history[1, 3] + 1)) theta <- mean(history[1:2, : missing value where TRUE/FALSE needed Does anyone have an idea on how to debug? Yours sincerely, Marco [[alternative HTML version deleted]] __ R-help@r-project.

[R] clustering of binary data

1 1 0 1 0 0 0 case12 0 0 0 1 1 0 1 1 0 1 0 0 0 . So, my questions are the following: Is the Jaccard index a good strategy for my kind of data? Is binary distance used in pvclust is theoretically more correct? Is there any alternative to pvclust for testing the significance of my clusters? Tha

[R] proportional weights

Hello all, can help clarify something? According to R's lm() doc: Non-NULL weights can be used to indicate that different observations have different variances (with the values in weights being inversely *proportional* to the variances); or equivalently, when the elements of weights are positiv

Re: [R] proportional weights

Thanks for the answers. Dear Marco and Goran, Perhaps the documentation could be clearer, but it is after all a brief help page. Using weights of 2 to lm() is *not* equivalent to entering the observation twice. The weights are variance weights, not case weights. According to your post here

Re: [R] proportional weights

I think we can blame Tim Hesterberg for the confusion: He writes " I'll add: * inverse-variance weights, where var(y for observation) = 1/weight (as opposed to just being inversely proportional to the weight) * " And, although I'm not a native English speaker, I think there's a spurious c

Re: [R] proportional weights

No, you are perfectly fine using WLS. The constant of proportionality is the estimated error variance, i.e., the square of the residual standard error (as I think I said earlier). John You're right. That was a little hard for me to grasp. Thanks for the patience.

[R] Exposures in bernoulli glm()

Hi, is it possible to add exposures to a glm with family=binomial()? It's easy to do it for a Poisson/negative binomial: just multiply the mean by the exposure, that is, offset(log(exposure)): but this obviously wrong for a binomial/Bernoulli since the mean must be no bigger than 1. My goal w

[R] [OT] Advice for medium size data management

me to bother at all. I'm not asking for prepackaged solutions, rather for help in documentation seeking and links to useful documentation or other threads (for example: is it worthwhile using parallel computing?) Thank you to anyone for reading this email. Marco Barbàra. P.S.: I work on a Debia

[R] Besag endive data

ataset, but I can no longer find it. Does anybody know where I can find the dataset? Thanks, Marco -- ***** Marco Bee University of Trento - Department of Economics and Management via Inama, 5 - 38122 Trento, Italy Phone: +39-0461-282296, Fa

Re: [R] Using "cpquery" function from bnlearn package inside loop

do so I suggest you should install the latest bugfix snapshot from bnlearn.com to avoid a few other bugs in cpquery(..., method = "lw"). Cheers, Marco -- Marco Scutari, Ph.D. Research Associate, Genetics Institute (UGI) University College London (UCL), United Kingdom __

Re: [R] Using "cpquery" function from bnlearn package inside loop

is easier to use in a programmatic way, but it is much more limited in the events you can use for conditioning. Logic sampling (method = "ls") accepts an arbitrary R expression as a conditioning event; likelihood weighting at the moment accepts only a single instantiation of the conditionin

Re: [R] Using "cpquery" function from bnlearn package inside loop

lihood weighting method will work. > Don't > hesitate to correct me if I'm wrong. >From your description, likelihood weighting should be fine. Cheers, Marco -- Marco Scutari, Ph.D. Research Associate, Genetics Institute (UGI) University College London (UCL), United Ki

[R] dynamic panelmodel pgmm

ot;,index=c("Aktie","Datum"),gmm.inst=~PB,lag.gmm=list(c(2,12)),transformation="ld") Fehler in FUN(X[[1L]], ...) : Indizierung außerhalb der Grenzen dim(panel) [1] 3408637 Best regards, Marco -- GRATIS für alle GMX-Mitglieder: Die maxdome Movie-FLAT! ___

[R] Adjustment for multiple-comparison for log-rank test

DeaR experts, I was asked for a log-rank pairwise survival comparison. I've a straightforward way to do this using the SAS system: http://support.sas.com/documentation/cdl/en/statug/63033/HTML/default/viewer.htm#/documentation/cdl/en/statug/63033/HTML/default/statug_lifetest_sect019.htm What

Re: [R] Bayesian Belief Networks

On Fri, Dec 24, 2010 at 12:01 AM, Michael Bedward wrote: > Hello Walt, > > Have a look at the bnlearn and deal packages. > > Michael Dear Walt, take a look also to the catnet/mugnet and pcalg packages, and if you have questions about bnlearn feel free to ask. Regards, Ma

Re: [R] help with a survplot

Thank you very much. Despite prof. Harrell's support (for whom I feel great esteem) I still remain doubtful about this feature. __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://w

Re: [R] help with a survplot

Il giorno Mon, 2 May 2011 07:30:10 -0700 (PDT) Frank Harrell ha scritto: > Please elaborate. It is simply that, generally speaking, i don't like adding numbers to a plot. I eventually realized that riskset cardinality may be a useful indication. However, do not discriminate between events and ce

[R] summing array elements

Hi to everybody I have an array with dimensions 2,4,3,3. Wanting to sum the matrices in the first two dimensions, I'm trying to use the "apply" function, but with no results. Have to say I'm quite new with R syntax. tx in advance marco [[alternative

[R] Why do not prevent users from creating a local copy of c() ?

doing such dangerous things? Thanks, Marco. __ R-help@r-project.org 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, re

Re: [R] summing array elements

Thanks a lot for your hint marco 2011/5/27 Erich Neuwirth > > arr<-1:72 > > dim(arr)<-c(2,4,3,3) > > apply(arr,1:2,sum) > [,1] [,2] [,3] [,4] > [1,] 297 315 333 351 > [2,] 306 324 342 360 > > apply(arr,3:4,sum) > [,1] [,2] [,3] > [

Re: [R] summing array elements

Thanks a lot for your help marco 2011/5/27 Thomas Lumley > On Fri, May 27, 2011 at 10:06 AM, marco milella wrote: > > Hi to everybody > > I have an array with dimensions 2,4,3,3. > > Wanting to sum the matrices in the first two dimensions, I'm trying to > use

Re: [R] about chisq.test()

> > Is this result reliable? No. chisq.test is a non parametric goodness-of-fit test for *counts*. In my opinion you need a one sample t.test. __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting

Re: [R] How to plot the sample size on top of the bars of a barchart plot?

Il giorno Wed, 19 May 2010 08:46:06 -0700 (PDT) Javier ha scritto: > > Dear users, > > Is there anyone how knows how to plot the sample size on top of the > bars of a barchart plot? Do you mean something like this? mp <- barplot( VADeaths[1,], ylim = c( 0 , 18 ) ) text( mp , VADeaths[1,] + 0

Re: [R] [Off topic?] Time dependent Cox model fitting and validation

this is self-reply only to insert my real name __ R-help@r-project.org 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,

Re: [R] Time dependent Cox model

eed time interaction... > > The cox.zph function is primarily graphical; I would respond to your > question with "is it good to look at scatterplots before fitting a > linear model?" My answer to this is emphatically yes. > > Terry Therneau > > Thank y

[R] Modeling time varying effects in with cph: how to ?

t this is generally not recommended, isn't it? Any sort of comment will be appreciated. Thank you very much. Marco Barbàra. __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide

[R] Partial correlation analysis for censored data in R?

Hi does anyone know if there is any R function that preforms partial correlations analysis for censored data? I found a fortran program on the PSU website at http://www2.astro.psu.edu/statcodes/cens_tau.f but I was wondering if there is anything like that in any R package. Thanks m. _

  1   2   >