Re: [R] How to import an excel data file

Dear David,  I also use excel files for R. This is what I do install.packages("readxl") library(readxl) mydata=read_excel('mydata.xlsx') Hope this helps.  Best, Maria Στις Σάββατο 23 Σεπτεμβρίου 2023 στις 09:23:43 π.μ. GMT+1, ο χρήστης Parkhurst, David έγραψε: I know I s

Re: [R] how to change the y-axis to logarithmic in a barplot ggplot

Hi Ivan,  yes you are right. I have many values close to zero and if I use log, they are instantly transformed to negative values and then the bars show above. The sqrt trans worked ok but when I used the pseudo_log_trans, it did the work! Many thanks.  Kind regards, Maria Στις Δευτέρα

Re: [R] how to change the y-axis to logarithmic in a barplot ggplot

Dear Bert,  Thank you very much for the tip. The reason that I took the initiative to write here is because I receive here as well questions about ggplot. Next time I will use the ggplot help.  Maria Στις Δευτέρα 17 Ιουλίου 2023 στις 12:20:25 π.μ. GMT+1, ο χρήστης Bert Gunter έγραψε

Re: [R] how to change the y-axis to logarithmic in a barplot ggplot

I will find the ggplot help.  But I have tried everything, including what you have suggested and nothing works. Kind regards, Maria Στις Κυριακή 16 Ιουλίου 2023 στις 11:22:36 μ.μ. GMT+1, ο χρήστης CALUM POLWART έγραψε: Try adding scale_y_log10() This is a general R help list

[R] how to change the y-axis to logarithmic in a barplot ggplot

Dark2") + scale_color_brewer(palette = "Dark2") I would very much appreciate your help. I am stuck.  Kind regards, Maria __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PL

Re: [R] [External] Error in percentage stacked barplot

= FALSE strip.right = TRUE likert(t(hellisheidi), ReferenceZero=.5, xlab="X-lab", ylab="Y-lab", main="Stacked bar chart", col=c("#E94E1B", "#F7AA4E", "#BEBEBE", "#6193CE", "#00508C", "#E94E1B", "#

[R] Error in percentage stacked barplot

omponent, instead of being the variable in the Samples. I was hoping if you could help me on that. Thank you very much in advance. Kind regards,Maria __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/l

[R] Unsubscribe from mailing list please

Hi, Could I please unsubscribe from the R-help mailing list? Thanks Eva [[alternative HTML version deleted]] __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the

[R] different dimensions in W and my data in moran and spatial model

w = listw1), R= 500, zstats = TRUE,zero.policy=TRUE) I don't know if it is for my panel data, the missing values or something else. I tried also in Stata having similar issues. thanks in advance and my apologies if it is very naive question, but I am stuck. all the best Maria Jesus

[R] Question about the package "MatchIt"

Hi! I'm trying to perform propensity score matching on survey data and so for each individual observation I have a statistical weight attached. My question is: is there a way within the package to consider these weights in the matching procedure? Thank you very much. -- Maria Cristina Maurizio

[R-es] OLvide contraseña tenai este problerma que queria publicar

