Re: [R] netcdf data precision or least significant digit

2016-07-07 Thread Ismail SEZEN
Thank you very much Roy and Jeff for your help. I contacted dear David, the author of ncdf4 package, at first and he lead me here. Sorry for any inconvenience. I wish you success in your work. Regards, On Fri, Jul 8, 2016, 04:28 Roy Mendelssohn - NOAA Federal < roy.mendelss...@noaa.gov> wrote:

Re: [R] netcdf data precision or least significant digit

2016-07-07 Thread Roy Mendelssohn - NOAA Federal
I have moved this over to the netcdf-group mail list, which I think is the more appropriate place at this point. You are copied, and hopefully someone from ESRL will see it and provide the proper response. HTH, -Roy > On Jul 7, 2016, at 6:02 PM, Ismail SEZEN wrote: > > Thank you Roy. If I u

Re: [R] netcdf data precision or least significant digit

2016-07-07 Thread Ismail SEZEN
Thank you Roy. If I use "round(uwind, digits = 2)”, all data will have 2 decimal places after decimal point. It’s ok. But How do you know you should round the number to 2 decimal digits? According to definitions of precision and least_significant_digit, should I round to 2 decimal digits or 1 de

Re: [R] netcdf data precision or least significant digit

2016-07-07 Thread Roy Mendelssohn - NOAA Federal
After looking at the file, doing an extract say into the variable uwind, if I do: str(uwind) I see what I expect, but if I just do: uwind I see what you are seeing. Try: uwindnew <- round(uwind, digits = 2) and see if that gives you the results you would expect. HTH, -Roy > On Jul

Re: [R] netcdf data precision or least significant digit

2016-07-07 Thread Ismail SEZEN
Thank you Roy. I use NCEP/NCAR Reanalysis 2 data [1]. More precisely, u-wind data of the year 2015 [2]. I am also pretty sure that the variables like scale_factor or add_offset should be precise like 0.01 or 187.65 but somehow (I hope this is not an issue originated by me) they are not, includ

Re: [R] r code for multilevel latent class analysis

2016-07-07 Thread David Winsemius
> On Jul 7, 2016, at 3:36 PM, Jim Lemon wrote: > > Hi Cristina, > Try this: > > names(mydata) > > It may be NULL or "ppitrst" may be absent. I've already suggested to Christina that she make sure the variables are spelled correctly and she reports they are all present in her dataset. So I t

Re: [R] netcdf data precision or least significant digit

2016-07-07 Thread Roy Mendelssohn - NOAA Federal
Hi Ismail: Can you point me to a particular netcdf file you are working with. I would like to play with it for awhile. I am pretty certain the scale factor is 0.01 and what you are seeing in rounding error (or mor precisely I should say problems with representations of floating point numbers)

Re: [R] netcdf data precision or least significant digit

2016-07-07 Thread Ismail SEZEN
Thank you very much Jeff. I think I’m too far to be able to explain myself. Perhaps, this is the wrong list for this question but I sent it in hope there is someone has deep understanding of netcdf data and use R. Let me tell the story simpler. Assume that you read a numeric vector of data from

Re: [R] r code for multilevel latent class analysis

