Skip to content
Back to formatted view

Raw Message

Message-ID: <CAOKbq8gqmC_nQHE=SAuDG34ehpqgaVBJvqc=A5wZfx1XAJV0ng@mail.gmail.com>
Date: 2026-01-25T23:29:07Z
From: Carlos Ortega
Subject: [R-es] aiken duda
In-Reply-To: <CANuuowOzSt9JDD6tNN4ELhgg0iaWP3WaUZseM3WMN78OTk-MVQ@mail.gmail.com>

Hola,

He probado con diferentes formas que me ha ido proporcionando la IA de
Claude y al final, esta es la que funciona.
Al parecer la función del paquete, tiene un bug (según Claude).

#--------------------------
calcular_aiken_ic <- function(datos, min_val = 1, max_val = 7, conf = 0.95)
{
  n <- nrow(datos)

  resultados <- sapply(datos, function(col) {
    S <- sum(col)
    V <- (S - n * min_val) / (n * (max_val - min_val))

    z <- qnorm(1 - (1 - conf)/2)
    se <- sqrt(V * (1 - V) / n)

    lower <- max(0, V - z * se)  # No puede ser menor que 0
    upper <- min(1, V + z * se)  # No puede ser mayor que 1

    c(AikenV = V, lowerCI = lower, upperCI = upper)
  })

  return(t(resultados))
}

# Aplicar
resultados <- calcular_aiken_ic(df[, 2:ncol(df)])
resultados
#--------------------------

> #--------------------------> calcular_aiken_ic <- function(datos, min_val = 1, max_val = 7, conf = 0.95) {+     n <- nrow(datos)+     +     resultados <- sapply(datos, function(col) {+         S <- sum(col)+         V <- (S - n * min_val) / (n * (max_val - min_val))+         +         z <- qnorm(1 - (1 - conf)/2)+         se <- sqrt(V * (1 - V) / n)+         +         lower <- max(0, V - z * se)  # No puede ser menor que 0+         upper <- min(1, V + z * se)  # No puede ser mayor que 1+         +         c(AikenV = V, lowerCI = lower, upperCI = upper)+     })+     +     return(t(resultados))+ }> > # Aplicar> resultados <- calcular_aiken_ic(df[, 2:ncol(df)])> resultados          AikenV   lowerCI upperCI
item1  0.9000000 0.7685216       1
item2  0.9166667 0.7955377       1
item3  0.8833333 0.7426415       1
item4  0.8750000 0.7300587       1
item5  0.8916667 0.7554547       1
item6  0.8916667 0.7554547       1
item7  0.9083333 0.7818709       1
item8  0.8750000 0.7300587       1
item9  0.9000000 0.7685216       1
item10 0.9083333 0.7818709       1
item11 0.9083333 0.7818709       1
item12 0.9000000 0.7685216       1> #--------------------------

>



Gracias,
Carlos Ortega
www.qualityexcellence.es

El dom, 25 ene 2026 a las 22:17, Jose Betancourt Bethencourt (<
betanster en gmail.com>) escribió:

> rm(list = ls())
>
> library(readxl)
> df <- read_excel("D:/Betan/Escritorio/encuestas/df.xlsx")
>
> df
> library(spheredata)
>
>
>
> attach(df)
> aikenV(df)
>
> ikenV lowerCI upperCI
> item1    1.80     NaN     NaN
> item2    1.83     NaN     NaN
> item3    1.77     NaN     NaN
> item4    1.75     NaN     NaN
> item5    1.78     NaN     NaN
> item6    1.78     NaN     NaN
> item7    1.82     NaN     NaN
> item8    1.75     NaN     NaN
> item9    1.80     NaN     NaN
> item10   1.82     NaN     NaN
> item11   1.82     NaN     NaN
> item12   1.80     NaN     NaNWarning message:
> In sqrt(4 * n * k * V * (1 - V) + (z^2)) : NaNs produced
>
> hay algún error que no veo
> saludos
>
> El dom, 25 ene 2026 a las 13:31, Carlos Ortega (<cof en qualityexcellence.es>)
> escribió:
>
>> Hola,
>>
>> Así, funciona:
>>
>> #----------------------------------------
>> library(spheredata)
>> library(readxl)
>>
>> aikenV(df[, 2:ncol(df)])
>> #----------------------------------------
>>
>> El problema que tenías era que la primera columna "Personas" que son los
>> individuos, no tiene que entrar en el análisis, solo las columnas de las
>> mediciones.
>>
>> Gracias,
>> Carlos Ortega
>> www.qualityexcellence.es
>>
>>
>> El dom, 25 ene 2026 a las 19:21, Proyecto R-UCA (<r-uca en uca.es>)
>> escribió:
>>
>> > Buenas Jose:
>> >
>> > ¿Cómo estás cargando el fichero de datos? Sospecho que el problema es
>> ese.
>> >
>> > ¿De qué paquete estás usando la función aikenV?
>> >
>> > Un saludo.
>> >
>> > El dom, 25-01-2026 a las 09:05 -0500, Jose Betancourt Bethencourt
>> escribió:
>> > > estimados
>> > > al tratar de calcularlo
>> > > aikenV(df)
>> > >
>> > > me da el error
>> > >
>> > >
>> > >
>> > > Error in colSums(data) : 'x' must be numeric
>> > >
>> > > adjunto excel
>> > >
>> > > aprciaria ayuda
>> > > saludos
>> > > jose
>> > > --
>> > > Dr. Jose A. Betancourt Bethencourt
>> > > Universidad de Ciencias Medicas Carlos j. Finlay
>> > > _______________________________________________
>> > > R-help-es mailing list
>> > > R-help-es en r-project.org
>> > >
>> >
>> https://urldefense.com/v3/__https://stat.ethz.ch/mailman/listinfo/r-help-es__;!!D9dNQwwGXtA!XsqN6H5xuanagnHR3-vU-dpF24ARnSKsScgtqztxZKYeQ9ObzrpDH8THhhgQ7Gazuw9msTKXtqNGJ9Y$
>> > >
>> >
>> > _______________________________________________
>> > R-help-es mailing list
>> > R-help-es en r-project.org
>> > https://stat.ethz.ch/mailman/listinfo/r-help-es
>> >
>>
>>
>> --
>> Saludos,
>> Carlos Ortega
>> www.qualityexcellence.es
>>
>>         [[alternative HTML version deleted]]
>>
>> _______________________________________________
>> R-help-es mailing list
>> R-help-es en r-project.org
>> https://stat.ethz.ch/mailman/listinfo/r-help-es
>>
>
>
> --
> Dr. Jose A. Betancourt Bethencourt
> Universidad de Ciencias Medicas Carlos j. Finlay
>


-- 
Saludos,
Carlos Ortega
www.qualityexcellence.es

	[[alternative HTML version deleted]]