estoy haciendo este ejercicio y en Crear una variable de rango basada en la importancia rankImportance% mutate (Rank = paste0 ('#', dense_rank (desc (Importance Me da este error Error en varImportance%>% mutate (Rank = paste0 ("#", dense_rank (desc (Importance: no se pudo encontrar la

[R-es] Olvide contraseña

e-mail: brew5...@gmail.com -- *María * [[alternative HTML version deleted]] ___ R-help-es mailing list R-help-es@r-project.org https://stat.ethz.ch/mailman/listinfo/r-help-es

Re: [R-es] Ayuda para campo con numerico y texto con Problemas

Por parte No puedo usar otra cosa que sea R, por motivos varios, Uno de ellos es tporque quiero hacer un script para introducir cada provincia conforme llegue, y es son muchos datos para hacerlo con excel, y tambnie que si lo hago mas año seria muy pesado. Primero que esta pasando que me

[R-es] Como se que he mandado algo correctamente a la lista

Hola Como se que lo que he mandado ha llegado correctamente a la lista y cuanto es la media de obtener ayuda. M- [[alternative HTML version deleted]] ___ R-help-es mailing list R-help-es@r-project.org

[R-es] Ayuda para campo con numerico y texto con Problemas

Hola a todos, Soy nueva usando R y tengo unos problema en procesado de datos que quiero hacer en R. Tengo un fichero de datos en formato texto separado con ; (punto y coma). Este fichero contiene las columnas: Poblacion;Sexo;Nota; EdadP, EdadM Tengo 2 problemas que no se como tratar. 1.-

Re: [R] External validation for a hurdle model (pscl)

tackexchange.com > >not here; I am pretty sure you'll receive lots of answers there. > > > >Cheers, > >Bert > > > > > >Bert Gunter > > > >"The trouble with having an open mind is that people keep coming along > >and > >sticking t

[R] External validation for a hurdle model (pscl)

Hi R-list, We have constructed a hurdle model some time ago. Now we were able to gather new data in the same city (38 new sites), and want to do an external validation to see if the model still performs ok. All the books and lectures I have read say its the best validation option but... I have

[R] pROC

Hello, I'm using the package pROC in RStudio to create my ROC-curves and I have patients in my data - healthy or sick - in the column "Status" and the value "SUVmax" to examine it. I used the following code:

[R] Σχετ: x-axis tick marks length in ggplot2

Dear Walter,  I tried to use scale_x_continuous but the arguments that I found was to change the labels, the limits and the breaks. I was only able to increase the number of the tick marks.  Best,Maria Στις 3:14 μ.μ. Πέμπτη, 28 Ιουνίου 2018, ο/η Walter Pina έγραψε: Dear Maria, you

[R] Σχετ: geom_text only in the first panel with facet_wrap in ggplot2

Dear Ulrik and all, Thank you all so much; well, Ulrik's suggestion worked better. Very much appreciated. Best,Maria Στις 11:17 π.μ. Δευτέρα, 25 Ιουνίου 2018, ο/η Ulrik Stervbo έγραψε: Hi Maria, you are on the right way. The data.frame with the text must have the same columns as you

[R] geom_text only in the first panel with facet_wrap in ggplot2

t an errorann_text <- data.frame(Q = 20, fit = 1.03, lab = "Text", asb = factor("ASB1",levels = c("1","2","3"))) #Error in FUN(X[[i]], ...) : object 'ASB' not found I would very much appreciate for your help. Thank you very much in advance. K

[R] interaction term by a factor group in gamm4

Dear all, I am writing as I would really need your help on the problem with gamm4. I have tried to find a solution online but I wasn't very successful. I am running a gamm4 model with an interaction between two variable using the tensor term, t2. I have a group variable (super end group) with

[R] svyglm

consdentcat2", "sexo", "estrato", "upa", "pesomorcc"), row.names = c(NA, 10L), class = "data.frame")## bteste1 <-bteste[1:10, ]#bteste1   consdentcat2 sexo estrato     upa    pesomorcc1             

[R] Σχετ: Σχετ: change colour in barplot

Dear Duncan Many thanks for this; actually the vector solution works much better.  Best, Maria Στις 12:27 μ.μ. Δευτέρα, 20 Νοεμβρίου 2017, ο/η Duncan Murdoch <murdoch.dun...@gmail.com> έγραψε: On 20/11/2017 7:09 AM, Maria Lathouri wrote: > Dear all, > > I am really so

[R] Σχετ: change colour in barplot

Dear all,  I am really sorry for this. I have attached the script and a .csv file with an example.  Hope this will help. Many thanks,Maria Στις 11:53 π.μ. Δευτέρα, 20 Νοεμβρίου 2017, ο/η Duncan Murdoch <murdoch.dun...@gmail.com> έγραψε: On 20/11/2017 6:38 AM, Maria Lathouri via

[R] change colour in barplot

colour changes in the legend but not the actual bars in the plot. I don't know what I am doing wrong. I know that I am missing something but I cannot figure it out.  I would very much appreciate for your help.  Many thanks,Maria [[alternative HTML version deleted]] ___

[R] Σχετ: show 0 at y axis in xyplot lattice

it doesn't go through.  I will know for the next time.  Many thanks. Regards,Maria Στις 2:11 π.μ. Τρίτη, 12 Σεπτεμβρίου 2017, ο/η Duncan Mackay <dulca...@bigpond.com> έγραψε: Hi Maria Rule 1 make sure your data is in the right format dat <- source("G:/1/savedat.txt") > d

[R] Σχετ: show 0 at y axis in xyplot lattice

"black", "red"), fill=c("red", "black", "red"),    key=dat_key, scales = list(x = list(at = sdate, labels = format(sdate, "%b-%y")), y=list(at=c(0, 20, 40, 60, 80))), xlab="Date",    panel = function(x, y, ...) { panel

[R] show 0 at y axis in xyplot lattice

Dear all I am trying to make a plot in xyplot lattice by groups. I would like to show "0" in y axis but I don't want to be aligned with the x axis. I want to be a little bit above. I have tried many options but I don't get what I want. I also tried ylim=c(0, 80) but both 0 and 80 are fully

[R] Σχετ: plotting gamm results in lattice

and the confidence intervals.  Kind regards,Maria Στις 4:26 μ.μ. Δευτέρα, 12 Ιουνίου 2017, ο/η Duncan Mackay <dulca...@bigpond.com> έγραψε: Hi Maria If you have problems just start with a small model with predictions and then plot with xyplot the same applies to xyplot Try library(gamm4)

Re: [R-es] Recta de regresión en un Scatterplot

int...@gmail.com>: 2017-03-13 13:39 GMT-03:00 Maria Rubio <maria.rub...@um.es>: > model=lm(d_curts ~ d_asims + I(d_asims^2)) > plot(d_asims,d_curts,xlab="Skewness", ylab="Kurtosis") ​Podrías añadir: abline(model)​ ​Seguramente te arroje un warning porque ti

[R] Negative binomial GAMM using 'by' in factor interactions

I am using a GAMM to model my data (this is as far as I know the only way I can use the negative binomial distribution AND a correlation structure within the model). I measured animal detections (including zero detections) per hour at 3 different locations in an area. location is a factor in my

[R] Σχετ: Σχετ: Σχετ: plotting dates in x axis

Hi Petr I followed your suggestion and I have attached a saved txt file with the data and the script using dput(). I hope this time will work.  But by the way your code worked as well.  Many thanks.  Kind regards,Maria Στις 1:06 μ.μ. Τρίτη, 28 Φεβρουαρίου 2017, ο/η PIKAL Petr <petr

[R] Σχετ: Σχετ: plotting dates in x axis

if I can show every two years instead. I think it would be much better.   Many thanks.  Kind regards,Maria Στις 11:56 π.μ. Τρίτη, 28 Φεβρουαρίου 2017, ο/η Maria Lathouri via R-help <r-help@r-project.org> έγραψε: Dear all, Sorry about that. My mistake. Here is an example of my data

[R] Σχετ: plotting dates in x axis

; 2012-12-10 0.484 I hope this helps. Maria Στις 10:00 μ.μ. Δευτέρα, 27 Φεβρουαρίου 2017, ο/η Jeff Newmiller <jdnew...@dcn.davis.ca.us> έγραψε: While humorous, the term "Mexican Wall" is unlikely to be clear to the OP. This is a reference to the mailing list anti-v

[R] concurvity

Hi, I'm using the concurvity function to check for concurvity in my model. The output I get when comparing to the rest of the model (ie full=TRUE) many of the variables have concurvity values higher than 0.9. However when comparing the terms pairwise most values are very small, less than 0.1

[R] plotting dates in x axis

ot(sdate, var1, type="l", xaxt="n")> axis(side=1, at=0:15, > labels=strftime(d2, format="%m-%Y"), cex.axis=0.8,las=2) But nothing happened. I cannot understand why it doesn't show anything.  I have attached the file as well in case you want to have a more cl

[R] Theta in Negative binomial GAM

Hi, I'm fitting a negative binomial GAM (using mgcv) to my data using family=nb() so theta is estimated during the fitting process. When I then extract this theta from the model and refit the same model with family=negbin(theta) it gives a much lower AIC. I know using AIC to compare negative

[R] scatter plot of numerical variables against different sample ids

ng.  I would very much appreciate if you could help me on that. Apparently, I am missing something.  Thank you in advance.  Kind regards,Maria [[alternative HTML version deleted]] __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more,

[R] change the colour line in gamm4 plotting

Dear all I am stuck probably in a simple plotting question My model is:model.1<-gamm4(y~s(x1, by=end.group)+Year+K1+k2+k3, data=.., random=~(1|WB_ID/Site_ID)) where y is my dependent variable, x1 is the smooth covariate and I use the by argument for the smooth term based on six different

[R] R version 3.2.5

Good night, How can I download and install R version 3.2.5 from cran mirror of University of São Paulo, São Paulo? Thanks for your support. Maria Alice P. Jacques [[alternative HTML version deleted]] __ R-help@r-project.org mailing list

[R] ZINB multi-level model using MCMCglmm

Hi, I am Olga Viedma. I am running a Zero-inflated negative binomial (ZINB) multi-level model using MCMCglmm package. I have a doubt. Can I use the "Liab" outputs as fitted data, instead of the predicted values from "predict"? The liab outputs fit very well with the observed data, whereas

[R] sorting of files using system

st; I was surprised by this result and not sure if it's expected? I am using R version 3.2.2 (2015-08-14) on Ubuntu 15.04 Thank you in advance, Maria __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listi

[R] Use of R for estimation of Discrete Choice Models incorporating latent variables

Dear all, I would like to know if it is possible to estimate an integrated choice and latent variable model with R. If it is so, please send me the name of the respective package. Best regards, Maria Alice P. Jacques Associated Researcher University of Brasilia - Brazil [[alternative

[R] Glmnet penalty.factor with multigaussian response

NA hsa-let-7c -0.06436899 -0.03362183 -0.007440406 -0.2581606 -0.01517728> dim(fullcoefs)[1] 723 7 More speciffically, in this case I would expect some value for "hsa-let-7a":"CYP1A2", as the penalty for it is 0. Many thanks! Maria [[alternative HTML version deleted]

[R] Multiple if function

Dear all, I am writing as I would like your help. I have a dataframe with two columns, ASB and Flow, where the the first one has values 1, 2 or 3 and the second flow data. Something like that: ASBclass    Flow1  11.51   9.2 2  10.5 3   6.7  ...   

Re: [R] clm funtion and CI

Data: Quinta-feira, 9 de Julho de 2015, 11:44 You need a reproducible example. On Wed, Jul 8, 2015 at 7:43 PM, Luciane Maria Pilotto wrote: Hi, I'm working with ordinal logistic regression and fitting the model with the clm funtion of the ordinal package and would like to get the CI

[R] clm funtion and CI

Hi, I'm working with ordinal logistic regression and fitting the model with the clm funtion of the ordinal package and would like to get the CI. According to the Tutorial on fitting Cumulative Link Models with the ordinal Package, Rune Haubo B Christensen (21 January 2015) you can run the OR,

[R] Output interpretation: standard error of lm dummy variable

Hi guys I have a statistical question to an analyse I ran in R. It is a dummy variable model with the 5 regions of Denmark as 4 independent dummy variables and price as the dependent variable: price = 10.325 - 0.176*Sjaeland - 0.368 * NJylland - 0.230*MJylland - 0.120* Syddanmark I

Re: [R] Urgent :) Output interpretation: standard error of lm dummy variable

and you can help :) Best Livia Fra: Livia Maria Vestergaard Sendt: 19. maj 2015 10:16 Til: r-help@r-project.org Emne: [R] Output interpretation: standard error of lm dummy variable Hi guys I have a statistical question to an analyse I ran in R

[R] lm model exported from R to excel

Hi all I all. I am wondering whether anybody know how to export an output of an lm model from R to excel in order to have excel recognize the table that comes and divide the numbers in the table into columns and rows? I really hope it is possible? :) Best Livia

Re: [R] lm model exported from R to excel

Hi Duncan Thank you so much - it worked :) Best Livia Fra: Duncan Mackay [dulca...@bigpond.com] Sendt: 6. maj 2015 14:26 Til: R; Livia Maria Vestergaard Emne: RE: [R] lm model exported from R to excel Hi Livia There are several html packages

Re: [R] help - hoslem.test

Maria Pilotto     r-help@r-project.org   r-help@r-   project.org     Data: Quinta-feira, 30 de Abril de 2015, 11:51       Kevin Thorpe   pointed out     to me that there is a   dropbox link at the very bottom of the     post that I missed. :(       I     just downloaded it, read

Re: [R] help - hoslem.test

=id3) tp1 library(ResourceSelection) hoslem.test(tp1$q131, fitted(tp1), g=10) dataframe: https://www.dropbox.com/s/9qrdf4mhd6tzypi/id3.rda?dl=0 __ Luciane Maria Pilotto Mestre e Doutoranda em Saúde Bucal Coletiva - FO/UFRGS  NDE Odontologia

Re: [R] LM() and time in R

() , but not with the time, how can I do that? Or do you guys maybe a third idea about what should make a regression about to improve the model? Cheers Livia :) Fra: Boris Steipe [boris.ste...@utoronto.ca] Sendt: 29. april 2015 20:25 Til: Livia Maria

[R] LM() and time in R

Hello, I need some help with a project that I’m working one. Im trying to make a l regression model (lm) in r with time as independant variable and gas prices as the depended. But It seems like everything im trying to run it, R freeze, I think that I need to tell R somehow that my time is

[R] help - hoslem.test

Hello, I'm working with ordinal logistic regression model (polr) and would like to test the proportional odds assumption. For this, I ran the binary logistic regressions with varying cutpoints on the dependent variable, as described in the following commands. When running the test of Hosmer

Re: [R] error MANOVA in R

Hi Peter, Thank you so much for your tips. Have a nice Easter. Gian Il giovedì 2 aprile 2015, peter dalgaard pda...@gmail.com ha scritto: On 30 Mar 2015, at 17:11 , Gian Maria Niccolò Benucci gian.benu...@gmail.com javascript:; wrote: Dear R-usrs, I am trying to perform a MANOVA

[R] error MANOVA in R

Dear R-usrs, I am trying to perform a MANOVA on a data frame with 31 columns about soil parameters and 1 column containing the explanatory variable (Fraction) that have three levels. my code is the following: datam - read.table(data_manova2.csv, header=T, sep=,) names(datam) manova_fraction2 -

[R] Fitting distributions in uneven parameters

Hi All, First of all, I would like to wish you all a Happy New Year, full of creativity, inspiration and prosperity. Now, I have a data set with two uneven size parameters  which is the following:     x1    x2 1  0.98    0.952  0.99    0.98 3  1.11    1.01

[R] Cox model with multiple events - Proportional Hazards Assumption

Dear all, I'm using the package survival for adjusting the Cox model with multiple events (Prentice, Williams and Peterson Model). I have several covariates, some of them are time-dependent. I'm using the functioncox.zph to check the proportional hazards. Due to the nature of the time-dependent

[R] taking daily modes from hourly data

Dear all, I have a data frame (datos) of hourly wind speed and direction with 4columns (1st date, 2nd hour, 3rd wind speed and 4rth wind direction). I have been able to do the daily mean of the wind speed, but when I try to get the more frequent wind direction of every day, it doesn't work. I

Re: [R] Cansisc: Error in eigen(eHe, symmetric = TRUE)

Of Michael Friendly Sent: Thursday, July 10, 2014 3:22 AM To: Maria Judith Carmona H Cc: r-help@r-project.org; John Fox Subject: Re: [R] Cansisc: Error in eigen(eHe, symmetric = TRUE) Maria The variables Araceae,Begoniaceae,Bromeliaceae,Clusiaceae,Cyclanthaceae,Ericaceae,Gesneriaceae

[R] Cansisc: Error in eigen(eHe, symmetric = TRUE)

) : infinite or missing values in 'x' In addition: Warning message: In sqrt(wmd) : NaNs produced Please help! -- Maria Judith Carmona Higuita. Estudiante de Biología - Universidad de Antioquia Medellín - Colombia La felicidad ocurre cuando encajas en tu vida, cuando encajas tan armónicamente que

Re: [R] Cansisc: Error in eigen(eHe, symmetric = TRUE)

produced You see I get the same error. Best regards, Judith On Wed, Jul 9, 2014 at 5:30 PM, John Fox j...@mcmaster.ca wrote: Dear Maria Judith Carmona Higuita, Since you didn't include enough information (such as your access to your data) to reproduce the error, one can only guess. My guess

Re: [R] Cansisc: Error in eigen(eHe, symmetric = TRUE)

again for several days. Best, John On Wed, 9 Jul 2014 17:54:56 -0500 Maria Judith Carmona H juditycarm...@gmail.com wrote: Dear John, I am including abundance values ??in my data set so obviously I have zero abundances. The problem is that if plot only the factors (biomasa, altdosel

[R] grouping explanatory variables into sets for GLMM

my model with these sets only instead of listing every variable. Also: is there a way of running all possible models with the different combinations of these sets and/or variables, sort of like running ordistep for ordinations? Thanks a bunch in advance for your help! Maria

Re: [R] FrailtyHL package does not work

Therneau De: Maria Helena Mourino Silva Nunes Enviado: sexta-feira, 14 de Março de 2014 19:02 Para: r-help@r-project.org Assunto: FrailtyHL package does not work Dear all, I want to describe the time interval between mammography screenings as a function

[R] FrailtyHL package does not work

Dear all, I want to describe the time interval between mammography screenings as a function of several clinical and demographic variables. For each woman (ID), I have multiple events and each event corresponds to the time interval between two consecutive mammography screenings. Also, women's

[R] Standardised Pearson residuals

Dear all, I'm using the package pscl for adjusting a Zero-Inflated Negative Binomial Regression to my data set. I would like to know if it is possible to compute the standardised Pearson residuals from the output of this package. Thanks in advanced. Best regards, Helena.

[R] merge multi part polygons

? I tried this with unionSpatialPolygons but another error appears: Error in unionSpatialPolygons (sp1, IDS, threshold = NULL,: input lengths differ Can anyone help me? Att, Maria Fernanda BonettiDoutoranda do Programa de Pós Graduação em Ecologia e Conservação - UFPRMestre em Ecologia e

[R] How do I get a a p-value for the output of an lme model with lme4?

Hi there, I just started using lme4 and I have a question about obtaining p-values. I'm trying to get p-values for the output of a linear mixed-effects model. In my experiment I have a 2 by 2 within subjects design, fully crossing two factors, Gram and Number. This is the command I used to

[R] Make cdf package error Human Exon array

or if there is a better method for making the package, so any help will be really welcome. Thank you, Maria -- The University of Edinburgh is a charitable body, registered in Scotland, with registration number SC005336. __ R-help@r-project.org mailing

Re: [R] column and line graphs in R

Thank you very much to you all, I'll play the codes and post my code once I have tested it. Cheers, -- Gian On 14 March 2013 16:27, John Kane jrkrid...@inbox.com wrote: The easiest way to supply data is to use the dput() function. Example with your file named testfile:

[R] how to import DNA sequences in R

Dear all, I was wondering how to read DNA sequences in R, is there a specific function and/or a specific package for that? Thank you very much in advance, Gian ** [[alternative HTML version deleted]] __ R-help@r-project.org mailing list

[R] column and line graphs in R

Hi all, I would love to plot my data with R. I have abundance and frequency of fungal taxonomic data that should be plotted in the same graph. In Microsoft Excel is that possible but the graphic result is, as always, very poor. Is there a function that may let me plot these data in R? I have a

Re: [R] how to import DNA sequences in R

Thank you all guys, for the useful information. Gian On 14 March 2013 14:14, Martin Morgan mtmor...@fhcrc.org wrote: On 03/14/2013 02:19 AM, Gian Maria Niccolò Benucci wrote: Dear all, I was wondering how to read DNA sequences in R, is there a specific function and/or a specific package

Re: [R] column and line graphs in R

Hi again, Thank you all for your support. I would love to have a graph in which two variables are contemporary showed. For example a histogram and a curve should be the perfect choice. I tried to use twoord.plot() but I am not sure I understand how to manage the the arguments lx, ly, rx, ry...

[R] Installation of rJava package fails

Hi, I am trying to install the rJava package in a Linux environment (Platform: x86_64-suse-linux-gnu (64-bit)), where I am facing an issue which I cannot solve so far and where I could require some help: The installation seems to go fine, however in the last step of the installation where

[R] dealing with different kinds of environmental data in one matrix

Hi, I have done a CCA with an environmental data matrix including both relative abundance of plant species, soil properties (continuous, with different units of concentrations), and factors that should explain my species matrix. Only after running the analysis have I realized that I should have

[R] CCA plot

Hi, I have a couple questions about fitting environmental (land use factors, plant species presence-absence, and soil variables) constraints to my CCA biplot. 1. After successfully plotting species and site scores in my CCA, I have been trying to insert the biplot arrows of the environmental

Re: [R] polycor package

/ On Fri, 16 Nov 2012 12:58:03 + Laura Maria Schwirz schwi...@tcd.ie wrote: I am currently working with R's polycor package and I have encountered a problem. I tried to follow the steps as outlined in the sem.pdf file where a CFA model is run using polychoric correlations. Every time I

[R] polycor package

I am currently working with R's polycor package and I have encountered a problem. I tried to follow the steps as outlined in the sem.pdf file where a CFA model is run using polychoric correlations. Every time I run the command sem(model, data, N=.), I get the following warning message: Warning

Re: [R] find the Best-ticker

Thank you On Mon, Oct 29, 2012 at 1:44 PM, R. Michael Weylandt michael.weyla...@gmail.com michael.weyla...@gmail.com wrote: On Oct 29, 2012, at 6:11 AM, sheenmaria sheenmar...@gmail.com wrote: i need to find the best ticker from the group of some tickers.? i also need to know on what

[R] Vars package: impulse response function

Hello, I'm using VAR models in R in order to obtain impulse responses of stock market shock on US economy. I have series of quarterly changes in real gdp, SP 500 and quarterly level of unemployment for 1985 - 2012 period. My series are stationary. So I did all the steps below. However I don't

[R] [Vars] package: impulse response function

Hello, I'm doing a research on the impulse responses in VAR models and I'm having troubles in interpretation of R results. My question is what is the shock of impulse variable that is produced to obtain the response? Is it one-standard-deviation positive shock? If it is so how can I obtain the

Re: [R] random sampling inside a dataset

= FALSE' Why it does not work? The whole dataset is composed of 536 rows and I just want to sample randomly 200 of them... Thank you in advance, Gian On 13 September 2012 14:01, Duncan Murdoch murdoch.dun...@gmail.com wrote: On 12-09-13 7:43 AM, Gian Maria Niccolò Benucci wrote: Thank

[R] Explore SPSS function in R

Hi everyone, Does anybody knows if there is an equivalent R function that gives the same outcome as in Explore function in SPSS ? (Analize-Descriptive Statistics-Explore) It does a categorical vs quantitative variables analysis. ( But not linear regression) I need to compare intragroup

[R] RV: problems for making grids from lmer models

De: MARIA OLGA VIEDMA SILLERO Enviado el: mié 26/09/2012 13:09 Para: ba...@stat.wisc.edu; bbol...@gmail.com; maech...@r-project.org Asunto: problems making grids from lmer models [[alternative HTML version deleted

[R] random sampling inside a dataset

Hi all, I am wondering if do exist a function in R that allow me to sample or choose randomly the rows (i.e., samples) inside a given matrix. Thank you very much in advance. Cheers, -- Gian [[alternative HTML version deleted]] __

Re: [R] random sampling inside a dataset

of Statistics University of Bologna, Italy via Belle Arti 41, 40124 BO. -- *Da:* Gian Maria Niccolò Benucci gian.benu...@gmail.com *A:* r-help@r-project.org *Inviato:* Giovedì 13 Settembre 2012 10:42 *Oggetto:* [R] random sampling inside a dataset Hi all, I am wondering

[R] anova for two mixture model objects

are not nested” I wonder if anyone can help Thanks in advance Maria Maria Manuel Angélico IPMA, Portuguese Institute for Sea and Atmosphere Lisboa , Portugal tel: ++ 351 213027068 angel...@ipimar.pt [[alternative HTML version deleted

Re: [R] cluster with mahalanobis distance

How can I perform cluster analysis using the mahalanobis distance instead of the euclidean distance? Thank you Maria Froes [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r

[R] Error on easy way for JoSAE Package

Dear Mr/Mrs, I have got the the data sets in the attachement (txt). The fact is that on applying the easy method for JoSAE package on these data sets I got the following error: result - eblup.mse.f.wrap(domain.data = d.data, lme.obj = fit.lme1) Error in solve.default(asympt.var.covar) :

[R] $ operator is invalid for atomic vectors

, Ana Maria Dobre [[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/posting-guide.html and provide commented, minimal

[R] help - history()

Hi Members, Do exist the possibility to delete a command line into the history? Fo example. If I' ve typed a line code that is wrong, Can I delete it from the hostory in reason to do not save it in the .Rhistory file? Thanks for helping, Gian [[alternative HTML version deleted]]

[R] Error message in gamm. Problem with temporal correlation structure

HELLO ALL, I AM GETTING AN ERROR MESSAGE WHEN TRYING TO RUN A GAMM MODEL LIKE THE ONE BELOW. I AM USING R VERSION 2.14.1 (2011-12-22) AND MGCV 1.7-12. M1 -gamm(DepVar ~ Treatment + s(Year, by =Treatment), random=list(Block=~1), na.action=na.omit, data = mydata, correlation =

Re: [R] Error message in gamm. Problem with temporal correlation structure

list is DepVar=H, Treatment=Tillage] best, Simon On 17/02/12 09:35, HERNANDEZ PLAZA, MARIA EVA wrote: HELLO ALL, I AM GETTING AN ERROR MESSAGE WHEN TRYING TO RUN A GAMM MODEL LIKE THE ONE BELOW. I AM USING R VERSION 2.14.1 (2011-12-22) AND MGCV 1.7-12. M1-gamm(DepVar ~ Treatment

[R] help with e+01 number abbreviations

Dear List, I will appreciate any advice regarding how to convert the following numbers [I got in return by taxondive()] in numeric integers without the e.g. 6.4836e+01 abbreviations. Thank you very much in advance, Gian taxa_dive Species Delta Delta* Lambda+

Re: [R] R-help Digest, Vol 108, Issue 16

Hi Petr, You advice options(scipen=20) gave me the expected result and fix the problem. Thanks a lot! Gian [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read

  1   2   >