2016-07-07 Thread Jim Lemon
Hi Cristina, Try this: names(mydata) It may be NULL or "ppitrst" may be absent. Jim On Thu, Jul 7, 2016 at 8:26 PM, Cristina Cametti wrote: > Dear all, > > I am not able to find a reliable r code to run a multilevel latent class > model. Indeed, I have to analyze how social trust (three vari

Re: [R] netcdf data precision or least significant digit

2016-07-07 Thread Jeff Newmiller
Correction: ?options (not par) -- Sent from my phone. Please excuse my brevity. On July 7, 2016 3:26:06 PM PDT, Jeff Newmiller wrote: >Same as with any floating point numeric computation environment... you >don't. There is always uncertainty in any floating point number... it >is just larger in

Re: [R] netcdf data precision or least significant digit

2016-07-07 Thread Jeff Newmiller
Same as with any floating point numeric computation environment... you don't. There is always uncertainty in any floating point number... it is just larger in this data than you might be used to. Once you get to the stage where you want to output values, read up on ?round ?par (digits) and don

[R] Fixed Effects in lme function

2016-07-07 Thread li li
Dear all, For the data below, I would like to fit a model with common random slope and common random intercept as shown below. I am interested in obtaining separate fixed effect estimates (intercept and slope and corresponding hypothesis test) for each method. Instead of performing the analys

[R] Revolutions blog: June 2016 roundup

2016-07-07 Thread David Smith via R-help
Since 2008, Microsoft (formerly Revolution Analytics) staff and guests have written about R every weekday at the Revolutions blog: http://blog.revolutionanalytics.com and every month I post a summary of articles from the previous month of particular interest to readers of r-help. And in case you

[R] Fwd: Re: lmer causes R session to terminate

2016-07-07 Thread Bert Gunter
I failed to c.c. the list. And it should be r-sig-mixed-models list. Bert -- Forwarded message -- From: "Bert Gunter" Date: Jul 7, 2016 2:06 PM Subject: Re: [R] lmer causes R session to terminate To: "David Kikuchi" Cc: Off the top, I would guess that you'd need 10 or more times

[R] lmer causes R session to terminate

2016-07-07 Thread David Kikuchi
Hi, I am working with a large dataset of neotropical birds, and am trying to partition the variance in log(body mass) within different taxonomic levels. To better explain what I mean, the taxonomic levels are species, genus, family, and order. Species are within genera, genera are within famillies

Re: [R] R code for multilevel latent class analysis

2016-07-07 Thread David Winsemius
> On Jul 7, 2016, at 4:31 AM, Yahoo Mail via R-help > wrote: > > Dear all, > > I am not able to find a reliable r code to run a multilevel latent class > model. Indeed, I have to analyze how social trust (three variables form the > ESS survey) might vary between countries (21 countries in my

Re: [R] Bessel function dll file use in VBA

2016-07-07 Thread Jeff Newmiller
This is not a VBA support forum. You need to be studying VBA linkage requirements and gcc linkage conventions, and neither of these subject areas are on topic in R-help. -- Sent from my phone. Please excuse my brevity. On July 7, 2016 5:27:02 AM PDT, "Mehta, Gaurang" wrote: >Hi Team, >I am tr

Re: [R] Performing Principal Cluster Analysis, k-means clustering etc. with PDB/DCD trajectory

2016-07-07 Thread David L Carlson
A Google search turns up a possible place to start. Trajectory Analysis with Bio3D Bio3D1 is an R package that provides interactive tools for the analysis of bimolecular structure, sequence and simulation data. The aim of this document, termed a vignette2 in R parlance, is to provide a brief tas

Re: [R] r code for multilevel latent class analysis

2016-07-07 Thread Anthony Damico
start at https://github.com/ajdamico/asdfree/blob/master/European%20Social%20Survey/structural%20equation%20modeling%20examples.R maybe? On Thu, Jul 7, 2016 at 6:26 AM, Cristina Cametti wrote: > Dear all, > > I am not able to find a reliable r code to run a multilevel latent class > model. Indee

[R] netcdf data precision or least significant digit

2016-07-07 Thread Ismail SEZEN
Hello, I use ncdf4 and ncdf4.helpers packages to get wind data from ncep/ncar reanalysis ncetcdf files. But data is in the form of (9.18, 8.78, 7.98, 3.08, -6.818, …). I’m aware of precision and least_significant_digit attributes of ncdf4 object [1]. For uwnd data, precisio

[R] R code for multilevel latent class analysis

2016-07-07 Thread Yahoo Mail via R-help
Dear all, I am not able to find a reliable r code to run a multilevel latent class model. Indeed, I have to analyze how social trust (three variables form the ESS survey) might vary between countries (21 countries in my database). I tried to use the poLCA package but I am not sure if my code is

[R] Bessel function dll file use in VBA

2016-07-07 Thread Mehta, Gaurang
Hi Team, I am trying to use Bessel function dll files in VBA. It would be great if someone can help. Regards, Gaurang Mehta This email is intended for the person or company named and access by anyone else is unauthorised. If you are not the person or company named, please delete this email and

[R] r code for multilevel latent class analysis

2016-07-07 Thread Cristina Cametti
Dear all, I am not able to find a reliable r code to run a multilevel latent class model. Indeed, I have to analyze how social trust (three variables form the ESS survey) might vary between countries (21 countries in my database). I tried to use the poLCA package but I am not sure if my code is

Re: [R] make a right subset!

2016-07-07 Thread Ulrik Stervbo
Hi Elahe, Use subset and the %in% operator: subset(df2, Serial %in% Matched)[["Count "]] Best wishes, Ulrik On Thu, 7 Jul 2016 at 13:37 ch.elahe via R-help wrote: > Hi all, > I have 2 data frames like the following: > the first one df1: > > 'data.frame': 141obs. of 1 variable: > $Ser

Re: [R] make a right subset!

2016-07-07 Thread ruipbarradas
Hello, Maybe something like the following (untested). idx <- Matched %in% df2$Serial MatchedCount <- df2$Count[idx] Hope this helps, Rui Barradas   Citando ch.elahe via R-help : > Hi all, > I have 2 data frames like the following: > the first one df1: > >    'data.frame':  141obs. of 1 varia

Re: [R] make a right subset!

2016-07-07 Thread Gerrit Eichner
Hello, Elahe, look at ?match and check if df2$Count[ match( Matched, df2$Serial)] does what you want/need. Hth -- Gerrit - Dr. Gerrit Eichner Mathematical Institute, Room 212 gerrit.eich...@math.uni-gie

[R] make a right subset!

2016-07-07 Thread ch.elahe via R-help
Hi all, I have 2 data frames like the following: the first one df1: 'data.frame': 141obs. of 1 variable: $SerialNum: int 41006 41013 41044 41046 41067 41166 41202 41262 41274 41290 and the second one df2: 'data.frame': 194 obs. of 2 variables: $Serial: int 41006 41013 41018 41