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 έγραψε:

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

e = "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

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

for Component, 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

[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 Je

[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 Mau

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

ackexchange.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 > >

[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 made

[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: plot.roc(daten$Status,daten$SUVmax,percent=TRUE,ci=TRUE,print.auc=TRUE,main="ROC-Kurve

[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 a

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

ASB, I got 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

[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 s

[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 έγραψε: On 20/11/2017 7:09 AM, Maria Lathouri wrote: > Dear all, > > I am really sorry for this. I have attached the s

[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 έγραψε: On 20/11/2017 6:38 AM, Maria Lathouri via R-help wrote: >  Dear all

[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

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

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

ed", "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, ...) {

[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 align

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

confidence intervals.  Kind regards,Maria Στις 4:26 μ.μ. Δευτέρα, 12 Ιουνίου 2017, ο/η Duncan Mackay έγραψε: 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) spring <- dget(file = &q

[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 mo

[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 έγραψε

[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 έγραψε: Dear all, Sorry about that. My mistake. Here is an example of my data > head(aa)       SDATE  va

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

; 2012-12-10 0.484 I hope this helps. Maria Στις 10:00 μ.μ. Δευτέρα, 27 Φεβρουαρίου 2017, ο/η Jeff Newmiller έγραψε: While humorous, the term "Mexican Wall" is unlikely to be clear to the OP. This is a reference to the mailing list anti-virus strategy of cutting out attachments

[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 (with

[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 mo

[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 binomi

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

ame thing.  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

[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 group

[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 the

[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/

[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 NA NA 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

[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

1, 1, 1, 1, 1, 1, 1)), .Names = c("regiao", "estado", "cod_mun", "setor", "cap_int", "idade", "getario", "sexo", "grp_etni", "quest_01", "quest_02", "density", "quest_03&

[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,

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

somewhere 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] 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 understan

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 that

[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] help - hoslem.test

gt; >  load("id3.rda") >  > >  attach(id3) >  > >  > >  #transformando q13 em binária (o/1) >  > >  q131<-ifelse(q13==1,0,ifelse(q13==2,1,ifelse(q13==3,1, >  > ifelse(q13==4,1,ifelse(q13==5,1,NA) >  > id3<-cbind(id3,q131) >  &

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

te...@utoronto.ca] Sendt: 29. april 2015 20:25 Til: Livia Maria Vestergaard Cc: r-help mailing list Emne: Re: [R] LM() and time in R Have a look at the help page for the function ts() (type ?ts at the R prompt). Other than that, you haven't provided nearly enough information for us to diagnos

Re: [R] help - hoslem.test

inomial(link = "logit"), data=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 Cole

[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 and

[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 time

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 ha scritto: > > > On 30 Mar 2015, at 17:11 , Gian Maria Niccolò Benucci < > gian.benu...@gmail.com > wrote: > > > > Dear R-usrs, > > > > I am

[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_fraction

[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 4

[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 function"cox.zph" to check the proportional hazards. Due to the nature of the time-depende

[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 ha

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

..@r-project.org] > On Behalf 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

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

much more that I can say. > Because I'm out of town and will be traveling tomorrow, I'm unlikely to be > able to respond again for several days. > > Best, > John > > On Wed, 9 Jul 2014 17:54:56 -0500 > Maria Judith Carmona H wrote: > > Dear John, > > >

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

ning message: In sqrt (wmd): NaNs produced You see I get the same error. Best regards, Judith On Wed, Jul 9, 2014 at 5:30 PM, John Fox wrote: > Dear Maria Judith Carmona Higuita, > > Since you didn't include enough information (such as your access to your > data) to reproduce the er

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

uot;std"), digits = 2) The error is: 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 -

[R] grouping explanatory variables into "sets" for GLMM

t;? I would like to build 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 i

Re: [R] FrailtyHL package does not work

ork with large data sets. Terry 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 ma

[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 da

[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

gnore that. Do I need to join these parts into a single polygon? 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 d

[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

ce or earlier problems? I have no idea if I am doing anything wrong 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. __

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 wrote: > > The easiest way to supply data is to use the dput() function. Example > with your file named "testfile": > dput(testfile) > Then copy

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... Any

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 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/

[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 mat

[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 https

[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 load

[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 d

[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 constr

Re: [R] polycor package

ress your >> question to the author of that package (or use an appropriate title >> line in your r-help posting). >> >> Best, John >> >> --**-- John Fox Sen. >> William McMaster Prof. of Social Statistics Depa

[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 mess

Re: [R] find the Best-ticker

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

[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, S&P 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 und

[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] "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 (cat

Re: [R] random sampling inside a dataset

n when 'replace = 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 wrote: > On 12-09-13 7:43 AM, Gian Maria Niccolò Benucci wrote: > >&g

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

8800 (ext. 5469) Tel movil: 675 880 297 Fax: 925 268844 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 model

Re: [R] random sampling inside a dataset

> > Anna Freni Sterrantino > Department of Statistics > University of Bologna, Italy > via Belle Arti 41, 40124 BO. > -- > *Da:* Gian Maria Niccolò Benucci > *A:* r-help@r-project.org > *Inviato:* Giovedì 13 Settembre 2012 10:42 > *Oggetto:* [R] ran

[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]] __ R-help@r-pro

Re: [R] INSTRUMENTAL VARIABLES WITH BINARY OUTCOMES

This is the named article: http://ije.oxfordjournals.org/content/37/5/1161.long maybe it can help you to help me... :-( -- View this message in context: http://r.789695.n4.nabble.com/INSTRUMENTAL-VARIABLES-WITH-BINARY-OUTCOMES-tp4642361p4642363.html Sent from the R help mailing list archive

[R] anova for two mixture model objects

del1, model2) : The models 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

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) : Lapa

[R] $ operator is invalid for atomic vectors

egards, 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, m

[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]] ___

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

rrelation = corARMA(form =~ > Year|Tillage/Block, p = 1, q = 0)) > > [I've assumed that the mapping between this message and the data you > sent off list is DepVar=H, Treatment=Tillage] > > best, > Simon > > > On 17/02/12 09:35, HERNANDEZ PLAZA, MARIA EVA wrote

[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 = corARMA(f

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 th

[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+ De

[R] interaction plot and x axis

6.59, 6.77, 6.89, 7.35 but it don't show the ticks in each value. However I would like to define the x values e.g. 5.0, 5.5, 6.0, 6.5, 7.0. I used the xlim but it says error as well >axis(1, at=c(5.0, 5.5, 6.0, 6.5, 7.0)) but still there is a problem. Does anyone know how I can defi

[R] urgent Help needed

Dear Sir, Please help me to open my rda files. It's urgent to see the content and even if I looked for solutions i didn't find one. When I try to load my .rda file I receive this message: "R> load(file = "SpecPhantomExp2WaterSupp Set2.rda") Error: bad restore file magic number (file may be corrup

[R] Specifying grouping factor for augPred in nlme

Dear list, I have a seemingly simple problem with plotting predictions from nlme for which I have not been able to find an answer in the archives nor do I see any mention of this in Pinheiro & Bates 2000. I am using the following model of chick growth for seven different species, with season, nest

[R] Question with RExcel

Dear all, I’m doing some simulation studies in order to compare the estimates (and estimated standard deviations) from the ARMA(2,1) Model with an estimator that I’ve constructed. For carrying out the simulations I created a VBA project within Excel. Now, I’m using the RExcel tool for running th

Re: [R] fitting a distribution to a ecdf plot

Hallo yes I tried it as well and it works; Thank you a lot Maria From: Dennis Murphy [djmu...@gmail.com] Sent: 15 March 2011 21:36 To: Lathouri, Maria Cc: r-help@r-project.org Subject: Re: [R] fitting a distribution to a ecdf plot Hi: The fitdistrplus package

  1   2   >