similar to: Tukey and extracting letters in multcomp

I wish to estimate sparse causal networks from simulated time series data. Although there's some discussion about this problem in the literature (at least a few authors have used lasso and l(1,2) regularization to enforce sparsity in multivariate autoregressive models, e.g., http://user.cs.tu-berlin.de/~nkraemer/papers/grplasso_causality.pdf), I can't find any R packages with these

I apologize in advance if this is posted all ready I have not been able to find any information about it. I have this data frame and I want to sort smoking by retlevel. Age Gender BMI Calories Fat Fiber Alc retlevel Smoking 1 64 Female 18.87834 1828.0 63.4 14.7 0.0 Normal Non-Smoker 2 25 Female 20.64102 1517.4 59.1 5.9 0.0 Normal Smoker 3

I dont get it: > is.vector(c(mydf[1])) [1] TRUE > unique(c(mydf[1])) Error in unique(c(mydf[1])) : unique() applies only to vectors > Is it a vector or not? This stuff is driving me nuts. I'm simply trying to convince R that my grouping vector is actually a vector so that unique will work. Its just a vector of numbers, so why shouldnt it work? --

I was wondering if there is a way of predicting factor scores of new data for factor analysis in R (similar to "predict" in S-plus). So far I have not been able to find it, nor found reference to it. Thanks, Chris Fonnesbeck -- ~*~*~*~*~*~*~*~*~*~*~*~*~*~*~*~*~*~* Christopher J. Fonnesbeck Ph.D. Student Georgia Cooperative Wildlife Unit University of Georgia Athens, GA 30602

Can one perform a QDA in R? I do not see it anywhere within the mda package. Any pointers here would be appreciated. Thanks, cjf -- ~*~*~*~*~*~*~*~*~*~*~*~*~*~*~*~*~*~* Christopher J. Fonnesbeck Ph.D. Student Georgia Cooperative Wildlife Unit University of Georgia Athens, GA 30602 Email: cjf at fonnesbeck.net Yahoo: fonnesbeck_chris ~*~*~*~*~*~*~*~*~*~*~*~*~*~*~*~*~*~*

I know this seems like a very easy question (and maye it is) but I've been trying to simulate nested data and been unsucessful so far.. I want to simulate a varying number of species within a group; and then create an array to store the results of my for-loop. For example: groups<-3 species$groups<-as.integer(runif(groups,1,10)) #species per functional group ###create arrays to store

I use mac os 10.3.9 and I've installed in my computer R 2.1.0 (I believe this is the latest R version). Although it works alright when I open R by clicking in the R icon, I cannot use the graphics facility when I open R from a X-terminal window (or x-11 window). In fact, when I try to open R I get the message that I pasted below. Is this a R bug? Is R assuming that I should have the

Dear list, I'm having trouble to see how my elements on a 4 dimensional array are listed. For example, I generated the following array: junk.melt=melt(occ.data,id.var=c("Especie", "Site", "Rep", "Año"), measure.var="Pres") y=cast(junk.melt, Site ~ Rep ~ Especie ~ Año) Now, I want to be able to look at how my species (Especie) are listed, in

Hi, Is there a routine available in R for the Newton-Raphson method for simulataneous equations in several unknowns? Thanks Robert -- Robert J. Chandran Department of Botany 3506 Miller Plant Sciences Building University of Georgia Athens, GA 30602 Phone: (706)-583-0943 -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read

Hi, I have an samba sevrver as PDC of an NT domain. I want to join another machine(samba) to the domain as a member of NT domain. How can I do it. Thanks -------------------------------- Jianping Zhu Department of Computer Science Univerity of Georgia Athens, GA 30602 Tel 706 5423900 --------------------------------

Perhaps I am missing something incredibly basic, but I can't get the wireframe function to plot anything other than an empty graph with axes (and the correct scale). Cloud() works just fine, as does scatterplot3d(), and wireframe() works fine for other data sets I have. Are there some data sets that can't be plotted with wireframe()? If so, what are the characteristics that prevent

Hi, I'm trying to perform a glm model on count data (poisson distribution of the errors) where data are nested. glmmodel<-glm(y~x/z,poisson) x and z are factors, z nested within x, y is count data. In that point the R just "stuck" and not respond anymore. I tried glmmodel<-glm(y~x,poisson) and there were no problems. Trying glmmodel<-glm(y~x/z) gave the same no response.

I am trying to do Tukey HSD comparisons on a repeated measures expt. I found the following example on r-help and quoted approvingly elsewhere. It is broken. Can anyone please tell me how to get it to work? I am using R 2.4.1. > require(MASS) ## for oats data set > require(nlme) ## for lme() > require(multcomp) ## for multiple comparison stuff > Aov.mod <- aov(Y ~ N + V +

R version 3.2.2. library(vegan) I was to look at community tables from my dendrograms and am trying out the vegemite command. This is the error I get: Error in vegemite(apst, apst.clusters) : Cowardly refusing to use longer than 1 char symbols: Use scale I thought the problem was that I was using the log transformed data, but I tried it on the raw (which is single digit numbers), and still

i know with icecast2/ices , karl heyes has some example config files for mixing live/static content ( if i recall correctly ) . however, what i had done in the past is add the static file to the .pls prior to the live mount . if you generate the .pls files on demand , this would work . ALSO , i'm quite sure that both icecast1/2 support intro/outtro files , etc .. this may achieve the same

I've been running classic icecast + liveice for some time now (several years) without problems. A couple days ago, the station manager (NPR affiliate) came to me and asked if we could 'inject' a static file into the stream. Basically what he wants is for the client to first receive a static file (requesting donations of course) and then on to the live stream (*not* intermittently).

Hi all, I would like to know whether there is any function in R were i can find the cross-correlation of two or more multivariate (time series) data. I tried the function ccf() but it seems like to have two univariate datasets. Please let me know. sincerely, sandeep -- Sandeep Joseph PhD Post Doctoral Associate Center for Tropical & Emerging Global Diseases Paul D. Coverdell Center,

HI, I''m getting some strange errors when I try to use a variable in a require => Package construct, namely: $packagelist = ["clamav", "clamd", "clamav-db"] package { $packagelist: ensure => latest, } # Above apparently works fine... # ...THEN... file { "/etc/rc.d/init.d/clamd.mimedefang": ensure =>

Hi, I have a few questions about glht() and the interpretation of output from Tukey's in multcomp package for lme() model. The main issue is that I noticed that a plot that I produced with code letters seem to contradict the graph itself. I provide data and code below. I end with my questions. A few things about data set. "LMA.vcp" is continuous response variable.

Hi, I am running R on my iBook (OS X 10.4) and whenever I try and plot a model that has been run, I get a grey background. This makes the plot difficult to read on b&w printouts. I have tried changing the background to white (or any other color) without any success. Here's a current example plot(augPred(m2corn.nlme, level=0:1), bg="transparent", lty=c(1,4), pch=19,